VOLUME 8 JANUARY 1959 


ADVANCES 
IN PHYSICS 


A QUARTERLY SUPPLEMENT 
of the 
PHILOSOPHICAL MAGAZINE 


EDITOR 
PROFESSOR N. F. MOTT, M.A., D.Sc., F.R.S. 


EDITORIAL BOARD = 
SIR LAWRENCE BRAGG, 0.B.E., M.C., M.A., D.Sc. F.R.S ~ 


_ SIR GEORGE THOMSON, M.A., D.Sc, F.RSS* > 
PROFESSOR A. M. TYNDALL, C.B.E., D.Sc., F.RS. ©) = 


wee 


C71 
‘= 


NUMBER 29 


PRICE per part £1 
PRICE per annum $3 15s. 0d. post free 


PRINTED AND PUBLISHED BY TAYLOR & FRANCIS LTD 
RED LION COURT, FLEET ST., LONDON E.C.4 


A Wlistory of Mathematics 


From antiquity to the early nineteenth century 


By J. ¥. SCOTT, B.A., DiSe;) PhD: 
Vice-Principal of St. Mary’s College, Strawberry Hill, Twickenham, 
Middlesex 


Author of The Scientific Work of René Descartes (1596-1650), 
Mathematical Work of John Wallis, D.D., F.R.S. (1616-1703), and 
other works 


Contents: Mathematics in Antiquity—Greek Mathematics—The Inven- 
tion of Trigonometry—Decline of Alexandrian Science and the Revival in 
Europe—Mathematics in the Orient—Progress of Mathematics during 
the Renaissance—New Methods in Geometry—The Rise of Mechanics— 
the Invention of Decimal Fractions and of Logarithms—Newton and the 
Calculus—Taylor and Maclaurin, the Bernoullis and Euler, Related 
Advances—The Calculus of Variations, Probability, Projective Geometry, 
Non-Euclidean Geometry—Theory of Numbers—Lagrange, Legendre, 
Laplace, Gauss. 


This volume is intended primarily to help students who desire to have a 
knowledge of the development of the subject but who have too little leisure 
to consult works and documents. The author has availed himself of the 
tacilities afforded by the Royal Society and other learned Societies to re- 
produce extracts from manuscripts and many scarce works. 


““...it has been written with clarity and balance, and the excellent printing helps to make 
it a pleasure to read.”’—The Times Educational Supplement, 21 March 1958. 


“,. his (Dr. Scott’s) wide knowledge of the material, his careful description of methods 


combine to provide an account which at times gives a sense of the excitement of discovery.’’— 
Nature, 26 July 1958. 

“This book will be useful to general readers and as a text-book for students, as well as to 
scholars who are already familiar with the subject.’’—Proceedings of the Physical Society, 
September 1958. 

“The printers and publishers are to be congratulated upon having produced such an 


attractive volume,.... We feel sure the book will be received with delight by all those 
interested in mathematical histories.”"—BEAMA Journal, August 1958. z 


93 x 63 Price 3 guineas plus postage 


Printed and Published by 
TAYLOR & FRANCIS LTD 
RED LION COURT, FLEET STREET, LONDON, E.C.4 


PPL PLP LPP Pd Pegg PP LP ict OOo 


INTERNATIONAL JOURNAL OF 


RADIATION BIOLOGY 
and related studies in Physics, Chemistry and Medicine 


Editor 
W. M. DALE, M.D., D.Sc. 
Dept. of Biochemistry, Christie Hospital & Holt Radium Institute, Withington, Manchester, 20 


Assistant Editor 
JOHN WAKEFIELD, B.A., F.Z.S. 


Editorial Board 


Charlotte Auerbach, U.K.; J. H. Baxendale, U.K.; D. W. van Bekkum, Holland; E. H. Betz, 
Belgium; H. A. S. van den Brenk, Australia; A. A. Buzzatti-Traverso, Italy; J. Coursaget, 
France; F. Devik, Norway; + J. Firket, Belgium; M. Haissinsky, France; A. Hollaender, 
U.S.A.; B. F. Kaufmann, U.S.A.; G. Klein, Sweden; H. Langendorff, Germany; M. Magat, 
France; + N. Miller, U.K.; M. Nakaidzumi, Japan; E. Paterson, U.K.; H. Ulrich, Switzerland; 
M. Westergaard, Denmark. 


Contents of January, 1959 


The Mechanism of Radiation Protection by Histamine and other Biological Amines. By C. van 
der Meer and D. W. van Bekkum, Medical Biological Laboratory of the National 
Defence Research Council T.N.O., Lange Kleiweg 139 Rijswijk (Z.H.), The Netherlands 

5-Hydroxytryptamine as a Radiation Protective Substance in Animals. By H. Langendorff, 
H.-J. Melching and H.-A. Ladner, Radiology Institute, University of Freiburg, Baden, 
Germany 

The DNA Content of the Small Intestine as a Quantitative Measure of Damage and Recovery 
after Whole Body Irradiation. By R. H. Mole and Diana M. Temple, M.R.C. Radio- 
biological Research Unit, A.E.R.E., Harwell, Berks. 

DNA Synthesis in the Thoracic-Duct Lymphocytes of Rats during Recovery from Sub-Lethal 
Irradiation. By E. H. Cooper, The Department of Biochemistry, The University of 
Oxford 

Hexokinase, Aldolase and ATP-Creatine Transphosphorylase in X-Irradiated Rats. By 
M. B. Sahasrabudhe, M. K. Nerurkar and A. J. Baxi, Biology Division, Atomic Energy 
Establishment, Trombay Indian Cancer Research Centre, Parel, Bombay, India and 
D. K. Mahajan, Human Variation Unit, Indian Cancer Research Centre, Parel, 
Bombay, India 3 

On the Relative Role of Mutation and Recombination in Responses to Selection for Polygenic 
Traits in Irradiated Populations of D. Melanogaster. By R. E. Scossiroli and 
S. Scossiroli, Istituto di Genetica, Universita, Pavia, Italy 

An Attempt to Assess the Activity of Tumour Cells by the Study of Nucleolar Morphology, 
Mitotic Index and Phosphorus with **P. By S. Mitra, P. De, R. Chatterjee, 
K. Bhattacharya, A. Bose and H. Dutta, Departments of Cell Research and Biophysics, 
Chittaranjan National Cancer Research Centre, Calcutta 

The Inactivation of SH Enzymes by X-Rays. By Rolf Lange, Alexander Pihl and Lorentz 
Eldjarn, Norsk Hydro’s Institute for Cancer Research, The Norwegian Radium 
Hospital, Oslo, Norway aes 

X-Ray-Induced Abnormal Differentiation of the Epithelium of the Small Intestine in the Mouse. 
By Richard A. McGrath and C. C. Congdon, Biology Division, Oak Ridge National 
Laboratory, Oak Ridge, Tennessee 


Price per part £1 Os. Od. plus postage 
Subscription per volume (4 issues) £3 15s. 0d. post free, payable in advance 


Printed and Published by 


TAYLOR & FRANCIS LTD 
RED LION COURT, FLEET STREET, LONDON, E.C.4 


Orders originating in U.S.A. and Canada should be sent to the 
Academic Press Inc., 111 Fifth Avenue, New York 3, N.Y., U.S.A, 


RRP PPP PPP PPP PPP PL PPP PPL PLP PP PIPPI PLA AAA AAA AAA AL 


A new International Journal 


Molecular Physics 


Editor: H. C. LONGUET-HIGGINS, F.R.S. 
Associate Editor: J. H. VAN DER WAALS 


Editorial Board: 

J. Bjerrum, Copenhagen; G.Careri. Padua; C. A. Coulson, Oxford; F, H.C. Crick, Cambridge: 
P. J. W. Debye, Cornell; D. Hadzi, Liubliana; O. Hassel, Oslo; W. Heitler, Ziirich; J. O. 
Hirschfelder, Wisconsin; D. F. Hornig, Princeton; J. A. A. Ketelaar, Amsterdam, J. G. 
Kirkwood, Yale; R. Kronig, Delft; J. W. Linnett, Oxford; A. Liquori, Rome; Dame Kathleen 
Lonsdale, London; P-O. Léwdin, Uppsala; M. Magat, Paris: R.S. Mulliken, Chicago ; 
A. Miinster, Frankfurt; L. J. Oosterhoff, Leiden; L. E. Orgel, Cambridge; J. A. Pople, 
Teddington; 1. Prigogine, Brussels; R. E. Richards, Oxford; J. S. Rowlinson, Manchester; 
G. S. Rushbrooke, Newcastle upon Tyne; L. E. Sutton, Oxford; H. W. Thompson, 
Oxford; B. Vodar, Bellevue, Paris. 


Contents of January, 1959 


Zur Theorie der generalisierten Gesamtheiten. Von A. Miinster, Physikalisch-chemisches 
Institut der Universitat Frankfurt/Main 


Pressure-dependent Partition Functions. By R. A. Sack, Department of Theoretical Chemistry, 
University of Cambridge and British Rayon Research Association, Wythenshawe, 
Manchester 


On Mayer’s Ionic Solution Theory. By Harold L. Friedman, Faculté des Sciences, Université 
Libre de Bruxelles 


Rotational Effect on the Therm +ynamic Properties of Isotopes. By Agnessa Babloyantz, 
Faculté des Sciences de |’Université Libre de Bruxelles. Belgique 


The Structure of Shock Fronts in Various Gases. By W. H. Andersen and D. F. Hornig. 
Metcalf Research Laboratory, Brown University, Providence, Rhode Island 


Relationship between Proton Shielding Constants and Electric Dipole Moments in the 
Hydrogen Halides. By H.F.Hameka, Department of Chemistry, Carnegie Institute of 
Technology, Pittsburgh, Pennsylvania 


Chemical Shifts of Thallium Resonance Spectra in Solutions of Thallous Salts. By R. Freeman, 
R. P. H. Gasser, R. E. Richards and D. H. Wheeler, Physical Chemistry Laboratory, 
South Parks Road, Oxford 


Correlations between the Electronic Spectra of Alternant Hydrocarbon Molecules and their 
Mono- and Di-valent Ions. 1. Benzene, Coronene and Triphenylene. By G. J. 
Hoiitink, Chemical Laboratory of the Free University, Amsterdam 


Absorption Spectra of Actinide Compounds. By Chr. Klixbiill Jorgensen, Chemistry 
Department A, Technical University of Denmark, Copenhagen, Denmark 


Price per part £1 5s. 0d. plus postage 
Subscription per volume (4 issues) £4 15s. Od. post free, payable in advance 


Printed and Published by 


TAYLOR & FRANCIS LTD 
RED LION COURT, FLEET STREET, LONDON, E.C.4 


Orders originating in U.S.A. and Canada should be sent to the 
Academic¢ Press Inc., 111 Fifth Avenue, New ORK Osa sSecas 


Journal of Electronics 
and Control 


A Philosophical Magazine Associated Journal 


Editor: 
J. THomson, M.A., D.Sc., M.I.E.E., F.Inst.P. 
Consultant Editor: 
Professor N. F. Mort, F.R.S. 
Editorial Board: 

Professor P. AiGRAIN (France) 
Professor H. B. G. Casimir (Holland) 
J. F. Coaces (U.K.) 

K. G. Emeveus (U.K.) 

D. W. Fry (U.K.) 

Dr. W. KLEIN (Germany) 

Dr. R. KOMPFNER (U.S.A.) 


Contents of January, 1959 


Electronics Section 

A Three-Frequency Electron Beam Parametric Amplifier and Frequency Converter. By 
W. H. Louisell, Bell Telephone Laboratories Incorporated, Murray Hill, New Jersey 

The Electrical Conductivity and Hall Coefficient of Bismuth Telluride. By B. Yates, Com- 
munication from the Staff of the Research Laboratories of The General Electric 
Company Limited, Wembley, England 

The Effect of the Inclination of the Focusing Electrodes on Electron Beam Formation. By 
R. J. Lomax, Peterhouse, Cambridge 

On Space Charge Waves. By D. H. Trevena, Ferranti Ltd., Wythenshawe, Manchester 

Discharge Modes using Thermionic Cathodes. By R. B. Cairns and G. C. McCullagh, 
Physics Department, Queen’s University, Belfast 

A Comparison of the Theory of Impact [Ionization with Measurements on Silicon 
p—n Junctions. By F. W. G. Rose, Research Laboratory, British Thomson-Houston 
Co., Ltd., Rugby 

Control Section 


Magnetization Reversal by Rotation and Wall Motion in Thin Films of Nickel-Iron Alloys. 
By E. M. Bradley and M. Prutton, Electronics Research Laboratory, British Tabulating 


Machine Company Ltd., Stevenage, Herts. 


Price per part £1 5s. plus postage 
Price per volume £7 post free, payabie in advance 
6 monthly issues per volume 


Printed and Published by 


TAYLOR & FRANCIS. LTD 
RED LION COURT, FLEET STREET, LONDON, E.C.4 


Orders originating in U.S.A. and Canada should be sent eS the 
Academic Press Inc., 111 Fifth Avenue, New York, 3, N.Y., A, 


A new International Journal 


ERGONOMICS 


HUMAN FACTORS IN WORK, MACHINE CONTROL 
AND EQUIPMENT DESIGN 


General Editor 
A. T. WELFORD 
University of Cambridge, Psychological Laboratory, Downing Place, Cambridge 


Associate Editor (Anatomy and Physiology) Assistant Editor 
W. F. FLloyD Miss H. M. CLAY 
Editorial Board 


H. Bastenier, Belgium; R. B. Bromiley, Canada ; R. Bonnardel, Bernard Metz, France; E. A. 
Miiller, Germany; M. G. Bennett, W. E. Hick, Sir Charles Lovatt Evans, L. G. Norman, 
Great Britain; F. H. Bonjer, Netherlands; S. P. M. Forssman, Sweden; E. Grandjean, 
Switzerland; H.S. Belding, P. M. Fitts, U.S.A. 


Contents of February, 1959 
Symposium on Training 

Introductory Address. By P. H. St. J. Wilson, C.B., C.B.E., Chief industrial Commissioner, 
Ministry of Labour and National Service : 

Physiological Bases of Training. By Albert Hemingway, Department of Physiology, School 
of Medicine, University of Leeds 

Training Operatives in Industry. By W. Douglas Seymour, Clent, Worcestershire : 

The Training of Shoe-Machinists. By W. T. Singleton, The British Boot, Shoe and Allied 
Trades Research Association, Kettering 

A Theory of the Acquisition of Speed-Skill. By E. R. F. W. Crossman, Department of 
Psychology, Reading University, England 

Practice and Knack : Some Comments on Learning and Training in Industry. By Isabel 
Blain, National Institute of Industrial Psychology 

The Operator as a Self-Regulating System: A Factory Experiment. By S. David M. King, 
Organization & Training (Consultants) Limited, 19 Bolton Street, London, W.1 

Note on Current Trends in Literature on Training. By Ezra V. Saul, Institute for Applied 
Experimental Psychology, Tufts University, Medford, Massachusetts 

Circuit Training. By G. T. Adamson, Department of Physical Education, The University, 
Leeds 

The Limitations of a ‘ Procrustean ’ Approach to the Optimization of Man-Machine Systems. 
By F. V. Taylor and W. D. Garvey, U.S. Naval Research Laboratory, Washington, D.C. 

Physical Training in Relation to the Energy Expenditure of Walking and to Factors Controlling 
Respiration During Exercise. By J. E. Cotes and F. Meade, Medical Research Council, 
Pheumoconiosis Research Unit, Llandough Hospital, Penarth, Glamorgan, Wales 

Problems of Training of the Cardiovascular System. By M. J. Karvonen, Institute of 
Occupational Health, Helsinki, Finland 

Training Muscle Strength. By E. A. Miiller, Max-Planck-Institut fiir Arbeitsphysiologie, 
Dortmund, Germany 

Abstracts of other Papers Delivered at the Symposium 


Price £1 5s. 0d. per part plus postage 
Subscription price per volume £4 15s. Od. post free, payable in advance 


Printed and Published by 
TAYLOR & FRANCIS LTD 
RED LION COURT, FLEET STREET, LONDON, E.C.4 


Orders originating in U.S.A. and Canada should be sent to the 
Academic Press Inc., 111 Fifth Avenue, New York, 3, N.Y., U.S.A. 


he Philosophical Magazine 


First Published in 1798 


Editor : 
PROFESSOR N. F. Mott, M.A., D.Sc., F.R.S. 


Editorial Board : 
Sik LawRENCE BRAGG, O.B.E., M.C., M.A., D.Sc., F.R.S. 
Str GeorGE THomson, M.A., D.Sc., F.R.S. 
PROFESSOR A. M. TYNDALL, C.B.E., D.Sc., F.R.S. 


Contents of February, 1959 


The Effects of Anisotropic Relaxation Times on the Hall Coefficients of some Dilute Alloys 
of Silver. By J. R. A. Cooper and S. Raimes, Department of Mathematics, Imperial 
College, London 


The Resistance of Diamond to Abrasion. By Eileen M. Wilks and J. Wilks, Clarendon 
Laboratory, Oxford 


Some Properties of Crystals of Silver Chloride containing Traces of Copper Chlorides. By 
ae S: eae and J. W. Mitchell, H. H. Wills Physical Laboratory, University of 
Tisto 
On the Knight Shift of Alloys. By A. Blandin, E. Daniel, and J. Friedel, Physique des Solides, 
60 Boulevard Saint Michel, Paris VI 


Measurements of the Optical Constants of Mercury and Mercury—Indium Amalgams in the 
Spectral Region 4000 to 17000cm~. By J. N. Hodgson, Physics Department, 
University College of North Staffordshire 

The Effect of Grain Size on the Tensile Properties of High-purity Molybdenum at Room 
Temperature. By A. A. Johnson, Physical Metallurgy Department, Imperial College, 
London, S.W.7 

Observations on Single Crystals of an Isotactic Polyolefin: Morphology and Chain Packing 
in Poly-4-Methyl-Pentene-1. By F. C. Frank, A. Keller and A. O’Connor, H. H. 
Wills Physics Laboratory, University of Bristol 

Nucleon Clusters in the Nuclear Surface. By D. W. Wilkinson, Clarendon Laboratory, 
Oxford 

On the Low Temperature Optical Absorption of CoCl,.6H,.O and CoBr,.6H,0, By 
R. Pappalardo, H. H. Wills Physical Laboratory, University of Bristol and Radiation 
Laboratory, University of Pittsburgh, Pa. 

Intensity and Fine Structure of the Main Optical Absorption Bands in Hydrated Cobaltous 
Salts. By Shoichiro Koide, H. H. Wills Physics Laboratory, University of Bristol 

Adhesion of Evaporated Aluminium Films with Underlayers of Nickel, Cobalt and Chromel. 
By C. Weaver and R. M. Hill, Department of Natural Philosophy, The Royal College 
of Science and Technology, Glasgow 


Price per part £1 5s. plus postage 
Price per annum £13 10s. post free, payable in advance 


Printed and Published by 
TAYLOR & FRANCIS LTD 
LION COURT, FLEET STREET, LONDON, E.C.4 


Physics in Medicine and Biology 


A Taylor & Francis Journal published in association 
with The Hospital Physicists’ Association 


Editor: J. E. Roperts, D.Sc. 
Consultant Editor: Professor N. F. Mort, F.R.S. 


Editorial Board 


R. Bonet-Maury, Paris; J. Dainty, Edinburgh; H. E. Johns, Toronto; W. A. Langmead. 
London; D. A. McDonald, London; J. S. Mitchell, Cambridge; G. J. Neary, Harwell; 
B. Rajewsky, Frankfurt; J. Rotblat, London; S. Rowlands, London; H. P. Schwan, Philadelphia: 
R. Sievert, Stockholm; F. W. Spiers, Leeds; J. F. Tait, Massachusetts; A. J. H. Vendrik, 


Nijmegen. 


Contents of January, 1959 


The Measurement of Exposure Dose for High Energy Radiation with Cavity Ionization 
Chambers. By T. E. Burlin, B.Sc., Physics Department, Hammersmith Hospital, 
London, W.12 


The Energy Response of a Victoreen Condenser Ionization Chamber. By A. J. D. 
Szilvasi, Ph.D. and G. N. Whyte, Ph.D., National Research Council, Ottawa, Canada 


An Analysis of Clasp Design in Partial Dentures. By J. A. Warr, M.A.(Cantab.), B.Sc. 
(London), B.D.S.(London), L.D.S.R.C.S.Eng., Department of Prosthetics, The 
London Hospital Medical College, Turner Street, London, E.1 


A Study of the External Gamma-ray Field Arising from a Radioactive Source Situated in 
Air and in a Scattering Medium. By H. A. B. Simons, M.A., Ph.D. and S. Chomet, 
M.Sc., Physics Department, The Royal Free Hospital School of Medicine, 8 Hunter 
Street, London, W.C.1 


A Wide Range Direct Reading X-Ray Dosemeter. By S. Stuart, Electronic Instruments 
Ltd., Richmond, Surrey and S. B. Osborn, B.Sc., A.Inst. P., University College Hospital, 
London, W.C.1 


Pulmonary Diffusion of Oxygen. By B. F. Visser, Ph.D. and A. H. J. Maas, M.Sc., Pulmonary 
Function Laboratory, State University Hospital, Utrecht, Netherlands 


The Influence of the Anomalous Viscosity of Blood upon its Oscillatory Flow. By M. G. 
yea M.D., Physiology Department, Medical College of St. Bartholomew’s Hospital, 
ondon 


Subscription price per volume £3 10s. post free, payable in advance 
4 parts per volume—£1 per part plus postage 


Printed and Published by 


TAYLOR & FRANCIS, LTD 


RED LION COURT, FLEET STREET, LONDON, E.C.4 


Orders originating in U.S.A. and Canada should be sent to the 
Academic Press Inc., 111 Fifth Avenue, New York, 3, N.Y., U.S.A. 


INTERNATIONAL JOURNAL 


OPTICA ACTA 


European Journal Journal Europeen Europaische Zeitschrift 
of Optics d’Optique fiir Optik 


English Editor: C. G. WYNNE, 13 Elwill Way, Beckenham, Kent, U.K. 
Rédacteur Francais: A. MARECHAL, 3 Boulevard Pasteur, Paris, France 
Deutscher Herausgeber: G. FRANKE, Laufdorfer Weg 2, Wetzlar, Germany 


Editorial Board 


J. M. Otero y Navascues, Madrid; R. W. Ditchburn, Reading; A. Arnulf, Paris; A. Biot, 
Gand; J. Cabannes, Paris; P. Fleury, Paris; F. Gabler, Vienna; G. Hansen, Oberkochen; 
E. Ingelstam, Stockholm; H. Korte, Braunschweig; W. S. Stiles, London; G. Toraldo di 
Francia, Florence; A. C. S. van Heel, Delft. 


Contents of January, 1959 


N.P.L. Colour-Matching Investigation: Final Report (1958). By W. S. Stiles and J. M. 
Burch, National Physical Laboratory, Teddington, Middlesex 


Optical Constants of Evaporated Films of Zinc Sulphide and Germanium in the Infra-Red. 
By Lennart Huldt and Torsten Staflin, Institute of Optical Research, Royal Institute 
of Technology, Stockholm 


Das Moiré-Gitter als vielseitiges Testobjekt: Photoelektrische Aberrationsmessung. Von A. 
Lohmann, z.Zt. Institutet for optisk forskning, Stockholm 70 


Matrix Formalism in the Theory of Diffusion of Light. By I. KuScer, Institute of Physics, 
University of Ljubljana and M. Ribaric¢, J. Stefan Institute, Ljubljana, Yugoslavia 


Lage und Sichtbarkeit von Keilinterferenzen bei Instrumenteller Beobachtung (Beitrag zur 
Kohiarenztheorie eines Unsymmetrischen Interferenzfalles). Von Rudolf Landwehr, 
Menden/Sauerland, Deutschland 


Essai de Détermination des Trés Faibles Diamétres Apparents en Radioastronomie. Par 
J. Arsac, Observatoire de Meudon 


Price per part £1 (1375 fr.), (DM 11.75), ($2.80), plus postage 
Subscription price per volume (4 parts) £3 15s. (5155 fr.), (DM 44), ($10.50), post free, 
payable in advance 


Printed and Published by 


TAYLOR & FRANCIS LTD 
RED LION COURT, FLEET STREET, LONDON, E.C.4 


A HISTORY OF 
MATHEMATICS 


from antiquity to the 
early nineteenth century 


byg er SCO Liab Ave ScerbnD: 
Vice-Principal of St. Mary’s College, 
Strawberry Hill, Twickenham, 
Middlesex 


Author of The Scientific Work of René 
Descartes (1596-1650), Mathematical 
Work of John Wallis, D.D., F.R.S. 
(1616-1703), and other works. 


CONTENTS: Mathematics in Anti- 
quity—Greek Mathematics—The Inven- 
tion of Trigonometry—Decline of Alex- 
andrian Science and the Revival in 
Europe—Mathematics in the Orient— 
_ Progress of Mathematics during the 
Renaissance—New Methods in Geo- 
metry—The Rise of Mechanics—the 
Invention of Decimal Fractions and of 
Logarithms—Newton and the Calculus 
—Taylor and Maclaurin, the Bernoullis 
and Euler, Related Advances—The 
Calculus of Variations, Probability, 
Projective Geometry, Non-Euclidean 
Geometry—Theory of Numbers— 
Lagrange, Legendre, Laplace, Gauss. 


This volume is intended primarily to 
help students who desire to have a 
knowledge of the development of the 
subject but who have too little leisure 
to consult original works and docu- 
ments. The author has availed himself 
of the facilities afforded by the Royal 
Society and other learned Societies to 
reproduce extracts from manuscripts 
and many scarce works. 


93x6} Price 3 guineas 


Printed and Published by 
TAYLOR & FRANCIS LTD. 
Rep LION CourT, FLEET STREET, 

Lonpon, E.C,4 


THE SCIENTIFIC 
WORK OF 
RENE DESCARTES 
1596-1650 


BY 
J-°Fy SCOTT, 3.4... MSC, 4PH-D. 


With a Foreword by 
H. W. TURNBULL, M.A., F.R.S. 


This book puts the chief mathe- 
matical and physical discoveries 
of Descartes in an accessible form 
and fills an outstanding gap upon 
the shelf devoted to the history 
of philosophy and science. There 
is to be found in this volume the 
considerable contribution that 
Descartes made to the physical 
sciences, which involved much 
accurate work in geometrical 
optics and its bearing upon the 
practical problem of fashioning 
lenses, as also the deeper problems 
of light and sight and colour. The 
careful treatment that Dr. Scott 
has accorded to the work of 
Descartes is very welcome. The 
book is well worth reading and 
will be an asset to all libraries. 
This publication is recommended 
and approved by the Publication 
Fund Committee of the University 
of London. 


212 pp. 7X10 amply illustrated 


PRICE £1-0-0 NET 


Printed and Published by 
TAYLOR & FRANCIS LTD. 
ReD LION CourT, FLEET STREET, 
Lonpon, E.C.4 


CONTENTS 


The Theory of Superconductivity. By C. G. Kuper, Department of 
Natural Philosophy, St. Salvator’s College, St. Andrews 


The Modern Theory of Superconductivity. By I. M. KaauatniKov and 
A. A. Aprikosov, Institute of Physical pbeae Kaluzhskove 
Shosse, 32 Moscow WES Sek: : 


45 


nde iad © oes athe 
Ay ake oe ey 


ee ee 


ADVANCES IN PHYSICS 


A QUARTERLY SUPPLEMENT 
of the 


PHILOSOPHICAL MAGAZINE 


VOLUME 8 January 1959 NUMBER 29 


The Theory of Superconductivity 


By C. G. Kuper. . 
Department of Natural Philosophy, St. Salvator’s College, St. Andrews 


CONTENTS 


§ 1. QuaLiraTIvE Discussion. 
1.1. Introduction. 
1.2. The lattice—electron interaction. 


§ 2. Tue Hamintrontan FoRMULATION. 
2.1. Outline of the Bloch theory. 
2.2. The quasi-free electron model. 
2.3. Perturbation theory. 
2.4, Renormalization of the phonon energy and coupling strength. 


§ 3. THe Grounp State and ENERGY Gap. 
3.1. Fréhlich’s one-dimensional model. 
3.2. The three-dimensional model of Bardeen, Cooper and Schrieffer. 
3.3. Bogoliubov’s ‘ quasi-particle ’ transformation. 
3.4. The role of Coulomb interactions. 
3.5. The energy gap at finite temperatures. 


. ELECTRODYNAMIC PROPERTIES OF SUPERCONDUCTORS, 
4.1. Supercurrents. 
4.2. The Meissner effect. 
4.3. Paramagnetism. 


ez) 
is 


REFERENCES. 


§ 1. QUALITATIVE DIscussION 


1.1. Introduction 


SUPERCONDUCTIVITY has, for many years, been one of the central unsolved 
problems of physics. The phenomenon was discovered by Kamerlingh 
Onnes as long ago as 1911; many metals and alloys, below a critical tempera- 
ture 7’, characteristic of the substance, have no measurable electrical 


P.M.S. B 
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resistance. At a given temperature 7’, below 79, there is a threshold 
magnetic field H,(7'), such that if H>H, the superconducting property 
is destroyed, The third fundamental property, discovered by Meissner 
and Ochsenfeld in 1933, is that a macroscopic superconductor is perfectly 
diamagnetic, i.e. no magnetic field can exist inside the superconductor. 

A very successful phenomenological theory was advanced by F. and ink 
London in 19357, but until 1950 the nature of the phenomenon remained a 
complete enigma. All other macroscopic properties of metals are, at 
least qualitatively, understood; in particular developments of Bloch’s 
(1928) theory of metallic conduction account for the electronic properties 
quite well. Nevertheless, the interaction responsible for superconductivity 
remained obscure. Most theoretical speculations assumed a mechanism 
depending on some of the factors left out of the Bloch scheme, e.g. Coulomb 
interactions or the details of the Brillouin zone structure. 

An important advance was made by Frohlich (1950), who suggested 
that Bloch’s electron—lattice interaction should be responsible for super- 
conductivity. Bloch had calculated matrix elements for emission and 
absorption of lattice quanta by the electrons in a metal and showed that 
(real) processes of this kind are responsible for the normal resistance of 
metals. Fréhlich pointed out that Bloch’s theory can be expressed as a 
field theory, and that it implies the possibility of virtual absorption and 
emission processes. These lead to an electron—electron interaction. 
As a consequence of this suggestion, the energy difference between the 
normal and superconducting states is proportional to M~-1, the inverse 
isotopic mass of the lattice ions. Soon after the prediction, the isotope 
effect was found experimentally, and confirms the above dependence. 
It is therefore quite clear that the electron—lattice interaction of Bloch is 
of basic importance for the understanding of superconductivity. Never- 
theless, subsequent developments have been slow, on account of the great 
mathematical difficulties encountered (Bardeen 1950, 1951, Fréhlich 1952, 
1954, Bardeen and Pines 1955). With the very recent work of Bardeen 
et al. (1957), Bogoliubov (1958a) and Valatin (1958), however, it has 
been shown that the lattice—-electron interaction leads to a ‘ condensation’ , 
and an energy gap in the excitation spectrum. The thermal properties of 
superconductors are accounted for quite well, and also the existence of 
curl-free persistent currents. But there is, as yet, no satisfactory theory 
of the Meissner effect (nor of the anomalous Knight shift). 


ee 
tA parameter A with the dimensions of a length is introduced, a * penetration 
depth ’, and the current is split into two parts, a ‘normal’ component j,, and 


- ‘ s > :, Tye a" * . 
a superconducting component j,. The superconducting current j, obeys 
the following equations in place of Ohm’s law : 


(—47A?/c) curl j, = B, 


(Arr d2/e2) ee (1.1) 
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Before discussing the successes and limitations of the theory in its 
present form, it is useful to bear in mind the following properties of 
superconductors : 


(a) the phenomenological theory of the Londons governs the electro- 
dynamic properties quite wellt. 

(b) the superconducting state is a distinct thermodynamic phase, 
and thermodynamic formulae of the usual type can be applied to the 
normal-superconducting transition, in agreement with experiment. 

(c) the actual energy difference between the normal and_super- 
conducting phases is exceedingly small (~4x 10?ergcm-*) compared 
with the electronic kinetic energy (~ 10'ergcm-*). Despite the smallness 
of the energy difference, many of the macroscopic properties are quite 
drastically altered. Apart from the electrodynamic properties, there 
are large anomalies in, for example, the specific heat and thermal 
conductivity. It is the smallness of this transition energy, compared 
with other relevant energies, that is responsible for many of the difficulties 
in constructing a theory. 

Experimental evidence is rapidly accumulating in favour of an ‘energy 
gap ’, i.e. there appears to be a finite separation between the ground state 
and all single-particle excitations. As a phenomenological model 
(Ginsburg 1946, Koppe 1947, 1950, Goodman 1953) the energy-gap can 
account for the specific heat and thermal conductivity, and for the absorp- 
tion edges found in acoustic and electromagnetic absorption experiments. 
Further experimental evidence for an energy gap is afforded by nuclear 
spin-lattice relaxation time measurements. 

It has been claimed (Bardeen 1956) that the energy gap model, with 
certain additional assumptions, shows also the correct electrodynamical 
properties (leading to Pippard’s generalization of the London equations), 
but this claim has been shown not to be valid (Buckingham 1957). 

The scope of this article will be restricted to an account of the develop- 
ment and present state of the microscopic theory. No attempt will be 
made to describe the phenomenological developments, nor the extent to 
which the phenomenological theories account for the experimental results. 


1.2. The Lattice—-electron Interaction 


Bloch’s (1928) theory of metallic conduction assumes that the valence 
electrons in a metal are nearly independent, and move in the field of the 
metallic ions. The Coulomb electron-electron forces are neglected or 


+ Detailed quantitative discrepancies exist however. To overcome them 
we need an improved phenomenological theory. Several have been advanced, 
by Schafroth and Blatt (1955, 1956), Ginsburg and Landau (1950) and 
Pippard (1953), the latter being probably the most successful. They all 
introduce an additional parameter, a ‘ range of order ’ or ‘ correlation length ’. 
Qualitatively, the spatial extension of the Cooper (1956) quasi-molecules (§ 3.2) 
probably provides the theoretical basis for the correlation length. This aspect. 
of the theory is, however, still insufficiently developed. 
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at best treated by a self-consistent field approximation (though the 
application of plasma theory (Bohm and Pines 1953, Pines 1953, Bardeen 
and Pines 1955) helps to justify this approach). In a perfectly periodic 
field, the single-electron wave functions are of the form 


w(x) =exp(ik.x)ju(x) . . . « a - « (1.2) 


where u(x) is a function with the periodicity of the lattice. The total 
wave function is a Slater determinant, built out of these one-electron 
functions. The electrons behave rather like free electrons, but with an 
effective mass m which may differ from the real mass. 

These ‘quasi-free’ electrons are scattered by departures from 
periodicity: (a) lattice defects of many kinds (impurities, dislocations, 
vacancies, etc.) produce an approximately temperature-independent 
scattering—the ‘residual’ resistance, and (b) the thermal motion of the 
ions leads to a temperature-dependent scattering, the ‘ideal’ resistance, 
which vanishes as 7'->0. The latter mechanism was described by Bloch 
in terms of the emission or absorption by an electron of a ‘phonon’, a 
quantum of excitation of one of the normal modes of vibration of the 
lattice. 

The realization that the interaction of lattice vibrations with electrons 
implies an interaction between electrons was the first essential step, made 
by Frohlich in 1950. Previously Froéhlich et al. (1950) had studied the 
behaviour of an electron in an ionic medium (‘polaron’). They formulated 
the interaction in a field-theoretical way—the electron produces a field of 
polarization of the lattice, and that field in turn reacts back on the electron, 
and leads to a self-energy. Provided the interaction is weak enough, the 
self-energy can be calculated by perturbation theory. The Hamiltonian 
for the polaron is quite similar in form to the Hamiltonian implied by 
Bloch’s theory, except in so far as there are many electrons instead of a 
single one. Fréhlich, therefore, by writing the Bloch interaction in 
Hamiltonian form, showed that it would also lead to an electron self-energy ° 
and an electron-electron interactiont. The latter effect arises from 
processes in which one electron emits a phonon, which another electron 
absorbs. 

In the Bloch model of a metal, the self-energy is by far the larger effect, 
but it is also the less interesting. It arises from the fact that for most 
lattice waves the electronic motion can follow the lattice nearly adiabati- 
cally (Huang 1951, Frohlich 1952). The electrons enhance the inertial 
properties of the lattice, leading to a renormalization of the velocity of sound, 
and the single-particle energy levels suffer a self-energy correction. For 
electrons near the Fermi surface, however, another possibility exists—the 


} Although the analogy is incomplete, a comparison with quantum electro- 
dynamics is illuminating. There the electron—photon interaction leads to 


charge and mass renormalization, and to the Coulomb interaction between 
electrons. 
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lattice displacements can follow the electron density fluctuations adiabati- 
cally instead and, by reacting back on the electrons, produce an interaction 
between electrons. 

One of the great difficulties confronting the development of the theory was 
the problem of isolating the small interaction energy, in the presence of the 
much larger self-energy. Application of second order perturbation theory 
(Fréhlich 1950) indicates the presence of the two effects, but implies that if 
the interaction is strong enough to be of interest, then the lattice is unstable 
(Wentzel 1951). 

In order to isolate the (small) interaction energy and study its properties, 
the (large) self-energy must somehow beremoved. Frohlich (1952) showed 
how this can be done by a unitary transformation of the Hamiltonian. 
The transformation leads to (a) a renormalized sound velocity s (from the 
inertia of the electrons which follow the lattice motion), (b) a renormalized 
coupling constant, (c) a renormalization of single particle energy levels, (d) 
an electron—phonon interaction restricted to transitions on an energy shell, 
and (e) an electron—electron interaction energy, positive when the electron 
energy difference = hkys, negative for energy differences Shkys. The order 
of magnitude of this energy is ~ F',/N,Q ~ F/Q& where Q is the normaliza- 
tion volume, V, the number of electrons per unit volume, F is defined in 
§ 2.2, and is a dimensionless measure of the strength of the lattice—electron 
interaction, and & is the density of levels, per unit volume, at the Fermi 
surface C,. 

If any ‘condensation’ in momentum space is to occur, it has to do so 
despite an increase in the electronic kinetic energy, and will therefore be 
restricted to an energy shell, of thickness ~ hk ys at the Fermisurface. This 
shell will contain ~fkys . QE electrons, and therefore the interaction energy 
will be roughly 


fhkys QS? ~ Fms2.N .Q 


~ FN Q(m[M)y, . . . .. (1.3) 


where y is independent of the isotopic mass M of the ions. The dynamic 
nature of the effect is apparent in the dependence of the energy on the 
isotopic mass, via the sound velocity s. 

The condensation energy is proportional to the square of the threshold 
field at absolute zero, and hence the theory predicts 


ERUVE A ced E sake re eee om re nea 
Since empirically H,(7' = 0) cT'y, the prediction can also be written 
GN COM UES ho SPS ial ea ie Le) 


Both these predictions (1.4) and (1.5) are in quantitative agreement with 
experiment (Maxwell 1950, Reynolds e¢ al. 1951, Allen e¢ al. 1950 a, b, Lock 
etal. 1951). These experiments indeed afford the most direct evidence that 
the electron—phonon interaction is responsible for superconductivity. The 


6 C. G. Kuper on the 


condensation energy (1.3) is however too large, compared with observed 
values, by a factor ~10%. Since hkys~«@, where « is Boltzmann’s 
constant and © is the Debye temperature, (1.3) would predict 7,)~0, 
while usually the values satisfy 10-°0 <7) < 1070. 

The low-lying excitations of a metal are highly degenerate ; consequently 
it is not satisfactory to use the Bloch wave functions as unperturbed states 
in a perturbation theoretical approach. The perturbation theoretical 
calculation of Frohlich (1950) served mainly to bring out the above orders 
of magnitude, and to predict the isotope effectf. 

For a system with an n-fold degenerate level, where n is finite, there is a 
well-known systematic procedure for selecting appropriate linear combina- 
tions of the degenerate wave functions, to lift the degeneracy. But fora 
continuum, there is no such systematic procedure. For further progress, 
clearly, what is needed is a new set of ‘ unperturbed’ wave functions, formed 
by appropriate linear combinations of Bloch functions. In 1953, Casimir 
commented: ‘“‘...I believe that one should try to make a good guess as to 
what the state of affairs is, and then calculate the energy.” Further pro- 
gress indeed awaited the ‘good guess’. 

One of the features which is lacking in the Bloch wave functions and 
which the improved basic functions should exhibit is the cooperative nature 
ofthe phenomenon. The specific heat at 7’= 7, shows a discontinuity, an 
Ehrenfest (1933) ‘‘ transition of second order”. Just as the Weiss (1907) 
theory of ferromagnetism assumes an ‘internal field’, which acts as an 
order parameter and vanishes at the Curie temperature, so we can demand a 
parameter here with a similar role. It is also desirable that the order 
parameter should have an obvious physical significance (as the Weiss field 
certainly has), so as to give as much insight as possible into the behaviour 
of the system. 

Up to about 1954, the sole success of the theory had been to predict the 
isotope effect, while in all other directions it had encountered such severe 
mathematical difficulties that it was doubtful whether the mechanism could 
produce a condensation. It was therefore desirable to show that a con- 
densation of some sort could occur as a result of the electron-phonon 
interaction, even in a simplified model, with unrealistic values of the para- 
meter. Frohlich (1954) studied a one-dimensional model with very strong 
coupling. In a one-dimensional metal the degeneracy is much lower than 
in three dimensions. By making the interaction sufficiently strong, he was 
able to neglect the recoil of the lattice modes when an electron is scattered, 
and to develop a Hartree self-consistent field approximation. The lattice 
mode of wave number wy = 2k, can become strongly excited, and acts as a 
spatially periodic field on the electrons, producing a gap in the single- 
particle energy spectrum (Peierls 1930). The electron density varies with 


ee ee 


} Starting from a knowledge of the experimental discovery of the isotope 
effect, Bardeen (1950) developed another formulation, which he later showed 
(1951) to be equivalent to Fréhlich’s. 
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the same periodicity, producing the periodic lattice distortion in a self- 
consistent way (§3.1). The approximations are valid only if the coupling 
is very strong (which implies 7) ©). 

The resulting energy gap is 


P= ¢, exp (—3/2F), salty is nr) carom 66s)! 
and the condensation energy per atom 
WE OGOX (ee BIE) ty on at fs eal l.'Z) 


where ()=?k,?/2m is the Fermi energy. The isotope effect has been lost, 
because the dynamic properties of the interaction have been ignored. The 
excitation of the lattice mode w,=2k, acts as an internal field, whose 
quantum properties have been neglected. This single mode is so strongly 
excited that the electronic motion follows it adiabatically ; this is just the 
condition we need in order that the interaction energy is as large as possible 
compared to the self-energy. In fact these two energies come out of the 
same order of magnitude. 

The amplitude of the lattice mode acts as an ‘internal field’, and provides 
the sought-after analogy with ferromagnetism at absolute zero. Moreover 
the self-consistent calculation can be extended to finite temperatures, 
giving a temperature-dependent gap width ['(7’) which vanishes at some 
temperature 7’), again in analogy with the Weiss field. The specific heat 
is exponential near 7'’= 0, and has a discontinuity at 7’) (Kuper 1955). 

Although the ‘ Hartree field’ method cannot be generalized, either to 
three dimensions or to weaker interactions, it is important in giving an 
energy gap proportional to exp(—3/2F). This function has an essential 
singularity at / = 0, and cannot arise in any finite order perturbation theory. 
If for some reason the exponential factor should prove to be correct, evenin a 
three-dimensional system with weak coupling, then one can understand the 
smallness of the condensation energy of the superconducting phase com- 
pared with (1.3), since exp (—1/F) is a very rapid function of F’, and, for 
example, is ~ 10-? for F~1/5. Also since F will vary a little under com- 
pression, theverysensitive observed pressure dependence of 7’, would follow. 

The reason the one-dimensional model shows no isotope effect is that the 
interaction has been made so strong that all electrons follow the lattice 
motion adiabatically—the lattice motion is treated quasi-statically. If 
we can restore the dynamic properties of the interaction, we might hope to 
get an energy containing both the factor ms? of the perturbation method and 
the exponential factor of the Hartree field method. 

By the choice of a good set of trial wave functions, Bardeen ef al. (1957) 
achieve just this desired combination. Cooper (1956) had noticed that 
under the superconductive interaction, two electrons near the Fermi 
surface could form a ‘molecule’. If indeed all (or at least many) of the 
electrons could become bound in pairs, and if these pairs did not overlap 
appreciably, then, as Schafroth et al. (1957) had shown, the system would 
exhibit superconducting properties by virtue of the Bose-Einstein con- 
densation of the (Boson) ‘molecules’, the condensed phase being made up 
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of pairs with zero total momentum. Actual Cooper pairs are not highly 
localized; the spread is ~ 10-4cm, so that the appearance of the condensed 
phase does not arise in quite this simple way. The essential features are 
still present in the Bardeen—Cooper-Schrieffer ground state wave function ; 
electrons are all correlated in pairs, with equal and opposite momenta. 

As in the one-dimensional model, the spectrum of elementary excitation 
energieshasagap. There are also collective motions of the whole electronic 
system, with low energies; these have a macroscopic non-zero current 
(again as in the one-dimensional model). The superconducting ground 
state is energetically more favourable than the Bloch ground state for all 
positive F, however small. Bardeen, Cooper and Schrieffer however take 
the screened Coulomb interaction into account in an approximate way ; 
the Coulomb repulsion makes a negative contribution —A to the ‘effec- 
tive’ coupling constant, /’ — A, say, and the criterion for superconductivity 
is that the ‘effective’ coupling constant F — A is positive. 

The condensation energy found by Bardeen, Cooper and Schrieffer is 
(neglecting the Coulomb repulsion) 


W~ms*exp(—2/F),. . . . 5g egg ama 


and should be compared with (1.7) and (1.3). The desired linking of the 
exponential factor with the isotope effect factor ms? has been achieved. 

The derivation of (1.8) by Bardeen, Cooper and Schrieffer requires many 
simplifications and approximations; in particular the lattice has only one 
vibrational frequency, and the matrix elements of the (renormalized) 
Hamiltonian are replaced by a simplified set. In consequence, the results 
require justification in a more rigorous fashion. A transformation of the 
Hamiltonian, proposed by Bogoliubov (1958 a,b) and Valatin (1958) 
provides the necessary justification. This transformation is an extension 
to a system of Fermions of Bogoliubov’s (1947) Boson transformation. It 
introduces a new set of field variables to replace the original electron field. 
The transformation exhibits, in a fairly explicit form, the (+k, + ; —k, 4 ) 
electron pair correlations in the ground state. The new variables can be 
thought of as describing a set of ‘ quasi-particles’. Although the number 
of real electrons is a constant of the motion, the number of quasi-particles 
unfortunately is nott. By a suitable choice of the parameters of the trans- 
formation, the ‘vacuum’ state with no quasi-particles becomes identical 
with the ground state of Bardeen, Cooper and Schrieffer. A simultaneous 
transformation of the phonon field variables exhibits Fréhlich’s 
renormalization (up to second order in the coupling constant). The energy 
spectrum of excited states, calculated from the transformed Hamiltonian, 
exhibits an energy gap . 

Dv nkys xp (UN) og ee eee ea 


+ The non-conservation of quasi-particles is linked with the difficulty in 
getting a clear physical picture of the mechanism of the ‘ internal field”, in 
contrast to Fréhlich’s one-dimensional case, § 3.1. 


Theory of Superconductivity i) 


at absolute zero. (It is interesting to compare (1.8) and (1.9) with (1.6) 
and (1.7)—see § 3.2.) 

At higher temperatures, a self-consistent calculation shows, as in the 
one-dimensional case, that the gap I'(7’) diminishes with increasing tem- 
perature 7’, vanishing at some temperature 7). The second-order transi- 
tion and all other properties which depend mainly on the energy gap and the 
density of states follow immediately, in good agreement with experiment. 
Examples are the specific heat (Zavaritskii 1957, Lynton et al. 1957, 
Corak et al. 1954, Corak and Satterthwaite, 1956), thermal conductivity 
(Goodman 1953, Laredo 1955, Montgomery 1957), acoustic absorption 
(Morse and Bohm 1957), electromagnetic absorption (Glover and Tinkham 
1956, 1957, Biondi et al. 1958, Richards and Tinkham 1958), and spin 
relaxation time (Hebel and Slichter 1957). 

For these problems, the size of the energy gap acts as the ‘order para- 
meter’ or ‘internal field’. To this extent there is a parallelism between 
superconductivity and other cooperative phenomena like ferromagnetism. 
The electrodynamical properties of superconductors are less well under- 
stood. The mere existence of an energy gap at the Fermi surface does not of 
itself imply superconductivity—it is equally characteristic of insulators. 
The vanishing of resistance is, perhaps, plausible, since collective states 
of motion exist in which the electron system, as a whole, has a non-zero 
momentum. These states can have energy as close as we wish to the ground 
state, since they are states of a system of macroscopic mass. These curl- 
free current distributions may be metastable, since no single-particle 
transition will be able to destroy the current. 

The derivation of the Meissner effect—i.e. the totally diamagnetic 
property, or the first of the London equations (1.1), is more difficult. 
Without a direct physical picture of the significance of the ‘internal field’ 
parameter, one cannot hope to get much insight into how the magnetic 
field is excluded from a body, and hence no guidance in making useful 
approximations. In this situation it may be dangerous to make any 
approximation at all; Schafroth (1951, 1958) has shown that unless the 
derivation is strictly gauge-invariant, the predicted Meissner effect can be 
modified arbitrarily by a suitable gauge transformation. 


§ 2. Tue Hamiirontan ForMULATION 
2.1. Outline of the Bloch Theory 


In Bloch’s (1928) theory of conduction in metals, the electron—electron 
interactions are neglected, and the lattice is regarded, in the first instance, 
as providing only a static spatially periodic fieldt}. The time dependent 


Rinse EE ee ee 

+ Since the electron density will also have the periodicity of the lattice, a 
Hartree self-consistent treatment of the electron—electron interactions is 
(Bardeen 1937), in principle, possible within the framework of the Bloch 
theory. The effect of electronic correlations cannot be handled thus, however, 
but have to be treated by, for example, the Bohm-—Pines plasma theory 
(Bohm and Pines 1953). 
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Schrédinger equation for the entire electron system separates into single- 
particle Schrédinger equations 
h? ; 

et Diy; V2 dari) + V (05) Pret) = Byhu (ti) + Sne large Sen) 
where r, is the coordinate of the ith electron, #, and H, are an eigenstate and 
the corresponding energy eigenvalue. Periodic boundary conditions will be 
assumed, in a box of volume Q, containing an integral number VQ oflattice 
cells. The Pauli exclusion principle is satisfied by writing, for the total 
wave function, a Slater determinant, 


Y= det #,(r;). . eee aed, eeirecaneme eenee 
In (2.1), V(r) has the periodicity of the lattice, i.e. 
Vir)=V(r+nyat+nob+ngce). . . . . . (2.3) 


where 7, %9, 3 are integers, and a, b, c are the basic lattice vectors. Using 
the periodicity of V, it is easily shown that the solution y, have the form 


o(r)=exp (ik. r)2,(r) | +995 eee 
where w,(1) has the periodicity of the lattice, and where the components of k 
are integer multiples of (27/Q)".  Sufficiently close to an atomic nucleus, 
u,(r) is, of course, asymptotically the wave function of an electron in a free 
atom. 

In the absence of any lattice irregularity, these solutions (2.4) are exact 
eigenstates, and the electrons suffer no scattering. The ‘residual’ resistance 
of metals arises from scattering of Bloch electrons by static lattice imper- 
fections, and will not be of interest here. The ‘ideal’ resistance, in 
contrast, is a dynamic effect, arising from the thermal motion of the ions. 

If the displacement field of the ions from their equilibrium positions is 
denoted P(r), the change in potential of an electron at a point r will be 
approximately 

Vis =P. VV (r)y) 2 eee eee 


the model has been implicitly extended to include an interaction between 
the electrons and the longitudinal displacement field of the lattice. The 
lattice displacement field in isolation may, in the usual way, be resolved 
into a set of harmonic oscillators. The interaction will then lead to 
transitions, in which a ‘phonon’ (a quantum of excitation of a single 
normal mode of the lattice) is absorbed or emitted by an electron. In these 
transitions, the total wave number is conserved, apart from integer multiples 
of the basic vectors of the reciprocal latticet. 
The matrix elements for transitions from a state with n, phonons to one 
with n, +1 phonons are 
(k,n, | Vi [k+w,n, —1)=20fhw,, /NQMs,2}2 n V2, 
9 
(k,n, [Vi |[k—w, 2, + 1) =30fhw,/NQMs,22(n, + 1)¥2 ay 


} Transitions in which this vector in the reciprocal lattice is different from 
zero are known as ‘‘ umklappprozesse ”’ (Peierls 1932). 
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(see, e.g., Sommerfeld and Bethe 1933), where WN is the number of atoms per 
unit volume, M the mass of anion, s, the sound velocity of the ion system and 


he 
o=3{ ——jeradu,(r)|?d’r. . . . . (2:7 

lattice cel 2 lg ("| eth 
C is only slightly dependent on k and will be treated as a constant energy : 
it is of the same order of magnitude as the Fermi energy (,. 


2.2. The Quasi-Free Electron Model 


In Bloch’s theory, although the transition matrix elements are calculated, 
no explicit Hamiltonian is ever writtendown. Frdéhlich (1952) showed that 
it is advantageous to rewrite the Bloch theory in Hamiltonian form; it isa 
field theory, in which the role of ‘ virtual’ transitions may be investigated by 
field-theoretical techniques, such as, for example, unitary transformations. 

Bloch’s matrix elements (2.6) follow from the Hamiltonian 


Fi Eon ting te 2 oom? (2.8) 
where the energy of the electrons is 
Ha= > & DB ko Uicas . . . . . . . (2.9) 


<, is the energy of a single Bloch electron of wave number k, at,,, a,, are 
canonical creation and destruction operators for an electron of wave number 
k, spin o, and satisfy the anti-commutation rules 


Ca Onrge® \4. = [ay Byer ge Ly. => 0, | 
en ss 


It is convenient to define 


(2.10) 


Ficc= Mo eg St A ee? eae) 


the number of electrons in the state k, o; /,, has eigenvalues 0,1. The 
lattice energy may be written 


Fy non= LhOw (Oy by + 3)s Ce Keo (2.12) 


where b,*, b, are canonical Boson creation and destruction operators, 


satisfying the commutation rules 
by +, Oye] = [By Oy] = 9, 
[Ow = [Ow, Ow? (2.13) 
Low, bw] = Ck w? 


and fia, = hws, is the energy of a phonon; 5, is the velocity of sound in the 
ion system in the absence of interactions (and may depend,on w). The 
number of phonons of wave number w is 


Roe Os We ae ae eee ees 
The lattice—electron interaction is 
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where c.c. denotes the complex conjugate of the preceding quantity, 
4C*hw 


te ee ta 
Pw 9NOQMs,? 
and p,, is the Fourier coefficient of the electron density, 
Pw = D> Ousw.en Oy eee eh 
k,o 
It is convenient to introduce the dimensionless coupling constant 
F=Q8D,3/hw,, 2 kai Se eieas Be 


where the number of energy levels per unit energy per unit volume at the 


Fermi surface is 
(Qn) =i /(S)| eae eae 
Ky 


Finally, in (2.8), the Coulomb term 
AA aca = > (27re?/q?) pe. Pq 3 . < : ° (2.20) 
q 


has been included for completeness. Its role will be discussed in § 3.4. 

A useful simplification is the ‘ quasi-free’ electron approximation. The 
single-particle energy ¢, is taken to be quadratic in k, e,=A+h?k?/2m, 
where m is an effective mass. Umklappprozesse will be neglected, and the 
Bloch wave functions exp (ik.r) u(r) replaced by free particle wave func- 
tions expi(k.r). Strictly, this simplification is not consistent, since if 
u(r)=1, then from (2.7) 

Gm). 4. Ue ee oe ek 


and there is no interaction. If however one regards C as a fixed constant, 
then the simplification may be made. In (2.9), (2.10), a,,*. a, are now to 
be regarded as creation and destruction operators for free particles of mass 
m, and €,=h?k?/2m. In the quasi-free electron model, the density of 
levels (2.19) takes the form 


E=3N JQ=iN/[G . - . . . . (2.22) 


where V,=number of free electrons per unit volume, y=number of free 
electrons per atom. ‘The problem of the interaction of electrons with. 
lattice vibrations has now been formulated completely as a field theory. 


2.3. Perturbation Theory 
The role of H;,, (2.15), in the first place, was to produce real transitions 
between Bloch states, necessary for the ideal resistance of normal metals. 
Besides producing real transitions, however, it will lead (through virtual 
transitions) to a modification of the energy eigenstates. If it were not for 
the fact that the Bloch states are highly degenerate (forming as they do a 
Ce at 


+ This /' differs from that of Frohlich (1950, 1952, 1954) by a factor v, the 
number of free electrons per atom. E 
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three-dimensional continuum), there would be no difficulty in treating the 
interaction as a small perturbation. In the actual system, with its 
degeneracy, one must first find appropriate linear combinations of the 
Bloch wave functions to lift the degeneracy+. 

Nevertheless, we might hope to get some qualitative and order-of- 
magnitude indications by simply writing down the usual second order 
energy perturbation, neglecting the presence of degeneracy (Fréhlich 1950): 


AH= — = Jee (Fae) |P o” I/(€, —€, +hw,,), pao ©4053) 

where w=k—q. The term in f,, 
AE, = — x jel Pw /(e,-egthw,) . . . (2.24) 
simply leads to a renormalization of the energy levels, but the term in f, Sg 


pea Freotao | Pw P/(€x— €g +h ,,) eee (too) 
,4,0 


represents a dynamic interaction between electrons. It will be shown 
that this expression merely represents the diagonal part of the effective 
electron—electron interaction, and that the non-diagonal matrix elements 
will be much more important. 

If the interaction is sufficiently strong it may be that the Bloch ground 
state is no longer the lowest state of the system; some sort of condensation 
may occur. Clearly, the larger the interaction constant F («| D,,|?), the 
greater is the likelihood of existence of a ground state differing from the 
Bloch ground state. In other words, the higher the normal resistance, the 
more favourable for the occurrence of superconductivity. One canestimate 
the order of magnitude of the energy of the condensed state by noting that 
small energy denominators, and hence large contributions to AH, occur 
only when |<,—<«,| </w). Moreover, only Bloch states within ~hw, of 
the Fermi surface can be involved in the condensation, since otherwise the 
increase in kinetic energy of the electrons will be unprofitable. Hence, if 
hie is an average phonon energy, 

AE,~ Dp | D,,? |/he@ 
le,— fol < ha 
leg— Col <a 
~ B20? (h@)?C?/NMs2Q 
~N Q(C?/fo)(m/M) eR es Ree 20) 


(remembering that the maximum wave number wp ~ ko, so that ha ~ hk gq). 
The expression (2.26) shows the isotope effect, since AH, cm/M, H,«M au 
The smallness of the superconducting energy compared with the Fermi 
Spee ee ete Ss en oe ee ee 

+ If the system does not become a superconductor, then perturbation theory, 
with the Bloch wave functions as the unperturbed states, should be reliable. 
The electron-phonon interaction makes a contribution to the linear term in 
the specific heat (Buckingham 1951, Buckingham and Schafroth 1954). ; 


14 C. G. Kuper on the 


energy arises also, from the factor m/M; AH,~ho per electrony. This 
prediction is still too large, by a factor between 10 and 100, compared with 
observed condensation energies. 

Bardeen (1950) proposed the use of a linear combination of a chosen 
Bloch wave function and all the other Bloch states that it can interact 
with in first order in the interval | ¢;,— ¢ | </i@ as a new basic wave function ; 
he hoped that such wave functions would exhibit superconducting 
properties. This approach was formulated variationally, and was shown to 
lead to energy expressions essentially the same as Fréhlich’s (Bardeen 1951). 


2.4. Renormalization of the Phonon Energy and Coupling Strength 


The indications from perturbation theory suggest that there is a critical 
value for F, below which superconductivity cannot arise. Although this 
suggestion turns out not to be true (see §§ 3.2, 3.3), it has stimulated some 
valuable investigations. If indeed F were large enough for the perturba- 
tion theory to give a new ground state, then (Wentzel 1951) application of 
perturbation theory to the energy of a phonon gives a negative total energy, 
implying lattice instability. It has been shown by Huang (1951) and 
Frohlich (1952) that this feature arises, in a perturbative treatment, from 
the neglect of a possible adiabatic adjustment of the lattice to the motion of 
the electrons. To allow for this possibility, a canonical transformation 
will be performed ; new field variables will be introduced, to decouple the 
electrons from the lattice motion as completely as possible (Frohlich 1952, 
Nakajima 1953, Bardeen and Pines 1955). 

Let 


Sy = > 9k, w Cre! Seow. clw fo.) ee nee 
kc 
and Kai vie he 8 9 io tee age ee ee Deas 


where gy, ,, is a real c-number, to be determined. The unitary operator 
exp isis used to transform the Hamiltonian, 


H*? = exp (—iS) H exp (iS) 


=f +14, 8) ((8; 8), S]/2i4e.4 a eee ee) 
To evaluate (2.29), the following commutators are required : 
[H.,,S]= : > (€,— €ik— wilt, witico® Oy, Uw OT - (2.30) 
[Anon S]= . 2 Rody Dew {Ow wot Ne g— 6.0}, « . . . (2.81) 
LP» SI= : 2 Jews taldee byt +¢.c.} 
+ terms non-diagonalin electron variables, —. (2.32) 


—_—_—_—— ee er 

} For a fuller discussion of the qualitative features of the perturbation 
theoretical approach, see Frohlich (1953) and the discussion on Fréhlich’s 
paper. 
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[D.,, Ppwowt —¢.c., S] 
=> D,, (bw Lew, 81+ [bw *; S] owt —e.0.) 
= > Dd, Cite pee, 5) Ix, w Ga Ow) 


k,w,o 


Be pw open. Oy pet ec.c.).. = = = . (2:33) 
With these commutators (2.29) gives 


HT? = A, a Le 7; i > {D., ts Ix, w (€; each icsw | a Teer.) {Ow qh _w, a Ae } 


“s 
»wW,o 


+2 (Dy Me, w + BIW" (Exe — Ei ew — H»p)} 
x iChye Ties, hae OV 0. byt) aaa (Pw Ak _w, ae Ong si c.c.)} 
+ higher order terms. PE eesti Beh. Fee Cit Pes fie )!s7 4 (2.94) 


The coefficients g, ,, are to be chosen so as to eliminate terms linear in 
Ome ine. 1.6. 


w ? 


Ii, w= Pu (Eewi— kN yy)- 5 ie: ee eT Sy) 


However, the series development (2.29) is only valid if the coefficients 
I, w remain finite—an ‘energy shell’ must therefore be excluded from the 
transformation, i.e. we must take g,, , =O if |e,._,,—«;|<A,say. Itcan be 
shown (Fréhlich 1952, Bardeen and Pines 1955) that A can be chosen both 
large enough for rapid convergence of the expansion (2.29) and small 
enough to have negligible influence on the renormalization. Henceforth 
the restriction that (2.35) holds only for |e, _,,,;— €;|<A will be weakened 
to the restriction that when a simple pole occurs in a summand or integrand, 
the principal value is to be understood. 
With this understanding, (2.34) now gives 
HH? = > €,4,¢,, + > 10, (Oytbyt+h)+i > Dy (pwOwt—e.0.) 
k,o w energy shell 
5 ee, 


a 
—4 lie Bt ORE La gee ae eS, 
k,w,o €lk—w 1 _ €, thw, 


(Pw%—w,o' Uo + C.-C.) +higher order terms, 
. (2.36) 


where the terms neglected in (2.36) are at least quadratic and non-diagonal 
in the b,,*, 6, and can be neglected, to within the validity of the harmonic 
approximation for the lattice waves. The prime on the summation denotes 
the principal value, and. 

a Oy pa ie w Cie they 4 


k,w,o 
My — 2 2D) Da" Ico! (1 x —wi— €) » + + (2.37) 
is the renormalized frequency. To within the present approximation, 
iw,, in the denominator may be replaced by /a,, everywhere ; in (2.37) it is 
good enough to neglect the fw in the denominator entirely, and replace f,., 
by its expectation value. The observed velocity of sound, s, is @,,/w. 


16 C. G. Kuper on the 
In terms of the renormalized coupling constant F’,,, defined by 
F LF = (0 ,=8l8 - ~~ — - +) Aeas8) 
= DEK, M gy to + Lheiw (bw* bw +2) 
ko w 


+i > (Dep Cee eo) 
energy shell 
F (hay)? TO _w o Pw - 
+ >! —— oe ee i  — +higher order terns: 
2, DEl(ckan (= ep en ee 
(2.39) 


Wentzel’s difficulty is resolved, for it is F which has to be less than a 
critical value ~ 1 for lattice stability, while ?, which appears in an electron— 
electron interaction term, may be as large as we please. 

This renormalization may also be derived by Bogoliubov’s transforma- 
tion. The advantage of that method is that it leaves the Hamiltonian 
in arather more convenient form for determining the ground state. Since, 
however, even the ‘weak coupling’ limit /’ +0 will lead to superconduc- 
tivity, the renormalization is not essential in principle. It will be necessary 
for a completely quantitative theory, since in fact for most metals F is not 
very small, compared with unity. In the weak coupling limit, the inter- 
action terms in (2.39), a 

po ’ Fhe ww) Bea "Ow, o Pw 
“ oe QE{(e1k—wi— G)ia (RG p)”} 


are just the second order perturbation matrix elements. 

The Coulomb interaction (2.20) has, so far, been neglected entirely. 
Following Bohm and Pines (1953), we may introduce additional coordinates 
for the collective (plasma) modes of the system, together with constraints 
to preserve the correct number of degrees of freedom of the system. Then 
a Bohm-—Pines transformation, to decouple the plasma and individual 
particle modes, leads to a ‘screened’ Coulomb interaction, of short range. 
Nakajima (1953) and Bardeen and Pines (1953) included the Coulomb 
interaction in this way, and subjected it to the transformation (2.29), 


resulting in the appearance on the right-hand side of (2.39) of the additional 
term. 


(2.40) 


> ae Soe Oe yt ae 
k>ke 
The Bohm and Pines approach is not rigorous enough for this result to be 
completely convincing, however. Further discussion of the Coulomb 
interaction will be given in § 3.4. 

The Hamiltonian (2.39) serves as the starting point for the work of 
Bardeen, Cooper and Schrieffer (§3.2). It exhibits an explicit electron— 
electron interaction. The diagonal part is just the energy AH of (2.23) 
but there are many non-diagonal terms of (2.40), of the same order of 


magnitude as the diagonal ones. The actual expressions are modified by 
Bardeen ef al., to simplify the calculation. 
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§ 3. THe GROUND STATE anp ENERGY Gap 
3.1. Frohlich’s one-dimensional Model 


The difficulty in the perturbation theory (§ 2.3) arises from the degeneracy 
of the Bloch states (Frohlich 1953). Perturbation theory alters the effective 
mass of electrons near the Fermi surface, but does not greatly affect the 
density of states for low-lying excited states; one would expect the lifting 
of the degeneracy to modify the density of states considerably. Experi- 
mentally the absence of any term in the specific heat linear in the tempera- 
ture shows that the density of low-lying states is greatly reduced. The 
most recent measurements of the specific heat of tin (Corak et al. 1954, 
1956, Corak and Satterthwaite 1956, Lynton et al. 1957, Zavaritskii 1957) 
confirm the suggestion (Ginsburg 1946, Koppe 1947, 1950, Goodman 
1953) that there is an actual energy gap between the ground state and the 
lowest excited states. 

An electron within ~fw, of the Fermi surface can emit a phonon and 
reabsorb it, returning to its original state. These are the only processes 
considered in the second-order energy perturbation expression. But 
having emitted this quantum it could, instead, absorb any of the other 
phonons present, ending in any of about 3Q.i& states. The probability 
that any given one of these phonons is excited is fairly small, ~ (ia@/C,)"?, 
but the probability that the electron absorbs a different phonon is over- 
whelming compared with the probability that it absorbs the phonon that it, 
itself, emitted (and the ratio is proportional to the volume (Fréhlich 1953)). 

If the coupling is strong enough, one might get very high excitation of a 
very small number of lattice modes; ifso the degeneracy is greatly reduced. 
To try to get some qualitative understanding of the ‘condensation’, 
Frohlich (1954) examined an extreme case—a one-dimensional metal 
with exceedingly large #. The difficulties are associated with the high 
density of states, which is greatly reduced in one dimension. Further, in 
view of the Pauli exclusion principle, the only possible resonances are 
between electrons at opposite ends of the Fermisea. Ifthe energy of recoil 
is neglected, then the lattice modes + wy (Wy) = 2k ) are strongly excited, 
since they are exactly the modes required to produce resonance. It is 
assumed that all other lattice modes can be regarded as a perturbation ; 
this assumption is only valid if the coupling constant F is quite un- 
realistically large. A particle representation, will be used in this section, 
since the physical mechanism will be clearer than in second quantization. 

If for the present all lattice modes other than + wp are neglected, then, 
with a Hartree self-consistent field approximation, the Schrédinger 
equation, separates into the single-electron equations 


(Dp; +i{D yy, (bp, + 5+ —w,) EXP (iors) — ¢-0.}) Pu(Xj)= A ppp(%j) (3-1) 
together with the lattice equations 
| Tedyy, Diy, Big, FIN QD y, (Ben, EXP (i0G@) — ©-0-)}X yy = At Xeoy? 
{fhe ,,, D_1p,* bw, FIN QD 1, (b_1,,* EXP (— 19) — C.C.)}X ay, =AL Xen," 
(3.2) 


P.M.S. .0) 
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Equations (3.2) are Schrédinger equations for oscillators referred to @ 
displaced origin, and it follows at once that 


b= +i(Dy, [hw »,) N Qexp (+ iW 9%) 


= Li (Dy, ary) [ dre*( exp ( F iw px)p,_ (2). a ath 


Equation (3.3) connecting the electron and phonon variables, provides a 
condition of self-consistency. In it the integration is over the normaliza- 
tion volume Q, and the summation is over all occupied states k. 

(3.1) describes a single electron in a sinusoidal field—the very situation 
which gives rise to the band structure in solids; an energy gap results 
precisely asin bandtheory. Ifk~ko, then the solutions of the Schrédinger 
equation are 


thy = Q7” (w,32 exp (ike) +0,42exp [i(k — wo) +i), re 
big = O72 (v,12 exp (ike) +43 exp [i(k— wo) +18) 


where f is a phase angle, k is to be taken greater than zero, u,, v;, real, 


a5 


WP =3 (1+ pel? +0)), . . - « « (8.5) 
by? Pup te= le ee er re 

and 
T=4$D,, |6,,,+5_.,+ |. < ua) Cee, ho ee Se 


The single particle energy, relative to the Fermi surface is (Peierls 1930) 


fo. = ene +o) PVA hw) ee 
(where €,=H;,,—[, n,=«€,—). There is clearly an energy gap 2 
between the ground state and all excited states. 

Although the expressions (3.4), (3.5) and (3.8) are valid only near k= kp, 
they can be extended to all k by perturbation theory (Frohlich 1954). The 
perturbative corrections are negligible in determining I, since it is the k 
near ky that are important. But when k is not near kp, either w;, or v; is 
nearly unity, and the other is small, of order ['/¢. Then the perturbations 
arising from transitions between k and k+w, or k—wy and k—2w, are 
comparable with the smaller term in (3.4). For k_ states, the perturbative 
corrections are 


4Q-1?DP exp [i(k + wy)a — iB] 43 
Ai OE i A eee 
by C5(1 + k]ko) h— Col +k/kp) < Ko, 
30-127 exp [i(k — 2wy)x +i] a aie 
oi, Paine Ueda Ry k> kg. 


(3 — ko) [o(8—K]keo) 


The quantity exp (iw x) appearing in (3.3) is 
eXp (iw r)=(N.Q)+ > pb, * exp (iw) ob, da 
occupied kh 


(i.e. al A —) 


ae 
= (V0) 35 NS, a 
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To evaluate this expression, use is made of the fact that in the ground state 
the k_ states are occupied and the k, states unoccupied. The one-dimen- 
sional expression for converting sums over occupied k to integrals is 


Ly 
‘ = 1,08, | OR ats ae eh 
occupied k 0 
with the aid of which (3.10) is 
exp (Iuyt)= $i (TG) exp (i8)In(4Q/T) . . . (8.12) 


(assuming the smallness of F compared to ¢). Substituting (3.7) and (3.12) 
in (3.3) gives 
Rita (Aca yes “ore = One % Dobes ¥.0'(3)13) 


The non-trivial solution is the one of lower energy, and gives for I’ 
Bere ALOR (i S/ 2H) oF va bey bee ees oe (3,14) 


The total energy, measured from the Bloch ground state is (from (3.3), (3.8), 
(3.9), (3.10) and (3.14)), 


We =2) y+ hier, (ry, * by, + b_o,* 54.) 
k 


SSN, Oxp(— ol). 4) ss as ee (8.15) 

In the calculation of W, there are two comparable contributions—a large 
shift (of order I’) of a small fraction (~ ['/f,) of the energy levels, and a small 
shift (~ T?/€,) of the other levels. For the latter, the perturbative correc- 
tions (3.9) are important. 

The calculation can be extended to a finite temperature (Kuper 1955), 
merely replacing expectation values by ensemble averages. Excited. 
electrons are in k, states, and reduce the amplitude of the electron density 
wave 

Ds EXP (iwpX;) 

and hence also of the lattice displacement,|6,,+6_,,*|. The gap 21'(7)isa 
monotonically falling function of temperature, and vanishes at a tempera- 
ture «7,=1:14T'; above this temperature the Bloch states are restored. 
The specific heat is exponential in 21\/«7' at low temperatures, and has an 
Ehrenfest (1933) second-order transition at 7’). The magnitude of I\(7’) is 
proportional to the intensity of the lattice displacements, and serves as 
an analogue to the ‘internal field’ of the Weiss (1907) theory of ferro- 
magnetism. 

The expression (3.15) for the condensation energy does not show any 
isotope effect. This is a consequence of the fact that the recoil energy of the 
lattice has been neglected—the b+, b.,,, have been treated as classical 
field variables. A condition for the validity of this approximation is that 
the gap must be so large that no other phonon mode w can cause an electron. 
to bridge it, i.e.%o,,,,<2I. But the total condensation energy per electron 
~xT,. The condition may therefore be written | 

TOMTOM eel eh aik, A, 2 aie (32 LO) 


C2 
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(where © is the Debye temperature), which is not satisfied by any real metal. 
The method cannot be generalized either to weaker interactions or to three 
dimensions, and its main value therefore lies in showing that the electron— 
phonon interaction can produce a condensation, and in suggesting the ex- 
ponential dependence of the energy on the coupling constant—a result 
which of course cannot arise in any finite order perturbation theory. 


3.2. The Three-dimensional Model of Bardeen, Cooper and Schrieffer 


The one-dimensional model has shown the possibility of finding a 
solution not expressible in a power series—exp(—1/f’) has an essential 
singularity at / =0 and cannot be expanded. In addition exp(—1/f) isa 
very sensitive function of F, and if it arises in the three-dimensional theory 
it would explain the wide scatter of superconducting critical temperatures 
(e.g. Shoenberg 1952), and also the large pressure coefficient (Kan et al. 1948, 
Chester and Jones 1953, Muench 1955, Liithi and Rohrer 1958) of the 
critical temperature; if the factor exp (—1/F) can be combined with the 
‘isotope effect’ factor (A)? of (2.26), many of the objectives of the theory 
are fulfilled. 

The renormalized Hamiltonian (2.39) exhibits an explicit electron— 
electron interaction which is negative (attractive) when the interacting 
electrons have energies differing by less than ~/0,,,,,,,, positive (repulsive) 
otherwise. Clearly the ground state should be one in which the attractions 
are as effective as possible, and the repulsions as ineffective as possible. 
One might even hope that in the ground state the repulsive interactions 
might be negligible. If this guess should prove correct, then it is possible 
to simplify the Hamiltonian considerably without loss of any essential 
features. This is the approach of Bardeen et al. (1957)—they replace &, 
by an average value # and simplify the interaction (2.40) 


F 
, x = 
(a= =O, > Hep ei pw). aw elena Gait eo Lees 
=F qa, w, 
0,0" 


by taking (broken curve, see figure) 


G=—1,|e,—€e |\<ho 

? C Ik—wl ? 
‘4 (3.18) 

< 

G=0, lex — in —wi | >a, 

instead of the actual expression (continuous curve, see figure) 
ho)" 

G= we) laut SATO} 


This inter-electron attraction makes it possible for electrons near the 
Fermi surface to condense into quasi-molecular pairs with zero total net 
momentum, as conjectured by Schafroth et al. (1957). Cooper (1956) 
calculated the binding energy and wave function of one of these pairs, how- 
ever, and found that the diameter of a quasi-molecule is ~10-*em; the 
overlap of quasi-molecules is important and prevents us from treating the 
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system of pairs as an ideal Bose gas. Let us consider a single state 4, of 
energy €o, of the electron pair, interacting through a constant matrix element 
V with NQ states ¢,, energies ¢,, in a narrow energy shell of thickness 
~haw. Then the Schrédinger equation for the state b= ¥c,,¢,, + Cob) may 
be written 

Co€oho + >CnV bn = Leody ; Mites Pees. 4003220) 


Cn€n Pn + Cy V by = H Cn dn: 39 ot Boon Ag he ea (3.21) 


a) — 


[&k — £1k-w) | —> 


There are two typesofsolution. Ifc)=0, then from (3.21), H =e, and hence 
from (3.20) 
SCV bi 0. eee i te eee) 


This is a single restriction on the NQ coefficients c,,, so that there will be 
NQ —1 solutions of this type with energies in the shell. The other two 
solutions are outside the continuum. If He, is assumed, then from 
(3.21) 

Cyr Coal 6, one een Sle ponte ean (3-23) 
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so that, from (3.20) 


e+thw 
=HOV? | «/(E—e) 
€9—tlio 
= —BOV? in |(B—e,— fiw)/(W—e+ dia), «6» (824) 


where 2 is the density of states (per unit volume). One of these states, the 
Cooper quasi-molecule, lies below the continuum (as may easily be seen in 
the particular case of an infinitely thin energy shell, when (3.24) reduces to a 
quadratic equation for /). 

The situation with many electrons is rather more complicated—the 
Pauli exclusion principle prevents one from simply pairing off all the 
electrons in Cooper pairs in a straightforward way. There is however one 
way in which it can be done fairly easily: all the pairs must have the same 
net momentum, and for the ground state this momentum must of course be 
zero. When exchange effects are taken into account, the antiparallel spin 
correlations are favoured compared to parallel spins. 

For the moment, only some of the non-diagonal terms of the interaction 
will be considered—those in which (see (3.17)) k—w=—q: a=o'. These 
terms scatter the pair (q, t ; —q, ) )into(k, + ; —k,\). This ‘reduced’ 
interaction Hamiltonian, of course, contains only a small fraction of the 
number of terms in (3.17). However all pairs (k, + ; —k, ) ) are connected 
by terms in. 77,4, ..q and all terms not included in H.,,,,,,,..g will break pair 
correlations. The simplest improvement on the Blozh wave functions 
would seem to be functions containing two-body correlations but no higher 
ones ; for such functions, H,,.4,,..q Will play an important role but H’ — H 
will not. 

It is convenient to consider a state which is not an eigenstate of the total 
number of particles VQ = S'a,,*a,,, i.e. the number of particles is not pre- 
cisely specified. Since }a,,*+a,, commutes with the Hamiltonian, the exact 
ground state is certainly an eigenstate of V,Q. IRfnecessary, we may always 
take the projection of the approximate ground state on to the subspace of 
N,Q particles. 

As a variational ansatz, consider the wave function of Bardeen ef al. 
(1957) (where |0) is the electron vacuum) 


reduced 


PHT (uy + M4 ty *)]0). eee eas oe) 


It exhibits the desired pairwisec orrelations, and the parameters w, UY, 
will be chosen to minimize the energy, subject to the normalization condition 
Uj,2+0,2= Le e e e e . ° ° e (3.26) 


Symmetry will be assumed between electrons above the Fermi surface and 


holes below it, 1.¢. #;, + Uj. = life, + €,,=2Cy. Interms of energies measured 
from the Fermi surface, 


T= tn— Cg: Ge cee, 
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the expectation kinetic energy for the state (3.25) is 
2 yy Nae + 2 > | pe |eene2s ° ° . ° e (3.28) 
n> 0 Ny <9 


and the expectation interaction energy is 


F 
On Me Hee Hn — Ie). ‘ S : . ‘ (3.29) 


Only those terms in H’ contribute which belong to Hyoajooq> i-€. in which 
the pairwise correlations remain unbroken; the remaining terms of H’ 
have zero expectation value. Minimization of the expectation value of 
the energy leads to an integral equation for w,, v, 


40 E+ Lhe #/QZ)GE(g,— Ne (Up? — 047) U1, = 0, 
i.e. replacing the sum by an integral in the usual way, 


Aaa = (22 VF | ere Horne) dn - . (3.30) 


This integral equation is difficult to solve as it stands, even with the 
simplified interaction (3.18). For that reason the further simplification 
of Bardeen ef al., 


G21 if* Inpl<%a, f,| <ha, 
(3.31) 


G=0 otherwise, 


will be made for the present. With this factorizable potential, the 
integral equation is 
~ ha 
realli? — 02) = (F200) | Upp Any  . » » (3.32) 
0 
which is easily solved. The solution is (where evidently the + sign is 
appropriate for n;, <0) 


P= 51 Enel (ne +1), 8... ee (8.38) 
[T=(F/QE) > u,v,-=he cosech (1/F), . . . . (3.34) 
a 


giving for the ground state energy (taking the normal state as zero) 
W = —28.0(fa)2/(exp (2/F)—-1).. . . . . (8.35) 


However weak the electron—lattice coupling (provided only > 0), 
the ground state of Bardeen ef al. has a lower energy than the Bloch ground 
state. (This result is to be contrasted with the perturbation theory, 
which only gave a distinct ground state if # was sufficiently large.) In 
the limit of weak coupling 


ener expi( 1 LE laren a cee ee 300) 
WA Ozia)? exp (32/F).4: 2 Gee e.37) 


Wave functions for low-lying excited states can be constructed by similar 
arguments, and Bardeen et al. showed that there is an energy gap 21° 
between the ground state and any ‘single-particle’ excited state. The 
derivation of the excitation spectrum is most easily performed by the 
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method of Bogoliubov (1958a, b), and will accordingly be deferred to 
Besides these particle excitations, collective (current carrying) excita- 
tions can be generated, by building up states in which the pairing is 
(k, + ; —k+q, ¥ ), where q isaconstant. These states will have essentially 
the same interaction energy as the ground state, and will differ from it 
only by the kinetic energy of the electron current. In contrast, any single 
particle excitation breaks the coupling of an electron to all the Cooper 
pairs, producing the finite energy gap. 

With the accurate form (3.19) for the interaction, the integral equation 
is more difficult to solve. What one can do fairly simply instead, is to 
use the solution (3.33), (3.34), in (3.25), as a trial wave function, with the 
accurate expression (3.19) for the interaction Hamiltonian ; [is regarded as 
a parameter to be chosen to minimize the energy (Kuper 1958). In the 
weak coupling limit, the exponential dependence of the ground state 
energy on F is again found, though with some slight numerical 
differences. At the same time, that method shows that the diagonal 
terms of (3.17) are not important. 

The results of the theory of Bardeen et al. (3.36) and (3.37) resemble 
(3.14) and (3.15) of the one-dimensional theory, with the ‘isotope-effect ’ 
factor h@ in place of ¢). This is just the combination needed to reconcile 
the isotope effect and the need for a very sensitive parameter. Moreover, 
if F is even moderately small, exp (— 1//) may be very small, explaining 
why 7')<© rather than 7, ~ ©, as one might otherwise expect. For real 
superconductors, the relation between 7’, and © suggests F~1. F is 
fairly small, so that the weak coupling limit should be a reasonable first 
approximation. 

In the one-dimensional theory, if 7) is not 2 ©, the approximation 
method fails. Qualitatively it can no longer be true that only one pair 
of lattice modes is excited; instead we expect a narrow range of wave 
numbers near ky to contribute. A finite range of coherence may be 
expected, since the resonant lattice vibration should have the form of a 
wave packet of finite extent instead of a plane wave. 

In the light of the results of Bardeen et al. one might make the ad hoc 
restriction in the one-dimensional case that only states within A@ of the 
Fermi surface should be taken into account in producing the gap (that is to 
say electrons not within ~ha@ of the surface follow the lattice adiabatically). 
In that case, all integrals should be restricted by the condition 


| — ko|/ky < $ha/ Co. 
Equation (3.10) is modified to read 


exp(iwyt)=(N.Q)* > dP exp (i8)/4/(T2+,2) 


Ing l<ha 


ko(1 + Na/2E5) : 
= (2t| iD exp (i8)dk/2./(T? + 7,2) 


ko 1 — Na@/2Eq) 


=$1 exp (i8)(I'/¢,) ln (20) een ee (3.38) 
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assuming D<2h@. Hence, in place of (3.13) we have 
2F ln (2ha/T)=1 or T=0, 
the former alternative giving 
R21) OP (i 3/22 aes...) ee. (3.39) 


Although (3.36) and (3.39) are formally almost identical, the wave functions 
are very different. In terms of second quantization operators the trans- 
formation leading to the wave function (3.4) of the one-dimensional model 
is of the form 

Op g = Uy Ons ay exp(if)v;,"7a;_ Wy, 8? 

(3.40) 

Ones V_ POs sit exp(iB)w;? ay w,, 3 J 
with no mixing of creation and destruction operators (cf. (3.41)). The 
wave function is always an eigenstate of particle number. 

Indeed the apparent non-conservation of particle number in the three- 
dimensional theory lies behind the difficulty in forming a physical picture 
of an ‘internal field’, such as one has in one dimension (a difficulty which, 
in turn, is responsible for the difficulty (§ 4.2) in constructing a theory of 
the electromagnetic properties). If we wish to form an excited particle, 
we ought, simultaneously, to produce an excited hole. The electrons which 
have to be displaced in k-space, together with the lattice displacement 
that they carry with them adiabatically, produce the internal field. In the 
wave function (3.25), with an indefinite number of particles, this picture is 
lost, since the wave function does not contain information concerning where 
an electron has ‘come from’. 

An intrinsic drawback of all variational methods is that there is no 
internal criterion for the accuracy of the solutions. Moreover the 
Hamiltonian (2.39) to which the method is applied is itself only correct 
up to order F. Since the solutions are not expressible in a power series 
in F, it is not obvious that they have any significance applied to the 
original Hamiltonian (2.8). The trial wave function (3.25) is designed to 
use only a small subset of the terms in the (renormalized) interaction, and 
it is at least possible that a more complicated wave function might have 
considerably lower energy. For all these reasons, a systematic method of 
approximation is desirable, in place of the ad hoc method of this section : 
the ‘quasi-particle’ transformation of Bogoliubov (1958a, b), to be 
described in § 3.3, provides this systematic approach. 


3.3. Bogoliubov’s ‘ quasi-particle’ Transformation 


The considerations outlined in §3.2 make the existence of a condensed 
state very plausible, and give the desirable exponential dependence of the 
condensation energy on 1/F. So many approximations and simplifica- 
tions are however needed that the method fails to be completely convincing. 
The best justification, as well as the simplest derivation of the excitation 
spectrum, follows from the ‘ quasi-particle’ transformation (Bogoliubov 
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1958 a, b, Valatin 1958). This is an extension of Bogoliubov’s (1947) 
Boson quasi-particle transformation (see (3.66)) developed originally for 
the study of the superfluidity of helium. Let us definet 
Opes = UWA, — VA y*s 
Oye = Wye pT — VA by» bo Je hee ee es 
Oy = Uy My +VEAn*s 
e+ = Ugly + VKA_KY 
with 
Uj, = U_ ke V7,=U_ yp. e ° . . . ° . (3.41’) 
It is easily verified that the a, .*, a, variables satisfy canonical Fermion 
anticommutation rules, provided the real numbers u;,, v, satisfy (cf. 
(3.6), (3.26)) 
U2 + Opa Ly ol YS) ee ee 
A particular choice of the ~,, v,, which gives familiar results is 


Up=1,0,.=0 if e.>G; 
k k k> So (3.43) 
U0; = 1 1h le, = G. 
With this choice, a, + represents a creation operator for an electron (k, * ) 
if €,,>¢,, and for a hole (—k J ) if €, <Q. 

The transformation (3.41) does not conserve particle number. To take 
account of particle number conservation, 

Ne) a, tac = N On ne ee 
k,o 

aterm —¢/ is included in the Hamiltonian, where the chemical potential 
¢ is determined by the condition that the expectation number of particles 


Mie N Di oe, 


In terms of these new variables, the Hamiltonian (2.8) becomes (still 
neglecting the Coulomb term) 


H=%7+H,.+H,+H,+H3 . .. . . (3.46) 
where 
P=23 (6 70s WE es are er es 
Ay= 2Aee—¢ (abe? — Vp? oye 4 Foye 4 + Me, FT %e, —} 
+ Yea (duty + 3), Menges we ee eee 
H,=25( (€,— Cpr fey, toy, + +014, oy, Js Care es 
Hy =i 22) Dy hU 1. tet wile, 4 jew, 1 + Oe ow, +7 %— *) 
“hy P(e Oly “FO ong gle es ep en ee 
H,= 1X Prcttlatt sw ile btw, + Oe, —F Ow, —) 
HOEY ew (paw, — 1%, — + jew, + Oa, 4) owt —0.c. . (3.51) 


I Oe 
} The notation differs slightly from Bogoliubov, and helps to make explicit 
the conservation of wave vector. 
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‘Two procedures have been used for the selection of the parameters wu Dea. 
One may again proceed variationally, selecting uw; ,v, to minimize the 
expectation energy—this is Valatin’s method. It depends on knowing 
the general character of the ground state, in the same way as the method 
of Bardeen et al. The alternative, due to Bogoliubov (1958 a), will be 
developed here. It has the advantage of yielding the form of the ground 
state. Valatin’s approach will prove useful in $3.5, in the study of the 
thermal properties of the system. 

If we attempt a perturbation calculation, difficulties arise on account 
of vanishing energy denominators. The ‘dangerous’ denominators are 
those in which two quasi-particles are created from vacuum and no 
phonons, and do not arise before transformation, below fourth order. 
Bogoliubov’s method consists in selecting the parameters w,,, v;, 80 that the 
‘dangerous’ denominators do not arise. In principle one would like to 
suppress them in all orders. In actual fact we will merely arrange for 
their cancellation in second order. This proves sufficient, since an energy 
gap appears in the excitation spectrum. Higher order processes are 
therefore unlikely to contribute ‘dangerous’ denominators after the 
transformation is completed. 

The terms H,, creating two quasi-particles without phonons, are 
‘dangerous’ already in first order. On the other hand, Hs conserves 
quasi-particle number, merely creating or destroying a phonon, and 
H, creates or destroys two quasi-particles with emission or absorption 
of a phonon. Hence in second order there is a ‘dangerous’ process in 
which H, creates two quasi-particles and a phonon, and H; absorbs the 
phonon and scatters one of the quasi-particles. Bogoliubov’s method 
arranges for the mutual compensation of these two ‘dangerous’ processes. 
Writing |V) for the quasi-particle ‘vacuum’ (i.e. the ground state of H,) 
we require 


{H, —H3H,“*H,}|V ) =0, eer Ae eas) 
i.e. 
2(€;,— C)upY x : . 
As D 2 Ut wiik+wi(Uy — Vz") — (Wier wi —Vik+wl Jee ie (3.53) 
wales Ey t+ Bi nswithory 
‘defining 
Ey,= (eg — (uj? — 072) = 11442 — Vn") 
EL, =FidUe —%2); ~ 2 6:64) 


Nk =9,—> DiP(hierw P= wil (Bet Ficewi +ho,,). 


Replacing the 7,,=«,—¢ in the denominator by 7%;, (making an error of 
fourth order only in the D,,), and converting the sum to an integral, this 
gives 

Dy? ew View \d3w 


5 A (Cane 
Bit Binrwithor 


QU IY = {Q] (227)? }(u,2 — 0,2) 
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an integral equation (cf. (3.30)) to determine Wes Vie Besides the trivial 
solution u,v, = 0, it has another one, the one leading to the superconducting 
state. Guided by the simplified solution (3.33), (3.34) we may define 


r(ey= 2 [ Dottinewitierwidly (3.66) 
(27)° Bit Bictwit ho 

to find om ‘at 
U2 = HL + Hy (AR? + P(K)”) ‘ ; (3.57) 

0,2 = HL — Fel (Ha + TRY?) } 

From (3.55) ['(k) satisfies the integral equation 

2 Zs 

Ck) = QO | aes) ; (3.58) 
(277)8 (Bit Bicswithoy)/(C ([k + W]) + Fi tezw! ) 


For weak coupling (70), this equation has an asymptotic solution 
exponential in F~!, as may be seen thus: It will be assumed that I(4) 
is a slowly varying function, and can be regarded as constant in the range 
0 <%,<hw,; moreover the hw,, will be replaced by an average value ha 
(cf. §3.2). Then from (3.58) 


- ha D(|k + w])di, kiwl 

0 EF, Fl E kiwit hoy, / (T?(|k 4 w]) + Thetw") 
The main contribution to the right hand side arises when 7,,.,,,;< 3 @ 
(since for Ay .ywiosh@ the 7,,,, in the denominator dominates). 


Within the important region, 7,,,,,, Will be neglected compared to. 
had y/(T? +77 ~1y 1), enabling the h@ to be cancelled, giving 


bla ~ bha/T 
r=r [ oe re | ds/,/(1 +82) 
0 » 


T(k)=F (3.59) 


VJ (0? +9) 0 
= PF Intel,” = 3% cee ee 
1.e. 
Detia exp (lil) a 2°. ee 


The chemical potential ¢ is determined by the condition that the total 
number of particles should have expectation value N .Q—this gives, in 
good approximation, €= Co. 

The integral equation for [ can be solved without the approximation 
of constant w, and gives the same result, with an explicit expression for: 
the average @. Apart froma factor 2 in (3.61), the resulting expressions 
for ['(k), w,, v;, are identical with those of Bardeen et al. (see (3.34) and 
(3.35)). The ground state, expressed in the new variables, is of course. 
the quasi-particle vacuum. Expressed in real electron variables, it is. 
just (3.25), as can easily be seen by verifying that a, annihilates the 
state (3.25). 

The energy of the ground state will now be calculated, without the. 
approximations of Bardeen, Cooper and Schrieffer. It is 


W=0—(V|H,H HIV)... «2 . . (8.62) 
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Hf; gives no contribution to the ground state energy in second order, as 
it creates and destroys a particle; while H,, which creates two particles, 
has already been compensated (to second order in the wave functions) by 
the consecutive action of H, and H;. The expression (3.62) is very similar 
in form to that of Bardeen et al., since the matrix elements of H,H,—H, 
between Bloch states are just the matrix elements of the renormalized 
interaction used by Bardeen. The evaluation of the energy is straight- 
forward, using the expressions (3.57), (3.61) and gives (Bogoliubov 1958 a) 


Washo Oe exp (— 2/9), 2) 2. 2 . (8/63) 
in agreement with (3.37), if @ is replaced by 2o (cf. (3.61)). 
The elementary excitations of the system are generated by the quasi- 


particle creation operators 4,,*, out of the ground state. The energy 
é;, of a quasi-particle with wave number k is, to second order, 


6, = « V |oy.,(H <3 LN )ou.5* | V) ae « V oj 3(Ho oa Ey) A 3%" | VO 


=; {1- ee eal 
w (Bingwithoy)?— By? 


$25 ab tlrvcrw Wi rexw Dy" repwi + RO w) _ .. (8,64) 
w (Hi idiwit 0@,,)°— 7" 


The first sum vanishes as k->k,, on account of the factor #,,, but the second 
remains finite and near the Fermi surface tends to 


2 >, Ur terw Mi teew wl (Fi rtwit 2Ow) 


~ 2u,0,1(k) ~ D(ky) =ha exp (—1/F). at Perel s.00) 
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Thus quasi-particle excitations have energies separated from the ground 
state by I'(k,); to create a quasi-particle pair requires at least 21(k,). 
As in the simplified theory, there are also low-lying current-carrying states, 
which will be discussed in § 4.1, since they are important in understanding 
the difference between a superconductor and an insulator. 

The above calculation has used Fréhlich’s unrenormalized Hamiltonian. 
Bogoliubov (1958 b) has shown that similar results may be obtained from 
the renormalized Hamiltonian, using the same transformation. The most 
satisfactory procedure, however, is to carry out a Boson transformation 
(Bogoliubov 1947) on the phonon variables 


Bye =A Ore + HOt» 
Beene + ppd_¢, ° ° . . e Py (3.66) 
Ay? — py = 1, 

simultaneously with the Fermion transformation (3.41) (Bogoliubov et al. 
1958 a). Using the condition that in second order the terms 1, + H,+ Hs 
do not contribute to the ‘quasi-phonon’ energy, the coefficients Ars Hr 
are determined, and the expression (2.37) for the phonon frequencies is 
obtained. The calculation of electronic energy levels follows as above, 
with F replaced by F. 
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3.4. The Role of Coulomb Interactions 


Qualitatively, it is very plausible that the Coulomb repulsions between 
electrons tend to prevent the formation of quasi-molecular pairs, and 
thus to hinder the occurrence of superconductivity, as suggested by Bardeen 
et al. A rigorous treatment of the Coulomb interaction is difficult, on 
account of the long range of the forces, but a rough order of magnitude 
can, be obtained fairly easily. 

It is well known (for example from the theory of a quantum mechanical 
plasma (Bohm and Pines 1953, Pines 1953, Bardeen and Pines 1955)) 
that the correlations between electrons lead to a ‘screening’ of the inter- 
actiont. For a preliminary discussion of the role of Coulomb interactions: 
in superconductivity, the screened interaction (which is of short range) 
will be used instead of the complete unscreened form. The structure, 
in Fermion field operators, of the Coulomb interaction (2.20) is the same 
as that of the ‘phonon’ interaction between electrons (2.40), but with 
the very different function 


G, = (35 /4F)27e2/w%, w>k, 
=()) w<k,, a, 3 tet Sel dee 


replacing G of (3.17). Gis a very anisotropic function of k, k+w but in 
evaluating the expectation energy of an isotropic state such as (3.25), an 
angular average isimplied. This average over angles is almost independent. 
of the energies 7;,, 7,4 (measured from the Fermi surface) (Kuper 1958), 
since for the important terms, we already know that 7 <a@. The intro- 
duction of screening before taking the average over angles simplifies 
the calculation greatly; in the sixfold integral analogous to (3.29), 
no divergences occur. The particular Bohm—Pines form (3.67) for the 
screening is not essential. Any modification of the factor w-, preventing 
the divergence for small w, will suffice. 

The constant (averaged) Coulomb matrix element can now be grouped 
with the interaction term (3.17) of the Bardeen—Cooper-Schrieffer 
calculation. A superconducting ground state will be expected only if 
G,.<1 (since only then will G+G@, be negative for sufficiently small 
i Nk+wi) The criterion for superconductivity is 


OAS 0 Oe cage cay ee yee ae 
where A is 


A= {$F /EO}G,) average 
~ e78/ky?, gh Se ee ee re ree RE 
since the cut-off £, is comparable with ko. 


eee 


} Chester and Houghton (1959) have recently recalculated the phonon 
frequencies and interaction matrix elements, using a Born-Oppenheimer 
expansion in powers of (m/M)!/4. The results are almost identical with those 


of Bardeen and Pines, and provide a rather better justification for the ‘ cut-off ° 
Coulomb interaction. 
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For a more rigorous treatment, it is clearly necessary to start from the 
unscreened Coulomb interaction (2.20) and to transform it, together with 
all other terms of (2.8), under all canonical transformations that are to be 
applied. In a Bogoliubov transformation, the cancellation of all 
‘dangerous’ processes must include those processes arising from the 
Coulomb interaction. 

Applying the transformation (3.41) to H 


Coulomb gSIves 
2rre? 
A comomb = 2 Ow PG is eee) caw, 5% eT Cl kw ikaw.) 
k, i,w, 
8, ye 
X (Ujoy gf + VpH_y — 5°)(Myy sw iM ew, Bete ewitet wean eae, (09 10) 


The part of Ao. diagonal in the a*, a operators (i.e. the terms arising 
when k=I+w, s’=s) is, expressed in terms of «+, «a, 


H,,=2 > {u2v2+ (wj2—v7,2)(uP — v2 )oy gt o44277€?/|k — 1]2.Q 
k, I, s, 
PenicperOrder terms... « . . « (osz4) 


The first term in (3.71) may be grouped with .¥ (3.47) and the second with 
HH, (3.48). Both these expressions are divergent, but the second diverges: 
even in the ‘normal’ state. If, ad hoc, the w~? is replaced by a ‘screened’ 
expression, these terms converge; moreover the value of .Y does not play 
any important role, while the contribution to H, is not very sensitive to. 
small departures of u;, v;, from the ‘normal’ values. This term will not 
be of importance in studying the difference between the normal and 
superconducting states. It will, rather, represent a further renormalization 
of the (normal) electronic energy levels. A full treatment is difficult, 
requiring a sum of a perturbation series, to infinitely high order (Gell-Mann 
and Brueckner 1957, Bogoliubov et al. 1958 a). There are no ‘dangerous’ 
terms in (3.71), however, since the summation over the two possible spin 
assignments cancels the quadratic terms non-diagonal in «*, «. The only 
‘dangerous’ part of (3.70) creating two quasi-particles from the ground 


state, arises when k= —I, s= —s’; this assignment leads to no diagonal 
terms, 

2rre* 2 + + 3.72 

v= p Owe tktw Pi tetw (UI? = 07.7) (Oey Fp T+e.0). . (3.72) 


When the. ‘uspctoel terms of H, and H,, are included, eqns. (3.47) (3.48) 
and (3.49) become (where the superscript c denotes the inclusion of Coulomb 


he 


Se = 2 aA C)v,2+2 2, Upev 2. eee Siegeienle ee (3573) 
Hy’ = S Heo (On Oa &) at i Oks Foy st (€ Siam 
w k, s 
+ > (uw? —v;/*) - 2rre?/|k = IPQ} a0,” —0;21, . e A ° ° . (3.74) 
! 


cee > Ale nnb)U 0, + > (2rre?/|k — 17.2) (uw, — 94") ur} 
k 1 


(i ee C.C.) ame We NR CN, ye UE te als (3.75) 
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and H, and H, are unaffected. The (divergent) contribution to H,° 
will be assumed to have been included in the energies #,. The compensa- 
tion of ‘dangerous’ processes in second order gives, in place of (3.55), 


CDs 
i et Nery a 9\—3 d3w 4 vee La 
U7.V0 7. = (UW, (Eh Met hos arr WU VY ! 
Wied ( : ; ( | LS as EF, Eix+wt hwy, 


ks Re 


we / 


where in (3.76), all integrals converge, for both large and small w. The 
solutions for w,, 0; are still given by (3.57), provided I is redefined, 


D(|k-+ whe OD 0° =a 

¢) = (277) | ee )/ ao eC 2 

re)=(@ry | (Tek + w))+finwe) | Bet Bierwit hey 0 
st (Beta) 


With the same assumptions that were made in (3.60), the new equation 
(3.77) is solved by 

Na/2 ti di 
TAC tat 


where A is an appropriate average of 27e?/Qw? and is clearly of the order 
of magnitude He?/k,? (cf. (3.69)). The criterion for superconductivity is 
once more given by eqn. (3.68). 

The assumption of smallness of the Coulomb energy is not essential, but 
if it is not made, the perturbation series must be carried to infinite order. 
This, indeed, is what has to be done to avoid the divergences in (3.73), 
(3.74). Bogoliubov, Tolmachev and Shirkov have shown how this may 
be done, by the methods of Gell-Mann and Brueckner (1957) and Sawada 
(1957); the actual criterion for superconductivity is somewhat modified, 
and is more complicated than (3.68). 


(4) 


3.5. The Energy Gap at Finite Temperatures 

When trying to calculate thermodynamic properties of a superconductor 
we can no longer use Bogoliubov’s method, the compensation of ‘ dangerous ’ 
processes. At a finite temperature, real transitions can and do occur, 
and in any perturbation theory lead to vanishing denominators. They 
occur with so small a statistical weight, in an ensemble average, that a 
statistical mechanical perturbation theory for the free energy can be 
developed (Schafroth 1951, Nakajima 1955) and converges more strongly 
than the usual perturbation theory for the internal energy. No processes 
are ‘dangerous’ in Bogoliubow’s sense. 

The quasi-particle transformation can however still be used, with 
Valatin’s approach. The coefficients w,, v; of (3.41) are selected to make 
the free energy a minimum. Thus at 7'=0, we recover the results of 
$3.2. For simplicity, we consider the ‘truncated’ interaction 
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Hamiltonian (3.17), (3.31) of Bardeen et al. (1957)7. The coefficients 
Uz, V, Of the quasi-particle transformation (3.41) are unspecified, for 
the present, but we note that the states 


{i a4: b IV) 
an arbitrary set of states (k, s) 


form a complete set. Provided the u,, v, are correctly chosen (e.g. 
Bogoliubov’s choice as 7'->0), these multiple excitation states will be 
regarded as statistically independent. Here a self-consistent approxi- 
mation has been implied, since the best choice w,, v,, and therefore the 
set of ‘independent’ quasi-particles, will depend on temperature. 

It is convenient to denote the expectation number of quasi-particles 
in (k, s) 


ite a (Oks Mes )Bnsemble average® ow cme Sakata ee (3.79) 


Then the entropy is 
pee eve In Te st US, s) In (1—fxs)}> . . . (3.80) 


on account of the statistical independence of the quasi-particle excitations. 
The contribution of the state (k, s) to the expectation kinetic energy is 
v;,°n,, if there is no quasi-particle in the state. This follows because the 


wave function, II «,,|0), expressed in terms of a,,+, a,, operators, contains 
k,s . ce ad . 
the factor (u;,d_;, .+Vx%, 5"), 1.e. the probability of a real electron in 


(k, s) is v,?. If there 1s a quasi-particle in the state, the wave function 
contains instead the factor (v;,d_,,_,—U,4,, 5"); the probability for a real 
electroniswu,?. These two configurations occur respectively with statistical 
weights 1—f,,, fs. Hence the expectation kinetic energy is 


Donat, eerie): = $ : Q 3 (3.81) 
With the factorizable potential (3.31) the pair of states (k, s), (I, s) contributes 
to the interaction energy an expression which is the product of — #/Q& 
and two factors, one from each state. If there is no quasi-particle in 
(k, s), the factor for state (k, s) (as in § 3.2) is u,v; and if there is a quasi- 
particle, the factor is —u,v,. The weights of these states being (1—f,) 
and f,,, respectively, the interaction energy is 


— > (F/QE)(1— 2fgs)(1 — fis) Oneu?- Bes oe KER) 


k, 1,8 


+ The following derivation and minimization of the free energy, due to 
Valatin (1958), is exactly equivalent to that of Bardeen et al. apart from the 
use of the quasi-particle notation. Bogoliubov et al. (1958b) have rederived 
these results, using the quasi-particle transformation together with statistical 
mechanical perturbation theory. 
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The free energy is then 


F x > nr lntes 2 v7 ( 1 —fxs)} 
k,s 


— > (F/QF)u,v,4y(1 - 2fus)(1 — 2h.) 
ae ny {fs In Fist (1 —Fys) In (1 —Fius)}s mk ate (3.83) 
k,s 
which has to be minimized with respect to uw, and f,. 
If 1'(7) is defined by 
T(T)= > (P/QE)u,v, (1-2, 5), - + - + (3.84) 


k,s 


then the minimization with respect to u,, (or v;,) proceeds as in § 3.2, with 
17’) in place of [ and gives 


uP= Hl = galt T Te ee 
Substituting for w,, v, in (3.83), and minimizing with respect to f,, , leads to 
«T{ln f,,—In (1—fy, )}+V (n+ 02(7))=0, . . . (8.86) 

1.6. 
A= st exp (7 [n.24+ T(L)/eT)4, . . . . (3.87) 


an ideal Fermi gas distribution, with an energy spectrum 4/(7,?+ [*(7)). 
With the solutions for w,,, v;,, and f,, (3.84) describes an explicit temperature 
dependence of the energy gap 217’). As in the one-dimensional case, the 
gap vanishes at a temperature 7’), given by «Z7')=1-14ho exp (—1/F). 
An Ehrenfest (1933) second-order phase transition occurs there. The. 
specific heat calculated from the above free energy is in fair agreement 
with experiment. More direct evidence for the temperature-dependent 
energy gap is obtained from the experiments of, e.g., Biondi et al. (1958 a, b), 
Morse and Bohm (1957), and again agrees with the theoretical prediction. 


§ 4. ELECTRODYNAMIC PROPERTIES OF SUPERCONDUCTORS 


4.1. Supercurrents 


We have seen how the phenomenological energy gap model has been 
justified by the quasi-particle transformation. Many of the observed 
properties of superconductors follow in a natural way, in particular the 
thermal properties. The property of superconductivity itself is however 
more elusive. The mere existence of an energy gap in the ‘individual 
particle’ excitation spectrum is not sufficient}; indeed it is equally 
characteristic of insulators. 

The properties most characteristic of superconductivity are the vanishing 
of the electrical resistance and the Meissner effect. It is here that the 


+ Bardeen (1956) has attempted to show that the Meissner effect does indeed 
follow from the energy-gap model with some apparently plausible additional 
assumptions, but Buckingham (1957) has disproved the consistency of these 
assumptions. 
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microscopic theory is still incomplete. The absence of resistance is 
partially understood; a demonstration follows that at zero temperature 
there is no scattering of an electric current, that is to say the residual 
resistance is zero. The vanishing of the ideal resistance (when T#0) 
still remains unproved. 

There is still no clear physical picture of how the Meissner effect arises. 
Mathematical derivations have been given, but Schafroth (1951, 1958) 
has pointed out that almost any error will lead to a spurious Meissner 
effect. A derivation in which any approximations are made can, therefore, 
be misleading, and great care must be taken at every stage of the calculation. 
So far no derivation appears wholly convincing. 

It can be shown (Fréhlich 1954, Bardeen et al. 1957, Bogoliubov 
1958 a, b) that a superconductor has stationary states with a uniform 
macroscopic current ; these states have energy as near as we please to the 
ground state. The density of such states is very small, so that they have 
no appreciable influence on the thermal properties. 

The total momentum operator 


Reese oe ey Bo 2 SE 4i1) 
ks 


commutes with the Hamiltonian (2.8) and is therefore a constant of the 
motion. One can therefore look for solutions of the Schrédinger equation 
which have a specified total momentum different from zero—their ortho- 
gonality to the ground state is assured. . 

In one dimension, it is easy to see that the lowest state of momentum P 
is not one with an electron and hole excited, but rather a collective excitation. 
of the whole system. Let the whole electron system have a drift velocity v. 
Then relative to the electrons, the lattice modes will suffer a Doppler shift 
in frequency ; in place of (3.2), we have 


{hi(@y, — Wo? )Dw, Ow, FIN QD y,(Oy,* EXP (iwor) — C.0.) Xu, =A4 Xew,2 } 
{ow + Ul [ow oa ar SU SID (OE OX — iw yr) —CC)} yw», =A_ Xu, 
(4.2) 
and therefore in place of (3.3) 
By, =iN QDy, exp ( F iwyx)/h(wy, F Wor). eae). 
The calculation follows as in §3.1, and leads to a half-gap 
Wa 4teexpy (04/52) 3/22 ites ecnecen \a:4) 
and to a total energy 
W,=W+4N Q(m + 31?/FCos?)v? pa he (4.5y 
where W is the ground state energy (3.15) (Frohlich 1954); the effective 


mass is macroscopic. ; 
In a similar way, supercurrents can occur in the three-dimensional 


case. The method of Bogoliubov (1958 a, b) will be followed here; the 
momentum conservation is treated, like the particle conservation law, by 


D2 
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introducing a Lagrangian undetermined multiplier v. The Hamiltonian 


is taken to be 
H =H-C > at .—V- > hha, tm, . « - - (48) 
k,s k,s 


a 2 (€ al ¢ —hv. k)ays* Us a Hon “7 Aint se Ta snents 
8 
a. 2 (€x—mv 1a ymu* << C) ays” Us Si FB. . Fin is A ooutomo: - (4.7) 
8 

With the substitutions 

k’=k—mv 
(4.8) 

C= C+ gm? 

(4.7) reduces to (2.8) with the extra term —{./”. The lowest state with 
velocity v has energy }./ mv*; the collective motion behaves like a single 
degree of freedom of mass. fm. Since there is still an energy gap separat- 
ing such states from quasi-particle states, it follows that a supercurrent 
at absolute zero is metastable (cf. Landau’s (1941, 1947) theory of super- 
fluidity). It has still to be investigated whether at a finite temperature 
the supercurrent can lose its energy to ‘individual* modes of motion. In 
view of its high mass, it will only give up ~ 1/./ of its energy in a single 
collision, but it should be expected to have many more collisions than an 
individual particle. The metastability of persistent currents is equivalent 
to the vanishing of residual resistance, i.e. of elastic scattering of electrons by 
lattice defects. The still undecided question is whether the ideal resistance 
can also be shown to vanish. (In this connection, the experimental 
results of Galkin et al. (1957) should be mentioned—it appears that 
fluctuations can occur in which the resistance R approaches 10-4R,, 
where F,, is the normal state value.) 


4.2. The Meissner Effect 


In §4.1, it has been shown that persistent currents can exist in super- 
conductors (with cyclic boundary conditions) i.e. the vanishing of the 
resistance at zero temperature has been demonstrated. For a complete 
microscopic theory of superconductivity, it would still have to be shown 
that a magnetic field is excluded from the body—i.e. that the first of 
London’s eqns. (1.1) (or something similar to it) holds. Unfortunately 
no satisfactory demonstration has yet been advanced. 

It is useful to reformulate London’s phenomenological theory in terms 
of the Fourier components j(k), A(k) of the current and vector potential 
(Schafroth 1951). 

The general relation between current and vector potential, in an arbitrary 
gauge (uv, v=1, 2, 3) 

j,n= > | d°x'K (x—-x’)A,(x’) 2... (4,9) 


becomes, in Fourier space, 
jp) => Kun (p)Ay(p).5 ee See naa} 


v 


Theory of Superconductivity | 37 


Gauge invariance and continuity require 


> K,,2,= 2 PK ee Vee ers et L TL) 
Let us write : ; 
EP Py Py toa Pp, soe vt (4.12) 
In an isotropic medium, (4.11) can be satisfied only if 
(oF A nee ee eerwed Challis) | 
giving, from (4.10) 
1D Diz — pore? (DP) eee.) ee (4.14) 


where #(p) is the Fourier component of the magnetic field #=curlA. 
If w,(p?) is regular for small p, 


Wo(p")=Ay tap +. ee ° y +; és 5 5 (4.15) 
then (4.14) gives, in configuration space, 
eurl j= —a,V°# -—4,Vi#...., 1 eendt 6 (4216) 


which is just the usual expression for the diamagnetic current (with 
higher order corrections for field inhomogeneity). On the other hand, a 
simple pole in w,(p7) gives the Meissner effect, since if 


w,(p")=a_ p*+a)+a p+... te rine te elWe) 
then 
ip, (p)=—a_1,4'(p)—a p°*#(p)..- - - - (4.18) 
i.e. 
Curl jen, oa... aue oe wt estan. (419) 


which is London’s first eqn. (1.1) (apart from higher order corrections). 
Hence, provided gauge-invariance is accurately preserved throughout, to 
decide if a model exhibits a Meissner effect, it is sufficient to show that 
@,(p?) has a simple pole. Violation of gauge invariance, however, can 
introduce just suchapole ; itis therefore likely to lead to a spurious Meissner 
effect. 

The simplified Hamiltonian of Bardeen e¢ al. (3.17), (3.18) is not gauge- 
invariant, although Fréhlich’s original Hamiltonian (2.8) is. The re- 
normalization, if carried through exactly, cannot of course, spoil the gauge 
invariance. It is the simplification or ‘truncation’ which does so. The 
difficulty arises in writing down the Hamiltonian in the presence of an 
electromagnetic field. In the untransformed Hamiltonian, no difficulties 
arise, since it contains only a term quadratic in the electron momentum 
operator p;=ihV,, and terms independent of p,, which commute with x,, 


FLEey pet 2a V (See K periee. Ow ky ads Oy, ols)s eee oy ens a2 20) 


Hence, in the presence of a field, we need merely to replace p by p—eA/c, 
and to write for the current density operator 


j(x) = (e/2m) ) > {pdx x) + 8(x;— x)p; — 2A(x,)(e/c)5(x;—x)} (4.21) 
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(again obtained by replacing p by p—eA/c in the expression for the current 
with no magnetic field). 

If however, in (4.20), V is a function of the momenta too (as it is in 
Bardeen et al.), then (4.21) does not automatically hold (Schafroth 1958), 


since then 
[eV] s# 0. ate ee (4.22) 


Bardeen et al. nevertheless assume (4.21) and Anderson (1958) has shown 
that this choice leads to an approximately gauge invariant expression 
for the current density in the ground state. That this is too weak a 
condition has been demonstrated by Schafroth—he shows that even if 
1(p) + w,(p?) = 0 (4.13), and w, w, are singular, so that 


lim p?ar,(p?)-—p*w,(p?)40 . . . . . (4.23) 
p>0 ¢ 
(respectively the ‘longitudinal’ and ‘transverse’ sum rules of Anderson 
(1958)), the arbitrariness in the choice (4.21) prevents one from drawing 
any conclusions about a Meissner effect. 

Since there is no unique procedure for defining the current operator, 
to the postulated current density (4.21) one may add a term 


i(x)= —(AJN Q)A(x) 2 2 2 ee 
where A is arbitrary, and where the only restriction on the gauge is 


divA=0.<>) 4. 215) x ee 
In Fourier space, 
t(P)= —A > (2, P,—8,,2)4,(p)/p? . . . . (4.26) 
so that if j(x) is gauge-invariant, so is j(x)+i(x). The addition of j to j 
will however alter (4.23), 


p*{wy( p”) —w,( p®)} > p*{ay( p?)—w,( p?)}- 2A... (4.27) 


Since A is arbitrary, the magnitude of the Meissner effect is arbitrary and 
in particular zero for a suitable choice of A. The good agreement of the 
predictions of Bardeen et al. with experiment, for these properties of super- 
conductors, thus appears fortuitous. 

The difficulty arises, of course, because there is no unique choice for the 
current density operator, when the Hamiltonian is not gauge invariant. 
If instead of the ‘truncated’ Hamiltonian one used the full renormalized 
Hamiltonian, then the expression for the current density would be unique— 
it would be got by transforming (4.21), under the unitary transformation 
(2.29). Unfortunately one cannot find a further unitary transformation 
leading from (3.19) to (3.31), i.e. a transformation which would ‘truncate’ 
the interaction Hamiltonian. 

A valid procedure would be, therefore, to start with the Hamiltonian 
(2.8), which has a unique current density (4.21) associated with it and to 
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subject both to the Bogoliubov transformation. Attempts along these lines 
have been made by Rickayzen (1958) and Wentzel (1958). Rickayzen, 
after performing the Bogoliubov transformation, makes a non-gauge- 
invariant approximation, and it appears that one could introduce 
Schafroth’s arbitrary Meissner effect (4.24) there at a later stage in the 
calculation than Anderson or Bardeen et al. Wentzel’s calculation is gauge- 
invariant at all stages, but he does make other approximations—namely a 
perturbation theoretical expansion of the magnetic terms in the 
Hamiltonian. It is not yet clear whether this type of approximation can 
also lead to a spurious Meissner effect (Schafroth, private communication). 

The need for a rigorously gauge-invariant formulation arises, of course, 
because we have not a clear physical picture of how the Meissner effect 
‘works’. If the programme of Schafroth e¢ al. (1957) could be carried 
through and the condensate really shown to have the properties of a 
condensed ideal Bose gas, then the Meissner effect would follow qualitatively 
already. Approximations could then be safely made, without introducing 
a spurious Meissner effect. But in the absence of a qualitative under- 
standing Schafroth’s work shows that mathematical approximations can 
be very misleading. 


4.3. Paramagnetism 


The Pauli spin paramagnetism of an electron gas is masked by the 
Meissner diamagnetism of a superconductor, and is not directly observable 
by static methods. However, Reif (1957) has observed the Knight shift 
(i.e. the shift of the nuclear magnetic resonance frequency in the medium 
compared with the free atom (Knight 1949)), which gives a measure of 
the paramagnetic susceptibility of the electrons—Reif’s experimental 
result is that the Knight shift in colloidal superconducting mercury is 
1:6%, about 2/3 of the value (2-5%) for the normal metal. (The experi- 
mental situation is not completely clear, however. Knight ef al. (1956) 
find a Knight shift of only 0-5%.) 

As the theory of superconductivity stands at present, it is incapable of 
accounting for this Knight shift. The quasi-particle excitation spectrum 
with an energy gap will give zero susceptibility for fields less than ~ 21T'/u 
(where » is the Bohr magneton), and therefore zero Knight shift. To 
account for the Knight shift, one requires the existence of excitations with 
small energy and with a net magnetic moment. Provided the density of 
such states is very small, the specific heat contribution can be small enough 
to preserve agreement with experiment. 

Heine and Pippard (1958) have suggested that the difficulty arises 
from the precision of the correlations between (k, t ; —k, ¥ ) electron 
pairs. If, in a fairly phenomenological way, the correlation is assumed 
less sharp, they are able to get a Knight shift of the correct order of 
magnitude. However, it is just the sharpness of the correlation that endows 
the ground state with its ‘ superfluid’ properties ($4.1). It is not obvious 
that the correlations can be ‘smeared out’ in a consistent way. 
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In conclusion, it appears that the nature of the ° condensation’ is now 
well understood, arising from the electron-phonon interaction. All 
properties which depend mainly on the density of states or energy gap 
follow from the theory, and agreement between theoretical prediction and 
experiment leaves little to be desired. Thus e.g. the specific heat, thermal 
conductivity and threshold magnetic field agree well (the latter 1s not 
primarily an, electromagnetic property, since it is derived from energetic 
considerations—H,2/87 is equated to the ‘condensation’ energy per unit 
volume). 

The electromagnetic properties are less well understood: although the 
ground state has been shown to have some ‘superfluid’ properties (§ 4.1), 
the Meissner effect has not been demonstrated in a wholly convincing way, 
nor the Knight shift. There is not a sufficiently clear physical picture of 
the process to give one confidence in the approximations, particularly as 
the subject has a long history of spurious Meissner effects arising from 


approximations. 
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§ 1. IytRopucTION 


THE phenomenon of superconductivity was discovered by Kamerling - 
Onnes in 1911. Since then, the subject has been studied intensively 
during nearly fifty years and a considerable amount of experimental data 
has been amassed. However, this has not led to any material progress 
in the understanding of the physical principles underlying the phenomenon. 
The first step on the road to its understanding was made in connection 
with the discovery of superfluidity and the construction of a theory for 
it (Landau 1941). Landau also drew attention to the analogy between 
' superfluidity and superconductivity, which allows superconductivity to be 
considered as the superfluidity of the electron gas in the metal. 

The second important step was the discovery (Allen ef al. 1950a) of 
the isotope effect, i.e. the dependence of the critical temperature 7’, on 
the isotopic mass of the superconductor. This effect directly showed the 
importance of the lattice in the phenomenon of superconductivity and 
in this connection Fréhlich (1950) and Bardeen (1950) pointed out the 
possibility of a new type of interaction between electrons by the exchange 
of phonons. However, until very recently this approach did not give 
satisfactory results and the cause of superconductivity remained essentially 
unknown. On the other hand, it was clear to all (Feynman 1957) that 
there were sufficient experimental data to construct a theory and that all 
that was lacking was some fundamental idea which would enable the 
orthodox methods of quantum mechanics to explain the whole phenomenon. 

This missing idea is contained in a short letter by Cooper (1956). Cooper 
showed that for a pair of particles which obey Fermi statistics and which 
have a mutual attraction, no matter how weak, there will exist a bound 
state. This extremely important result is the key to the understanding 
of the phenomenon of superconductivity. In the ground state of a metal 
without taking account of the Cooper effect the electrons fill (in the 
isotropic case) the Fermi sphere in momentum space. If however we 
assume that there is in metals an effective attraction between the 
electrons, pairing will occur and the ground state will then be lower than 
for free electrons by the binding energy of these pairs. The electron pairs 
have integral spin and therefore obey Bose statistics. Now it is known 
that a Bose gas is superfluid at absolute zero and applied to a Bose gas of 
charged particles this property manifests itself as superconductivity. The 
above considerations do not pretend to be rigorous but they do 
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indicate that superconductivity can be completely explained on the 
basis of the Cooper effect. 

A complete and consistent theory of superconductivity was given by 
Bardeen et al. (1957, referred to here as BCS 1957). They put forward a 
technique which enabled them to examine the ground state of a super- 
conductor (i.e. that at 7’=0) and to obtain the excitation spectrum. 
On this basis they calculated the thermodynamic properties and deduced 
the electrodynamic equations for superconductors in weak fields. In the 
field of thermodynamics the fact that a finite jump is obtained in the 
specific heat at the transition between the normal and superconducting 
states deserves mention. For the description of the behaviour of super- 
conductors in a magnetic field, a quantity characterizing the size of the 
electron pair turned out to be an important parameter. Ifthe penetration 
depth of the field is large compared to this parameter, there is a local 
relation between the current and the field potential and the old London 
equations are then valid. If however the penetration depth is small 
compared to the above parameter, there is no local relation between the 
current and the potential and the properties of the superconductor are 
then described by other equations, obtained earlier by Pippard (1953) 
on, semi-intuitive grounds. ; 

Bogolyubov (1958, see also Bogolyubov et al. 1958) put forward several 
original methods which lead to the results of BCS (1957), and one of these 
is very similar to the technique used by BCS. Another, which seems to 
us more natural, is similar to the method which he used in the study of the 
excitation spectrum of a Bose gas. 

The authors of the present article, using Bogolyubov’s technique, 
have carried out calculations analogous to those of BCS (1957), and an 
exposition of the modern theory of superconductivity is given on the 
basis of the Cooper effect and the above approach. Some results are 
also given which are not contained in the work of BCS. In particular 
the réle of the electronic Coulomb interaction was not explained by 
these authors. It should be pointed out that that there is as yet no rigorous 
solution to this problem, though there is an interesting paper by Tolmachev, 
(see Bogolyubov ef al. 1958) which showed that the effect of Coulomb 
repulsion is weakened in the presence of electronic interaction through a 
phonon field. The results obtained by L. P. Gor’kov and the authors of 
the present article on the behaviour of superconductors in high frequency 
electromagnetic fields are also given on the basis of the new theory 
(Abrikosov et al. 1958). 


§ 2. Tur Cooper Errect 
To understand the basic content of Cooper’s idea, we consider the 
following example: a gas of Fermi particles with a direct interaction 
between particles with energies ¢(k) within a small energy range w in the 
neighbourhood of the Fermi sphere —w< «(k)—e¢(k,)<w and constant 
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inside this region} (Abrikosov and Khalatnikov 1957). The Hamiltonian 
of such an interaction will be of the form: 


=9 Lat ky’, s)at(ky’, —})a(ko, 2)U(ky, —$) ([E(k)| <e). (2.1) 
where the summation is over values of k, etc. such that 
k,+k,=k,’+k,’. 


The operators a+(k, c) and a(k, c) correspond to the creation and annihilation 
of particles with given momentum and direction of spin; €(k) is the energy 
of the particle, reckoned from the Fermi surface, and the summation in 
(2.1) is made with due account of conservation of momentum. We shall 
assume that the constant g is small enough to permit the application of 
perturbation theory. 

We consider the scattering of two particles with nearly opposite 
momenta, i.e. such that |k,+k,|<k). To a first approximation the 
scattering amplitude is simply equal to g if the particles have opposite 
spins and zero if the spins are parallel. This is of course also true for an 
arbitrary angle between k, and k,, but the next approximation already 
shows the singularity of the case examined. Although the second approxi- 
mation gives only a small correction to the scattering amplitude for an 
arbitrary angle between k, and k,, the correction may increase without 
limit for almost opposite momenta. Thus let us consider the terms of 
the transition amplitude corresponding to the transition of two particles 
outside the Fermi sphere (€(k) > 0) without participation of the particles 
inside the Fermi sphere. Such a term arises only for particles of opposite. 
spins and takes the form (setting #=1 in this and subsequent formulae) 


= Age. | d3q Pap tes Ayia | aset EM se 5i°6. 
(27)? E(heq) + & (Fea) — &(Ky + Ka — q) — €() 
where the integration is in the regions 


0<&(q)<w, 0<&(k,+k,—q)<w. 


(2.2) 


In this formula, since all the momenta & are near to ky, we may put 
E(k) =v(k—ky), v being the velocity of a particle at the Fermi surface. 
If we also introduce the notation A=£,+€, and p=v|k,+k,|, then since 
and p» are small we get 


Qrgrke? [° 1 1 
- ae | aca) | OW ore penn menriety a(S) 


This integral is logarithmic in €(q). For the case A<w and «<w it is not 
difficult to see that to the first order this expression is equal to 


gKo?, 20 Ss ITO) 
4n2y = max {A, wu} 
Sen eee ee 6 eee 


+ Later we shall see the relation between this and the real interaction between. 
electrons in a metal. 
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It is easy to see that such a quantity contains a term which corresponds to 
a transition of two other particles which are inside the Fermi surface 
to a final state k,’k,’ and the subsequent transition of the initial particles to 
the positions ‘til ade vacant. If we consider transitions which concern 
one of the initial or final particles and one of the particles from the Fermi 
surface the logarithmic factor is absent. (Only this type of virtual 
transition may take place for parallel spins.) 

Consequently, with an accuracy including small terms, the scattering 
amplitude for particles with antiparallel spins and almost opposite 
momenta is to the second order of the form 


g? ko” 2w 


————. 2.5 
Qn "max {Ap i 


A@S9 4 
This expression is however, not correct when p» and A are very small, in 
which case the second term of formula (2.5) may become larger than the 
first, and higher terms of the series given by perturbation theory must 
be included. In this it is sufficient to limit ourselves to those which for 
each power of g contain the corresponding power of the large logarithm. 
It is easy to see that these terms arise for such virtual processes in which 
the initial pair is not broken up, i.e. both particles simultaneously transform 
either inside or outside the Fermi surface. The series which arises as a 
result of this is a geometric progression, the first two terms of which are 
given by formula (2.5). The sum of the whole series is evidently 


gk gko? 2w 
y. 1— : oes sa eee 
= a/| Inn 2 EEE ae 


In the Hamiltonian H’ the constant g can be either negative or positive. 
In the first case (repulsion) the scattering amplitude tends to zero for 
small A and p. However, in the second case (attraction) the position is 
completely different and in fact, for sufficiently small A and » the scattering 
amplitude has a pole, which corresponds to the appearance of a bound 
state. We note that the formation of such a state would not be possible 
without the presence of the other particles, which, together with those 
considered, are in the ground state and form the Fermi sphere. The 
quantity A corresponds to the binding energy of this state. Different 
values of A will correspond to different values of j, i.e. different angles 
between k, and ky, but the calculation of the dependence is beyond the 
limits of our accuracy. From considerations of symmetry it is natural 
to conclude that the ground state is that in which k, and k, are opposite. 
In this case the value of A corresponding to the Pele ected by Aj) will 
be equal to 


2ar2y 
Ay=2wexp{ — = ooh SE oh, ORT EET 


Thus we arrive at the conclusion that for any attraction between the 
particles in a Fermi gas, no matter how small, there will arise bound states 
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_ of pairs of particles with opposite momenta and spins. This is the essence 
of the so-called Cooper phenomenon. Cooper himself obtained it by 
different and simplified arguments which, it should be noted, lead to 
quantitatively incorrect results since they do not take into account virtual 
transitions of the initial pair with other particles of the Fermi sphere 
(the second of the processes examined here). 


§ 3. THE GROUND STATE OF SUPERCONDUCTORS 


The existence of the isotope effect in superconductors is evidence for 
the important réle played by the lattice in this phenomenon. The 
influence of the lattice shows itself in an electronic interaction through 
the phonon field and this interaction has the form of an attraction; at 
the same time there is also the Coulomb repulsion between the electrons. 
The electronic interaction through the phonon field can be described using 
the Fréhlich Hamiltonian, though a consistent account of the Coulomb 
interaction is as yet an unsolved problem. The difficulty in taking account 
of the Coulomb interaction may be avoided by the introduction of a certain 
effective direct interaction between the electrons. We will consider that 
this direct four-fermion interaction leads to a certain effective attraction, 
which is represented by a small binding constant g. Such an attraction 
only acts in the electrons in the neighbourhood of the Fermi surface, in 
a certain narrow momentum range the width of which is of the order 
of the maximum Debye frequency. Ifthere were no Coulomb interaction 
this could be considered rigorous, since the phonon interaction concerns 
only electrons in this momentum band}. 

A direct four-fermion interaction is described by the Bardeen 
Hamiltonian 


H= » & = 1) at(k, o)a(k, o) 


+g dat(hy', d)at(ky’, —4)a(ky, $)a(ke—3) . - . - (3.1) 
where the second summation is over values 
k, —k,’=k,’—k, 40, 


and where »=k,2/2m and ky is the momentum at the Fermi surface. In 
order that in further calculations we should not need to consider 
conservation of the number of particles in the system, the term 


—pN = —p Yat(k, o)a(k, o) 
has been induced in the Hamiltonian (3.1). The difference 
k?/2m — w= 0(k — ko) re ie eae ae LO ]9")) 
(for values of k near to ky) will in future be denoted by the letter ¢,. The 
condition for the momenta in the second sum in (3.1) is due to the fact 


that, as the energy of the phonon tends to zero, so does the phonon inter- 
action. Instead of the operators at(k,o), a(k,o) we introduce new ones 


+ The Coulomb interaction will be dealt with in more detail below. 


P-M.S. 
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a(k, 0), «(k, 1) and defined by the equations 
a(k, 4) =u,,0(k, 1) +0,0+(—k, 0) 

a(k, th 2) = Uk, Oo} Vat ( ate k, 0) } 

where u,, and v;, satisfy the condition 
tyt+vugali gcc ae oe ee 

When these conditions are obeyed the commutation relation for the new 
operators a(k) will be the same as that for the previous ones a(k). According 
to Cooper (1956)the Hamiltonian (3.1) led to the possibility of the formation 
of bound electrons. The operators a(k) characterize the reconstructed 
system in which electron pairing has taken place. In the ground state 
of the new system the creation of pairs with opposite momenta and 
spins must be forbidden, although this was possible in the original system. 
Therefore in the Hamiltonian (3.1) we express the operators a(k) in terms 
of a(k) and obtain the conditions which forbid the creation of pairs of 
excitations a(k,1) and a(—k,0). Hence the Hamiltonian (3.1) acquires 
the following form 


(3.3) 


H=H,+H,+H, 
where 
Hy=2 > 042+ > (uj,2—0;,2){at(k, la(k, 1) + «t(k, 0)a(k, 0)}, 
k k 


H, =2 DE urearrctor* 1)at(—k,0)+a(—k, O)a(k, 1)}, 


Hy=g > (u(ky' u(ke’)at(hy’, Lat (he’, 0) — v(ky’)o(ke’) 
x o( — key’, O)ou( — keg’, 1) — u(y’ )v(hee’ art (hey’, Lor — ha’, 1) 
+ v(ke’)u(Keg’ )ou( — ky’, O)a*(ky’, 0)}{(w(ky)u(he)a(k, 0)x(k a, 0) 
— (ky )o(ky)ort( — ky, L)at(— ka, 1) — w(hy)v(hq)a(hy, L)at( — ke, 1) 
+ v(ky)u(ky)at(— ky, 0)ax(ko, 0)}. Me rrr ee 4 
Here u(k) is written for u, and the summation in H, is over 
k, —k,’=k,’—k, 40. 

Let us now consider a superconductor at a 7’'=0 and let us calculate the 
properties of the system in the ground state. The operator at(k), acting 
on a vacuum, will create a quasi-particle, but the action of the operator 
a(k) on a vacuum gives zero, since at 7'=0 all the states are full. The 
condition forbidding the creation of pairs is obtained by equating to zero 
the coefficients of all the terms which contain the product of the two 
operators at(k,), «t(—ky). Thus we get 

2E UV =G(Uy2 — 0,2) > UpYyr =. « ss (3.6) 
. k 
We introduce the notation 
A(0)}= 9.) Unie cur) pena (oenk 
k 
and substituting for w;, and v;, in (3.6) and (3.7) in terms of &, and A(0), 
we have 


my? = ELL + Ep] {A(0)2 +82}, . eS Berta, 
v2 =4[1 —E,//{A(0P2+E,2}]. ee ee ay Cy 
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Substituting these expressions for wu, and v, in (3.6) we obtain an equation 
which defines the constant A, namely 


a7 DW AAO) ee ee on  . $(8.10) 
After transforming the sum in (3.10) to an integral and integrating we get 
Pe ovine 2p) eee ee Soe (3.11) 


where v=k,?/7*v and (3.11) defines the constant A(0). As we shall see 
below the constant A(0) plays the part of an energy gap in the spectrum of 
excitations ; from (3.11) its magnitude is 


DOy=2aexp(—2/gv) % 5.3. 2 « % (3.12) 


which is the same as Cooper’s binding energy (2.7) of an electron pair. 
From the form of the expression it is easily seen that no perturbation 
theory could produce this effect, since the function exp (—1/x) cannot be 
expanded in a power series in the neighbourhood of x =0. 

The energy of the ground state is obtained by the diagonalization of the 
Hamiltonian (3.5) 


SY ATES DIAN ge SE BIEY 
Varying this expression with respect to n and subject to the supplemen- 
tary condition (3.4) we obtain (3.6) as an extremum condition. Thus 
equation (3.6) is simultaneously both the condition forbidding the creation 
of pairs and the condition for an extremum of the energy. 
Using eqns. (3.8) and (3.9) for w, and v;, we write the energy of the ground 
state in the form 


By=2 D Su BAO 2 UY viA(O+ £4'}- - . » (3.14) 


The first summation is taken only over negative &, (inside the Fermi 
surface), since v;,, tends to 0 very quickly outside the Fermi sphere (€;, > 0). 
From formula (3.14) we see that because of the second negative term the 
energy obtained for the ground state is lower than the zero-point energy 
corresponding to a full Fermi sphere for free electrons. We note that 
(3.11) has a solution only for g > 0, i.e. for an attraction between the particles. 


3.1. The Réle of the Coulomb Interaction between Electrons 
(Bogolyubov et al. 1958) 


In order that superconductivity should be possible it is necessary that the 
total effect of the phonon and Coulomb interactions should be attractive, 
leading to the formation of electron pairs; we shall therefore consider the 
role of the Coulomb interaction in more detail. Since the problem of an 
exact treatment of the Coulomb interaction has not yet been solved, we 
shall try to take account of it by means of a mode (Abrikosov et al. 1958). 


E2 
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The Hamiltonian of the initial system previously included a certain 
effective attraction which was effective in the momentum range —w to w, 
where w corresponds to the limiting Debye frequency of the lattice. This 
interaction model was equivalent to an interaction through phonons. 
The Coulomb interaction can also be accounted for approximately by a 
model if the Hamiltonian includes a constant repulsive interaction which 
acts over the momentum range —@, to w,, which is larger than the first 
rangey. Thus the initial Hamiltonian now has the form: 


H= atk, o)a(k,o)+ ¥ glk, ky) a*(k, 3) a* (ky, — 3) a(k, 3)a (hy, — 2) 
k 


k, ky 
where 
_{g+9v  |E(k)|, E(k] <@ 
aks) =] ts lé(k) |, [eC ky) | < a (3.15) 
and g>0, 9, <9. 


The further calculations are carried out as before, i.e. the operators a(k, 1) 
and a(k, 0) are introduced as in (3.3). The Hamiltonian (3.15) is diagona- 
lized and the transformation coefficients uw, and v; are found from the 
condition for an energy minimum. We then find 


Ue =F +E V (Ee +A), « - - + ~ (3.16) 
oP =F l-EnV (EP + Aw). . . « «+ (3-17) 

The quantity A, is now a function of k, defined by the equation 
A,= 2X gk, ky) Uy, Vx,+ A Eero ees 


Using this equation we get two equations which define the quantity A 
in the two regions of momentum space : 


A -{. |E(k)|<w, 
= 
yma |E(*)| <4, 
Ao=3(9 +91) Ap In (2@/ Ap) + $91 Ay In (,/e), 


A, =39, A ln (2w/A,) +49, A, In (w,/w). . . . (3.19) 
Eliminating A, from these equations we get an equation for A, : 

1 In (2w/A moe td eel Ge 

3 In (2«/Ao) Gi eer Ta; | 1, sm cee ean 


This equation has a solution only when the terms in the square brackets 
exceed unity and this condition is the criterion for the possibility of a 
bound state, and consequently for the occurrence of superconductivity. 
It is interesting that the interaction g(’, k,) may be repulsive over the whole 
momentum region (g+g,<0) but that nevertheless the condition for 
superconductivity may be satisfied. This is related to the fact that the 
Coulomb interaction which acts over a wider momentum range is renorma- 
lized, and the renormalization contains the large logarithm In (w,/w) 


Oe  —— ————————SSSSFSSSSSSFSFFFfsmsmsMHshFeFeseseseF 


t This corresponds to the fact that the Coulomb interaction is screened at 
distances of the order of the lattice spacing. 
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§ 4. Tue THERMODYNAMICS OF SUPERCONDUCTORS 


We will now examine the properties of superconductors at temperatures 
other than the absolute zero. At non-zero temperatures systems of 
electrons are not in the ground state, and consequently the average values 
of the operator products will no longer be zero (as for 7'=0), but will equal 
the corresponding occupation numbers n,. 

We define the operators u;, and v,, analogously to the previous case, so 
that in the Hamiltonian all terms which lead to the creation of pairs should 
cancel out and the corresponding condition has the form: 


285, Un V7, = G (Uj? — V;,”) Pai Ae — NK — Np). - + (4.1) 
Introducing the notation 
A=g > Uy, 0; (1 — Mp1 — Mo) ° ° . ° ° (4.2) 


we express uw, and v, in terms of €, and A and obtain the following 
expression for A: 


ea (l= tye —Me)i4/ (Ae, )- = (4.4) 
For u,, and v;, we use the earlier formulae (3.8), (3.9) : 
Uy, OF V;,7 = F{1 + E,/4/(A? + &;,2)}. a) eee (as0) 


We now calculate the energy of the system in the excited state as a function 
of the distribution of excitations. Taking the diagonal matrix element of 
the Hamiltonian (3.5) we get 


H=2 Dene + DE (63.2 — V4?) (My + Nxo) 
ee, 2 ME Vy, Wye Vie (1 — pq — Myo) (1 — Ny Ny) cee (4-0) 


The system of excitations can also be considered as a system of quasi- 
independent particles, so that we may use the well-known formula for the 
entropy: 
S= —} [Mp91n Myo + (1 — Mo) Mn (1 = Mp9) + My 1 pq + (L — Mya) In (1 — 4). 
i aay: (4.7) 
From the condition for an extremum of the entropy for a given energy we 
may derive, by a method analogous to that used in the theory of a Fermi 
liquid by Landau (1956), the distribution functions nj, and jo. They 
turn out to be the Fermi function 


1, =lexp hee) t+ lye) mee. ey (4.8) 
where «(k) is the derivative of the function 
AV9) AOI 1 Bed Ba be ab ice ME) 


The variation of the energy EH with respect to n, may be carried out at 
constant u, and v,. The fact is that the condition (4.4) may be derived 
from (4.6) as a condition for an extremum of the energy with respect to 
u,; then, varying EZ with respect to w;, and using the condition w,,2%+0,2=1, 
we easily obtain 


ae {QE Up Vp — (Uy? — 1° )(1 — Mea — Mr0)9 DY eK Pn (Ler — Mire) } = 0 
Vv Kk’ 
ke (4.10) 
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which evidently coincides with (4.1). Taking this into account, we find 
for the excitation energy «(k) : 
c(h) = PE = 8 (2-018) +2 Oe tee (Meo Ms) (41) 
k Ke 

and from (4.1) and the (4.5) we finally obtain 
e(k)=4/(A® 42). % oi 3 ee ees 
In form this is similar to the excitation energy at 7’=0, but it must not 
be forgotten that the quantity A, which plays the role of an energy gap, 


depends on the temperature and is given by eqn. (4.4). Substituting 
(4.8) in (4.4) we get a relation which determines the magnitude of the 


energy gap: 


t 


kg Y{tanh(€/2T)/e}=1.  . . . . . (4,18) 
k 


4.1. The Temperature Variation of the Energy Gap 
We now consider eqn. (4.13) in more detail. For 7’=0 the tanh tends to 
unity and integrating then with respect to dé between the limits —w and w 
we get 
tqvin(2m/Ajj=De Se sa eee ee 
where v=k,?/7%v. This is the condition (3.11), which was obtained 
earlier as determining Aj. Equation (4.13) can now be written as 


” tanh (¢/2k7') 2w ® tanh (€/2k7')—1 
In(2a5] Ay) se ee ee spams Boek SB | 
a (old) | eon, aon ee een 
SO ae 
In the last integral the upper limit may be replaced by oo. Introducing 
a new variable by the substitution «=A cosh ¢ we get 


Ao _ dg 
In5= eaeccrvarst oe ee oe (SEG 


From this formula we may find A over the whole temperature range. 
To do this it is convenient to express the right-hand side as a series of Bessel 
functions ky. 


in(Ge)=2 S (1 (FP). th ah 6a (Sik) 


Near 7'=0 we may take just the first term of this series and use the asymp- 
totic expansion of the function ky. This gives 


A=A,—/(2rTAy)exp(—A/T). . . . (4.18) 


The transition to the normal state corresponds to the vanishing of A and 
therefore (4.15) may be used to find a relation between A, and the transition 
temperature 7’,.. For small values of A/7 it is convenient to use the formula 
for the summation of the series (4.17) : 


A al ca 1 1 
In—2 =In—- + 27 1) ee eee eee 
n A Syl z 2 lo 1)a /{(A/T)? + (21— aad . (4.19) 
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sea Iny is the Euler constant 0-577. Hence at the transition point we 
ave 


Vea el (Gg ef . ~ (4:20) 
Expanding (4.19) in powers of A/’, we find 
In(T/T,) = — 5 00)(q) PEO (ZY. aero) 
Fig. 1 
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To the first order in (7'— 7',) 
A/T -33-06(1-—T)/T 2a... 1 (4.23) 


The function A(7') may then be calculated in detail from (4.21) and (4.19) 
and is shown graphically in fig. 1. 

Tn (4.14) the coefficient w is of the order of the maximum Debye frequency, 
which is inversely proportional to the square root of the atomic weight (JZ). 
Thus the energy gap in the excitation spectrum is inversely proportional to 
4/M, and since the energy gap differs from the transition temperature only 
by a constant factor we therefore have 


P~i/M, 


i.e. the transition temperature 7', is inversely proportional to the square 
root of the atomic mass (isotope effect). We note also a possible relation 
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between the energy gap and the Debye temperature, namely, as follows 
from (4.14): 
A,/w = 2 exp (— 2/gv). 
An investigation of the stability of the lattice (Tiablikov and Tolmachev 
1958) leads to a criterion which limits the exponential by the following 


inequality : 
2/gv> 2. 


Thus the energy gap must be about an order of magnitude smaller than the 
Debye temperature. 


4.2. The Specific Heat 
We begin with a calculation of the entropy. From (4.7), on integration 
by parts and the use of (4.8), we find : 
S=- 7\, 4/(e — A) Mae, « “we Gri ee 
and substituting «=A cosh ¢ this becomes 


ie Seah 12 K, (Gm). gike woeteeal 


m=1 
For 7'=0 we get the asymptotic value: 
S=4/(27A3/T)vexp(—Agil)- 5 eee eee 
Near the transition point 7’, the expansion of formula (4.23) leads to 
3 S/N 3q /A\4 ‘) 
Se OT a eet 
3 | 27? (7) cs 87? @ | 
where <i, . (420) 
q=1| x—* tanh? x daz ~ 0:8. 
0 


To calculate the specific heat, we differentiate (4.25) and, using the 
expansion (4.21), we find: 
T)/O,(T) 244-753 (TTT ae eee 

Here we have expressed the specific heat in terms of the value for the 
normal state at 7'=T7',, equal to C,(7',.)=47?/v 7... From (4.27) it follows 
that at 7'=7', the specific heat has a finite jump and that its derivative to 
the left of this point remains finite. A more detailed analysis shows that 
the next term in the expansion in 7'— 7’, is proportional to (7',—7')32, 

Differentiating (4.25), we obtain an asymptotic expression for the 
specific heat for 7’>0 


3 /2/A,\32 A 
OLE | 7? yet “Ss 30 = Ape 0 
(T)/C,(T,) ; J=( 7) exp ( 7). ees) 


The full curve for the specific heat is shown in fig. 2. 

In spite of the fact that the model considered, which neglects anisotropy, 
is a fairly rough approximation, it is interesting to compare the theoretical 
results obtained with the experimental data for a number of superconduc- 
tors, as is done in fig. 3. The data for Sn, Zn and Al were obtained by 
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Zavavitskii (1 958), and those for V and Nb are taken from Corak et al. (1956). 
From this diagram it is clear that in most cases the experimental data are 
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quite well fitted by the theory. The worst fit is for zinc in the low tempera- 
ture region, and this could be due to the fact that in this region the specific 
heat is very sensitive to the value of Aj. Since the magnitude of the jump 
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in the specific heat does not show up well on the logarithmic scale; we give 
the experimental values of C,(7',)/C,,(7".) separately in the table. 


Element CAT ies ta 
aH 2.82 
Al 2-5 
Zn 2-32 
V 2-52 
Nb 3-07 


It is evident that, despite the reasonable roughness of the model, the 
values which it predicts agree quite well with the experimental data. 

Here it is apposite to note a characteristic peculiarity of the real aniso- 
tropic case, in which the energy gap will in general depend on direction ; 
however at 7’, the separation must tend to zero simultaneously for all 
directions. In the measurement of the specific heat in the low temperature 
region it is natural that the minimum gap should be dominant. On the 
other hand the value of Ay, found from (4.20), which is valid for the 
isotropic case, gives some average value of the gap. Hence it is clear that 
the theoretical curve for the isotropic model will give a reduced value of 
specific heat in the low temperature region, and this is indeed observed 
(see fig. 3). 

4.3. The Critical Magnetic Field 

The magnitude of the critical field can be found by using the well- 

known thermodynamic relation 


H2/8a=F.— fF, nist ade eee ees 
and we now calculate the free energy F’, in the superconducting state, equal 
to F,=E-TS. se 4 seen SEOs 


The energy # is given by formula (4.6), and the entropy by (4.7), but 
before substituting the expressions for H and S in (4.30) we change their 
form slightly. On substitution of the expressions (4.5) for w;, and v;, and 
on simplification using (4.3), the expression for the energy becomes 


E= dE — A) (Et + AP) 4- GAME Aa) ey 
ae DA2mV/ (€,2 + A?) — A2n,(&,2+ A?)-}, 2... (4.81) 


After a simple integration by parts the entropy (4.7) can be written ast 


2 
ee 7 2 [21,4/(E,2 + A?) — A2n,/4/(A2+2)].. . . (4.32) 
Substituting (4.31) and (4.32) in the expression for F’,, we find 
Pe Ditbe— V(b + A*) + AP V/ (E27 + A*)}-3S,T. . . . (4.38) 


The summation can be rewritten in the form 


Dae 2 Sut 2 {bul —V (Ex? + A) + $A2(A2 + é,2)-19}, (4.84) 
k<O 


ee eee ee 
} Here we denote the entropy in the superconducting state by the index s. 
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The first sum in (4.34) is independent of the energy gap and since it is a 
constant common to both the normal and superconducting states, it may 
be omitted. The second sum may be written in the form of a convergent 
integral, the magnitude of which is 


—}vA?, 
Thus we finally get 
PavAs Sore ek ek el}. (4,35) 
and since the free energy in the normal state is given by 
eee wd ee 2 ks he eae ud cee - 4.06) 
we find, using (4.29), 
2/8 =4fvA2?+4(S,—S,)T. . . . . . (4,87) 


This formula relates the critical magnetic field to the energy gap and the 
entropy difference between the superconducting and the normal states. 
Hence at T'=0 we obtain a relation between the energy level separation A, 
and the critical field H,,: 


Ld ta homp ye te te se (4.38) 
For T near to zero S, is exponentially small, while S,, is given by 
S,=4rvT 
So for the temperature variation of H, in the neighbourhood of 7'=0 we 


find 
inline (Wey Part A) cee TAT) x9. (4.39) 


The temperature variation of H, near 7’, is found by substituting (4.26) 
in (4.37). The quadratic terms in A cancel and taking the first non- 
vanishing term we get 


H2/8r=qvAt/16T 2, . . . . . . (4.40) 
and using (4.22) and (4.38) this gives 
ipa 3 (lett Pe) eae. Siu. so (4:41) 


as the temperature variation of H, near 7',. 

This temperature dependence of H, near 7’, is in agreement with the 
experimental data (Shoenberg 1955) which are usually described over the 
whole temperature region by the formula 


He Be ee 8) ee raed ee ee (44) 


It is evident that these two formulae do not give widely different values for 
H, when 7 is near 7’, and for temperatures close to 7’=0 there is also 
good agreement between theory and experiment, as can be seen by 
comparing formulae (4.41) and (4.42). 

To conclude this section we give some relations which are obtained by 
differentiation of (4.37): 


Se (Cte | SOA OT, cegee: MMe: eae ree ea EAS) 
H2/80=(v/4T) (TA2)(@T+47(C,—Cp). . » + (4.44) 
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The relation (4.43) is noteworthy in that it enables us to express the 
entropy in the superconducting state as a certain integral of the derivative 
dA/d7'; thus writing Cin (4.43) as 7 0S,/07 and integrating we find 


ak ee 
= peopel / Up IT As es 9 
De a 7 dT dT ( ) 
4.4. The Number of ‘Normal’ Electrons 


The so-called ‘number of normal electrons’ is important in the descrip- 
tion of the behaviour of superconductors in magnetic fields in certain cases 
(in the so-called London region). This quantity may be defined as a 
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coefficient of proportionality between the total momentum of the excita- 
tions and their macroscopic velocity. The total excitation momentum is 
equal to 

P= 2k (ag Ny) . . . . . . . (4.45) 


In the presence of macroscopic motion the argument of the function n, 
will not be «(k), but ¢«(k)—ku. Treating u as small, and expanding 
n(e—ku) in a series in ku, we find 


2 On). ky? o 9 
P=— = dk(k Bes ae da 
Tt (ku) 5) aye Per: . . (4.46) 
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It is now simple to find the ratio of the number of normal electrons in the 
superconducting state of the total number of ‘free’ electrons 


wenn f sae. UW Alor as (207) 

a the new variable = A cosh ¢, this formula can be rewritten as 
cosh ¢ dd 

xo aa (exp [A cosh ¢/7']+ 1)(exp [— Acosh d/7]+1)° * 3?) 


Differentiating (3.16) with respect to the temperature, we note that V,,/N, 
may be expressed in terms of dA/07' as 


Nena OAG hoe (A | 
a par] arr) Re een 474.9) 


Hence we find the following simple relation for the number of ‘super- 
conducting’ electrons, V,=N,—N,, where 


For 7'+0 the limiting form for V,,/N is 
N,/No=(27A,/T)¥2exp(—A,/T), . . . . (451) 
while for 7’>7', we have 
Dee Ngee ee idee kw ee ws | (4.52) 


As for the specific heat the next term in the expansion in powers of 7',— 7 
will be proportional to (Z',—7')?. The full curve of the temperature 
variation of V,,(7')/Ng is shown in fig. 4. 


§ 5. Tue ELECTRODYNAMIC EQUATIONS OF SUPERCONDUCTORS 
5.1. The Absolute Zero (T' =0) 


We consider a superconductor in a constant (time-independent) 
magnetic field, the vector potential of which is 


A(r)=>A,exp(iqr). . . . . . . (6.1) 


and calculate the current caused by this field. Todo this we find the change 
5H in the energy of the system due to the presence of the field, and then use 


the relation 
6H = — joA/c AM sig rite any eee Oe) 


to find the current. 
The Hamiltonian for the system in the presence of the field is 


1 y € : e? 
Hy=H+ 2 arth, o)alt’e){—Sterk A+ Gael. . (63) 


where H is the unperturbed Hamiltonian of the problem (3.1). Next 
we transform from the particle operators a(k, c) and a(k, a) to the excitation 
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operators a(k,1) and a(k,0) using (3.3), and we find the transformed 
Hamiltonian H,;=H +H’, with 
ae | —S(kK+k A+ =A | [1 ety, {ort (He, 1)(’, 2) + oc*(k, 0) « (t, 0)} 

+0; V4 {a(—k, 0) at(—k’, 0) +a(—k, 1) a+(—k’, 1)} 

+ Uy, Vqfat(k, 1) a+(—k’, 0) —at(k, 0) at(—k', 1)} 

+0; Up {a( —k, 0) a(h’, 1) —a( —k, 1) a(k’, 0)} J. «4S he onal ene ree 
It is evident that the change in energy of the system will be only second 
order in A since A has no diagonal matrix elements. To find the change in 
energy it is necessary to take the first approximation of perturbation theory 
for the quadratic terms in A, and the second approximation for the linear 
terms in A which appear in H’. At 7’=0 the non-vanishing diagonal 
matrix elements contain only terms with the operator product a«*. Thus, 
taking the diagonal matrix elements of H’, we find the energy change in 
first-order perturbation theory 


? e2 = 
SH;=) A’ ,,= ne aan é haut ee eT 
2 


In the calculation of the second approximation of perturbation theory there 
are possible transitions in which two excitations with opposite spins are 
first created, and then annihilated. The matrix element of H’ correspon- 
ding to the creation of a pair of excitations in accordance with (5.4) is 


(24, Vig — Uy, Vp—g)(2k ra 3 Aq), 
and the energy difference between the ground state and the state with 
two excitations is equal to 


— (ext €xq)- 
Using this, we find for the second approximation of perturbation theory 
5E ae es e >> (UP p—g — VEU Rg) (2k — q, A,)(2k — q, A_,) (5 6) 
oe 4m?c* Ex + €x—g : ; 


and using (5.2), and varying the sum 6#,+6#,,;, we find the Fourier 
component of the current 

jg= ee { — ¥20,7A,+ is 2) (2k — q)(2k — 4, Ag) (Yea) oe : (5.7) 
2m * Ext €x—g 

Thus we have a relation between the current j and the vector potential A 
which shows clearly that there is in general no local relation of the type of 
the London equations between the current j and the potential A. The 
expression derived for j, has one main shortcoming, namely that it is not 
gauge-invariant, as can be seen if we calculate div j, which according to the 
continuity condition should be zero. In terms of the Fourier components 
this condition becomes qj, =0, and it is easy to see that in the general case 
(4.7) does not satisfy this condition. This point, which was not noticed by 
BOS (1957), is not surprising, since the technique of perturbation theory 
which was applied is itself not gauge-invariant. Thus, in (5.7) we must 
understand by A, only the transverse part of the potential. In fact the 
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transformation of (3.3) from the operators a(k,c) to the pair of operators 

a(k,1) and «(—k, 1) is carried out in such a way that the pair of bound 
electrons created always has zero momentum, and this last property of the 
theory is evidently not invariant since the pairs may move in the field. 
However, a more detailed analysis of this question shows that only the 
transverse component of the vector potential, which satisfies the condition 
div A=0, makes a contribution to the current. The transverse nature of 
the first term in j,, which arises from SH,, is evident, and as regards the 
second term, after multiplication by q and transformation from the 
variable k to the variable k —$q this becomes 


20% Ss (Kq)(KAQ)(Y,— 39M e440 — Un—sa b-+40) 

mc?“ Entiat €e—4g 
It is easy to see that as a result of integration over the angles of the vector k 
this expression will contain the vector potential only in the product qA, and 
consequently, if (qA,)=0, then j,.q=0 also. 

The expression (5.7) for the current has a property which is important in 
the further discussion. This is that for large values of the wave vector q the 
energy gap is not of importance so that in this energy region the electrons 
must behave like free electrons. Thus, since free electrons are not super- 
conducting, the current j, must tend to zero as the wave vector q tends to 
infinity. We will show in fact that for large values of q the first and second 
terms in (5.7) cancel out. To do this we use the symmetrical form for 
the second term in the curly brackets in (5.7) : 


. (5.8) 


suede 53 EKA )(Op—19 +40 — Vet yon—4a) 

m €x—tq T Ex+h9 

We now express v and w in terms of the eer, 
E,, €=v(|k + 3q|—h 
making use of (4.5). For large values of q, 
E,,€ =&+4u,c08 0. 

0 is here the angle between the vectors k and q; q is large but |q|<). We 
find for these large values 


2 ss ee ea 
E> FEMA ey Hl ete vee) 


In this sum the only non-zero terms are those in which €, and ¢_ have 
opposite signs. Transforming the sum to an integral, we have 
2 

yo= © 4 (0088) | a Gann RMA = NA . (5.10) 
the limits of the second integral being ky +4qcos@. Here N is the total 
number of electrons (per cm*), given by 

Niza Sak t/3( 2m) ee ees aay A eA EL) 

The first term in the curly brackets in the expression (5.7) for j, is evidently 
equal to — VA, and therefore for large q both terms in j, cancel out as stated 
above. 
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In future the relation between the current and the transverse component 
of the potential A, will be written in the form 
dg=cK(q)A,/40 . . « . + 8 = 40.12) 
with a kernel K(q) which by (5.7) and (5.8) is equal to 
Ne* 1 ko" | sin? 8 cos? 6 d(cos 0) dd dé ( LaSaeae =) 
me? (277) Nmv €,+e€_ €,€_ 
(5.13) 
In this expression we may integrate with respect to dé from — 00 to + co and 


over the angle 4, and we get, setting ,.=cos 0, an expression for the kernel 
K(q) at 7=0, 


K(q)= 


(5.14) 


Rie 4nNe? 3A, |i sinh (quv/2A,)(1— pn?) du 

me 2\glo) 1 V+ eur/4h9 
which as stated above tends to zero for sufficiently large values of g. The 
field distribution in the superconductor (Meissner effect) may be found if 
we apply to equation (5.14) the Maxwell equation (div A=0) 


Me —4ifin 8 eee 


Equation (5.12) gives a relation between the Fourier components of the 
current j, and the field potential A, for arbitrary values of the wave vector q. 
In the general case the relation derived is non-local, i.e. the coefficient 
K(q) for A, depends on the magnitude of the vector q, and in coordinate 
space this corresponds to some integral relation between j(r) and A(r) 
of the form 


j(r)= | K (r—r’)A(r’) dr’. 


Thus the value of the current at any point is determined by an integral 
of the values of the field potential in the neighbourhood of the point. We 
call such a relation non-local. In the old London theory the function 
K(r—r’) had the character of a 5-function (A(q) =const), so that j(r) was 
proportional to A(r), we shall but see that this proportionality holds only 
for small values of the wave vector. The non-local character of the 
relation between j and A was predicted by Pippard (1953). 


5.1. The Penetration Depth 

We shall now apply eqns. (5.14) and (5.15) to a semi-infinite region, 
i.e. we shall consider a semi-infinite metal. It is natural to question the 
validity of (5.14) in this case, since the equation was derived for an 
infinite medium. Fortunately, however, the situation is favourable as 
regards this question and the presence of the boundary turns out not to 
affect the kernel K(q). The fact is that the wavelength of the excitations is 

small and is of order of magnitude 
A=1/K, Ppt na sais) ius gen eg 
i.e. of the interatomic distances. Thus the presence of the boundary will 
affect only those excitations within a narrow layer at the boundary, of depth 
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of the order of the interatomic distances. Since, the penetration depth 
is much greater than this, it is natural that the expression (5. 12) with K(q) 
defined by the relation (5.14) should be valid in all cases, i.e. for any spatial 
configuration. 

The second question which arises in the presence of a boundary is con- 
nected with the character of the reflection of the excitations at the boundary. 
It is known, that in some problems, for example, in the anomalous skin 
effect (Reuter and Sondheimer 1948), practically the same results are 
obtained in the two limiting cases of specular and diffuse reflection. We 
may hope that in the present case too the character of the reflection of 
the excitations at the boundary does not materially affect the penetration 
depth. We shall therefore first discuss the case of specular reflection since 
it is the more convenient in connection with this problem. The presence of 
a specularly reflecting boundary enables us to reflect the potential at the 
boundary, i.e. A(z)=A(—z) (z is in the direction of the normal to the 
boundary), and in terms of the Fourier components of the potential this 
condition becomes : 

Ape eee WN. see!) ans” 6 (G17) 


If there is a magnetic field at the boundary of the metal it follows that the 
derivative 0A/dz will not be continuous across the boundary z=0. Thusif 
0A/0z were continuous, then for an even function at the boundary (z= 0) the 
derivative and consequently the field also would be zero. Since the de- 
rivatives 0A/dz at z= +0 and z= —0 are not equal, the second derivative 
0?A /dz? at z=0 contains a 5-function of the coordinate z. Therefore for a 
half-space we must add on the right-hand side of eqn. (5.15) a term cé(z) 
with a coefficient c, the magnitude of which can be easily determined from 
the boundary conditions. In fact the derivative 0A/dzatz=0is equal tothe 
value of the magnetic field at the boundary, so that the second derivative, 
must be 2/,8(z), and it follows that c=2M,. In the case of a half-space 
all the quantities depend only on the z coordinate, so rewriting (5.15) in 
Fourier components and taking into account the presence of the 6-function 
on the right-hand side, we get 


PA =4nj jet Holm <i --- + + » (5.18) 


Substituting into this the expression (5.12) for j, and solving the resultant 
equation for A,, we find 


A, =Hylmg+ Kg}. 9. - -2- . -, (6-19) 
or, in the coordinate representation 
2H, ( ° cosqzdq : 
A@)= = ped bles wee ears (nt) 
OT Jo PHEW) 


Differentiating this with respect to z and taking into account the boundary 
condition H =H,, for z=0, we find an expression for the magnetic field in 
the semi-infinite metal 


te 2H, (° singz qdq 5.21 
H(@)=Hy- =| +k - - . . ( . ) 


P.M.S. F 


66 I. M. Khalatnikov and A. A. Abrikosov on the 


Finally we calculate penetration depth 5, which is defined by the relation 

seal 0 de ee 

°F II tae h, =i ¢+K(q) wae, 

and we shall now analyse this expression for the penetration depth in more 

detail. First we examine the form of the function K(q) for the limits of 

small and large values of g, and find from (5.14) after some fairly simple 
calculations, for the two cases 


47Ne? 1 gv 5.23.) 
K(q)= 22 (3 30 aa) Re etsy 
_ 4n Ne 3x* Ao 1 ee) | (5.23) 
mce* 4qv nm quv/Ao 
The two expressions for K(q) are approximately equal for 
Yo ¥ 377A, / 40. 


When the penetration depth 6 is small compared to 1/q) (Pippard case), 
large values of g are important in the integral (5.22), and using (5.23 b) we 


get 
sa 2 4nNe? 3n®Ag\-¥8 (° dq _ 4 (3n*Ne* ay (5.24) 
a\ me 4 v o @+1 3/3\ me v 


In the other limiting case (London case), when 6<1/q¢) we use (5.23 a) and 
obtain 


S=(4nNe2{/me?)2?, =. 2... (5.25) 
which is the well-known formula of the London theory. It should be noted 
that for NV ~ 107? and A, ~5°K both formulae give approximately the same 
numerical results for ihe penetration depth. It turns out that supercon- 
ductors of both the Pippard and London types occur in practice; thus 
aluminium and mercury are Pippard metals, while lead and indium are 
London metals and there are also intermediate cases such as tin. 


Dial en U 


We now consider a superconductor in a magnetic field at a non-zero 
temperature. As before the magnetic field is represented by the vector 
potential A, and:as for the case when 7'=0 we calculate the change in the 
energy, using the Hamiltonian (5.4). However, the contribution to the 
second approximation of perturbation theory will now give not only terms 
containing the two operators a+ or a, but also terms containing the product 
ata. Inthe present case the product at(k)a(k’), which corresponds to the 
creation of an excitation in the state k and the annihilation of an excitation 
in the state k’, can give rise to the necessary matrix elements in the trans- 
formation since the system is not in the ground state. Here the action of 


the operators a*(k) and a(k) on the ys-function of the unperturbed state is 
defined by the following relations 


at(k) yl.» <%yeu any Leam,) i Cane eee 
(5.26). 
ak) f(... Mp...) =a/ mg (...m,+1...). 
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In the first approximation of perturbation theory the change in energy is 
e 
0b =>’ .= nee A, A_4[vj2+ (uj2—¥,2)n,]. . (5.27) 


Since, however, the difference wu, —y ; by (4.5) is an odd function of €, the 


sum >(u,?—v,?)n;, is zero. Hence for 5H, we get the same formula as 
before (see (5.5)) 


Oi (ene) AA ha 4 sy. (5.28) 
The calculation of the energy change 6#,, in the second approximation 


of perturbation theory is also easily carried out. The matrix element of H’, 
corresponding to the creation of pairs of excitations is, by (5.4) and (5.26) 


(Vy, Ugg — Ux Ve —g)(2k—G, Ag) V(1—nz_9) V (1-1). 
The matrix element corresponding to the annihilation of pairs in the 
unperturbed state is (u,0,_¢—UxUx—g)(2k-G@, Ax) V(1—"z_¢) Vy, and 
finally the matrix element corresponding to the simultaneous creation and 
annihilation of excitations is 
(Wy Up—g t VE VE—g) V(1— My) /MK-g- 

We may now without difficulty write down the second approximation to 
the energy change as 


: —V,U,_,)2(2k—q, A,)(2k —q, A_ 
Vise e€ > Pima Vby—g)"( q, A,)( q a 


4m?2c? Gas €x—q 
(Uy Urq am VE Vp—q) (2k a q,A,)(2k Tia q,A_,)( Ny, ue Me—a) 


Oa Ex—g 


+(1—1,—My_9) + 
(5.29) 


Transforming from the variable k to the variable k—3q and taking into 
account the expression (4.8) for n,, we get a symmetrical formula for 6#: 


2 
bf = 68, + dBi = = | 32m beA 


a 3 (Un— sq Veda — Vk—30 UR +40)" 


uh €xp—49 1 Ede 


_ (Un-sq Uetagt Vn—40 Un+sq)® (kA,)(KA,) 
€x+ta— ©k—44 
27 


Varying this with respect to A_,, we find an expression for the Fourier 
component of the current 


x {tanh Sk —tanh EH} |. ae ee eS 


Jaa On OAL tee eet) 
Next we use the expressions (4.5) for u, and v;, and take into account the 


PMS. G 
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fact that &, $= EF +Q? cos@ and after some fairly simple transformations 
we find, where A, is the transverse component of the potential 


K(q)= aaa Ay - 1+ raz [sintodeose | ag { ae} 
x (sate (1- aan (tonh $5 a +tanh =) 


é, +e €,€- oT 27 
L £,€ + A* ene 5.32 
+o (14 BES tanh — oT tanh = . (5.32) 


in which the argument k + }q is denoted by the index +, and the argument 
k —3q by the index —. 

As for 7'=0 the relation (5.32) obtained between the vector potential and 
the Fourier component of the current is of non-local character, the value of 
j(r) at any point being determined by the values of the field potential A(r) 
over a certain region around this point. We now transform the expression 
obtained for K(q); first we collect the terms with tanh ¢, and tanh e_ in 
(5.32), and then integrate with respect to dé and then integrate a part of the 
terms with respect to dé. As a result some of the terms cancel and we find 


Cie 47Ne? (Ff lee sin? 6d cos @ 
me? \ 8qv ~ €cosd 
fh eS _ tanh (e ee) dé. (6.38) 
€y es 

The integration of (5.33) cannot be carried further in the general case, 
so we now consider only the limiting cases of large and small values of gq. 
For large values of g the region of small values (near to zero) of € and 
1.=cos6@ are important, and we may therefore put € and yu equal to zero 
in the slowly varying function tanh (e/27') in the integral (5.33). We 
then get 


ro) 5 (8) (f° Daan), . a 


and further integration may be carried out in terms of elementary functions. 
After a change of variable from & to 1/7, the integration with respect to 
dy iselementary. In the resulting integral with respect to dy we make the 
substitution qvu/2A =sinh y, and find 

127 A x dy A 


h—. Py comet (0, 
mez qu J o Seer tanh 7 Soa 


K(q)= — 


Finally, carrying out this simple integration, we get 


K(q)= — — 
K(q) ae riety ~SGres Nevill OOOr 


This formula differs from (5.23) only in the replacement of A, by 


A tanh (A/27'). 
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We now obtain a formula for the limiting case of small q; this can be 
obtained most easily from the initial eqn. (5.32), which gives 


Raye ar {i ; | Ss (‘anh 5§ ar) a} .. (8.37) 


But }$(0/d¢) tanh (¢/27') is simply —d@n/d« and by (4.47) —f(dn/dc)dé is 
equal to the ratio of the number of normal electrons to the total number of 
electrons. Thus for g>0 we have 


K(q) =4nN ,e?/mc?, SO Aes ry eA PD toN 
where the number of superconducting electrons NV, is given by 
N,=Ni(l+[(enjde)dé). . . . . . . (5.39) 


Near 7’=T,, N, changes linearly with temperature, and by (4.52) and 
(5.38) K(q) is in this region given by 


4nNe? (7 .—-T 
K(qg= 2—£ 
NE err ( T. ) (5.40) 
while, for large values of g, K(q) by (4.22) and (5.36) is 
,_ sNer A Dead 
K(g)= Oe : : 
WD a { } (5.41) 
The value qo, for which the two formulae coincide is 
Gea LUNG, Warmest iri! tee coy \(Oee) 


while for q) at 7’=0 the corresponding value of q is from (5.23) approxi- 
mately 7-5A,/v. Thus we arrive at the conclusion, which as we shall see 
later is of fundamental importance, that the magnitude of the limiting 
wave vector q), separating the regions of large and small values of q, 
approximately is independent of temperature. 

The magnitude of the penetration depth at T=0 is as at 7'= 0 determined 
by (5.22). 

When the penetration depth is small in the sense 6<1/q), we get from 


(5.36) 
4 (308N2@A, , A) A A) 78 
vate aa ean ten © 4(b.43 

hs ORE me = tanh sr} Broom. Moor Co 


In the opposite limiting case > 1/qo, using (5.38), we get the well-known 
formulae of the London theory 


ag jel Pe yiuSet seit 4:8 (5144) 


mc 
Evidently for 7 near to 7’, this last limiting case will always apply, i.e. 
the London formulae for 5 will be valid. The observed temperature 
dependence of the penetration depth is often described by the empirical 


formula i 
ee ae 4irNe*| — 
5=By1{] is rit Eye {=} sts CAO) 
while near 7’, the theoretical formula is of the form 
§=S85,/V{2(1-T/L.)}} «© + et (5.46) 


which is 2-1/2 times as large as the limiting form of (5.45). 
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Since the binding energy of an electron pair is approximately A,, the 
width of the region in momentum space which characterizes a pair is 
A,/v and consequently the linear dimensions of a pair must be of order 
of magnitude v/Ay. In fact the limiting wave vector can be considered as 
providing a more accurate measure of the dimensions of a pair; thus 
corresponding to the wave vector q, the length 


A= lig oiBAy + 2 es ee Oe 


is a characteristic parameter of the theory. 

The old London theory is valid in the case when the penetration depth 
is greater than A). In this case only the wave vectors q < q» are important 
in the kernel K(q), and the relation between j and A becomes local (the 
coefficient of proportionality between j, and A, is independent of q). 
Thus the London theory is valid when the dimensions of a pair are small 
compared to the penetration depth. On the other hand, when the dimen- 
sions of the pair are large compared to the penetration depth (Pippard 
case), the new theory is valid, in which the relation between j and A is 
non-local. It is interesting that the characteristic length ‘the size of a 
pair’ is approximately independent of the temperature. 


5.1.2. Diffuse scattering 


Up to now, in the calculations of the penetration depth, we have assumed 
specular reflection in the scattering of electrons from the surface of a metal. 
We now consider the other limiting case in which the electrons are reflected 
diffusely ; this can be easily calculated if we use the results of the paper by 
Reuter and Sondheimer (1948) on the anomalous skin effect. If we were 
to rewrite (5.18), which defines the potential A in the coordinate representa- 
tion, then, as we have already said, it would become 


fa 


AA(r) = — | K(r—r')A(r’)dr’ . 2... (5.48) 


with a kernel K(r—r’), the Fourier component of which is K(q). For 
specular reflection the region of integration is from — 0 to o (A in this 
case is reflected in the second half of space). For diffuse scattering, 
however, the region of integration evidently extends from 0 to oo; and the 
Fourier method is no longer applicable in the solution of the equation, 
and it is necessary to use the Wiener-Hopf method. We use the results 
of Reuter and Sondheimer (1948) directly and write the formula for the 
penetration depth in the diffuse case 


2 i 
saa/ gm {1+ 5K}. Rete wn Ae) 


In the London region K(q) is given by formula (5.38) and on integration 
by parts we get 


8={40N e2/me%}-2, |. |... (6,50) 
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In the other limiting case (Pippard region) K(qg) is given by (5.36). 
Substituting for K(q) in (5.51) and integrating by parts, we get 

_ /3. / (37 Ne A Nees 

=*+ / — = tanh ah Me ac ( 5151) 
As was expected the penetration depth in the London region does not 
depend on the type of reflection. In the Pippard region (small penetration 
depths) we get, by comparing (5.43) and (5.51), the relation 

OU See) Glee, als ane we (B52) 

so that even in this case the dependence on the type of reflection is small. 


5.1.3. Remarks on the properties of superconductors of finite size 


In the case just considered, of a superconductor filling half of space, 
the vector potential A automatically satisfies the condition div A= 0, and, 
as followed from (5.7), the condition div j=0, necessary for solutions in 
stationary conditions is satisfied. This point needs some discussion for 
finite conductors of arbitrary form, where the condition div A=0 is indeed 
sufficient, but not necessary for the fulfilment of the conservation condition 
div j=0. This in the general case of superconductors of arbitrary shape 
it is necessary to solve the complete system of equations, consisting of the 
Maxwell equation and eqn. (5.32) for the superconducting current 


curl curl A= — 47j/c 
| dq=cK(q)A,/4z. 


In agreement with (5.53) we have at the conservation equation for the 
current 


(5.53) 


divi j= 0; J PS Pee Sar ee i 9) 


In the general case the vector potential for this will no longer satisfy the 
condition div A=0. 

As an example of the problem of a superconductor of finite size we will 
calculate the magnetic moment in a constant magnetic field Ho, of a 
superconducting sphere of radius compared to the penetration depth. 
For the London region the solution, as is well known, is 


MSO | S005 oe Pees ee. (0.00) 
where a is the radius of the sphere. 


We shall now solve this problem for the Pippard region. In the 
coordinate representation (5.32) may be written as (Bardeen et al. 1957) 


d(r)=B, [eae Pee i) 


where R= |r— r’ |, nis a unit vector in the direction of R, and Bs is a constant 
given by 
2 
he. Gee oon) 
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Since the body is small compared to the penetration depth, we may avoid 
solving of the Maxwell equation since the vector potential A(r) differs 
little from its value in the absence of the sphere. In a constant magnetic 
field H, the vector potential is 


A(r)=4[H).r]. :. ay ak en tend 
For this special case n(A(r’) —A(r))=0 and therefore A(r’) in (5.57) may 
be replaced by A(r) and 
3 
j(r)=B » | sel) jae Mee tds et) 


The conservation law div j = 0 can also ro written in the form of a boundary 
condition for the normal component of the current at the surface of the 
body, 

74330). Sy ee eo 
and it is evident by symmetry that this condition will be automatically 
satisfied for a sphere. 

We now use the following simple transformation 


nie) = 4 {aaiv 3 (avy 3h bet) See 


and in (5.59) transform to an integration over the surface. Thus we then 
find 


i(r)= 48, {Ale ny ~ {(Aansh. . - . (6.62) 


The two integrals in (5.62) may be evaluated by introducing the angles y 
between r and dr, 0, the polar angle with respect to r as polar axis, and 4, 
the azimuthal angle. In the second integral, since it is evident that j 
will have the same direction as A, only the component of n in the A 
direction need be taken, i.e. sin(@+<2) sind, so that 


j=4B,A(r 4 [| $e0s xaleos 8) ap 


2a 772 
-{ [3 © sin ( (0+ yx) sintsin®g dag aon (iia Oent 


Making use of the relations 
cos 6 = (r?+ a? — R?)/2ar, cos y=(R?+a?—r?)/2aR 
sin (9+ y)=asin 6/R, (5.64) 
the integrals reduce after integration over ¢ to simple integrals over R 


between the limits a—r and a+r. Finally we find, on substituting for A 
from (5.58), 


(5.65) 


at 7B, ane aaah + 378) ate 3ar = 


4 r 4r? a—r. 2. oF 
The calculation of the magnetic moment of the sphere 


M=4 lirlav Wet) Rae lye 2 gle BR) 
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is now elementary. We substitute (5.64) for j in (5.65) and after a simple 
integration find 


Mizar Hondas BalOa) ie? ei et) th (8.67) 
Making use of (5.51) and (5.57), we have that 
B, = 3/3/3273 5p, 
and so finally, for the Pippard case and for a<6,, we havet 
ase Hen 3a°|96a0' pe... i, anes = «* (6.68) 
While in the London case the magnetic moment of the sphere is proportional 


to the square of a/6,, in the Pippard case it is proportional to the cube 
of the ratio a/6,. 


§ 6. THE BEHAVIOUR OF SUPERCONDUCTORS IN A PERIODIC 
(HicH-FREQUENCY) Fretp (Landau 1956) 


The behaviour of superconductors in a high-frequency field can be 
calculated in the same way as for a constant field case. Thus considering 
a superconductor in a magnetic field which varies with frequency w, we 
write the vector potential of the field as 

A(r't)=) A,exp(iqr—tt). <-. . . . (6:1) 
We next find the current operator; from (5.4) we have 
y=(—e/2mc) > (k+k' —2eA/c){ujuy[at(k, lo(h’, 1)+at(h, O)a(k’, 0)] 
+00 ,[a(—k, O)at(—k’, 0)+a(—k, 1)at(—k’, 1)] 
+0, [at(k, lat(—k’, 0)—at(k, O)at(—k’, 1)] 
Oth lat— Kk; Ojala l)=a(ie, Lo(k',.0)]}. e248 oe (6.2) 
The matrix element of the current operator may be found according to 


the standard formulae of perturbation theory for a time dependent 
perturbation (Landau and Lifshitz 1948). We thus find 


J e 
LE = Cees k(kA,) > {(a-paPesta~Pi-salies ta) 


1 1 
it ( pean SiS PR Ie Sie Ce oe 
Exthqt €e-tg bt €x+tgt €x-a— 


x (ton a +tanh ara +a similar term}. coe (0.0) 
Using now the expression (4.5) we finally get 
jqg= cK (q, w)A,/4a, & ieee La ae 
47Ne? 1 €éé 4+? on eo 


A eRAa 
x emigets 4 —) + (1+ oS) 
é, + elPo lene — wa emer 


Gee JSi, 1 een eslinres ) |} ou (02d 
x (tanh 5a tanh) (— Seer) (6.5) 


NN 
+ This result was first obtained by Whitehead (1956) ; note that dp is the 
same as Pippard’s (and Whitehead’s) Ad. 
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Further we take into account the fact that (qA,)=0 and calculate from 
K(q,w) the value K(q,0) for g=0, which is given by formula (5.32). We 
note also that the function K(q, w) is complex, the imaginary part of this 
function may be found from the damping condition, which is taken into 
account in, (6.5) by a small positive imaginary increment 76 to w. We 
then obtain 


K(q, ») = K(q,#)—K(q,°) 
A2 
= ee | sind cos | dé | (1 a 
4mc* J €,€_ 
x {tanh (e,/27') + tanh (e_/27)} 
x {(e,te_+w+id)1+4 (e, +e —w—t6)4—2e, + ej} 
2 
a0 (: + Ht) {tanh (¢,/27') — tanh (¢_/27')} 
: 


+€— 


x {(e,-—€_+w+16)1+4 (€,—€_—w—18)1— 2(€,— <3 | : 


(6.6) 
and using the relation 

1/(a +76) = 1/4 +28(2), 2h. oteleed ae eee ae 
we can find the real and imaginary parts of K(q,w). If we introduce 


the following notation 
K(q, w) = (8? NV A/4emvq)Q(q, w) Seed de Geel 


and transform from the variables € and @ to the variables «,=«,/A, 
€,=e¢_/A, we find for large values of q¢ 


Q(w)—Q(0)=4 | (2-1) de, |” (€,2— 1)" de, 


1 1 

x [(1 — €,€,){tanh (Ae,/27') + tanh (Ac,/27')} 

x {(ey + eg + @— 15) 14 (e, + €.-—@—8)4 

— 2(€, + €,)}— (1+ e,e,){tanh (Ae,/27') — tanh (Ae,/27)} 

x {(€y— €g + @ + 18) 1 + (Ey -—€g —G—18)4 

— 2(€,— €9)}] « . O52 pee nye aes 
where &=w/A. 


ioe) 


6.1. The Pippard Limiting Case 
It is not possible to carry out further calculations without more concrete 
assumptions as to the values of the important parameters of the problem. 
Here the most important circumstance is that apparently for most super- 
conductors the penetration depth is much less than v/A, ie. vg> A. 
This inequality becomes invalid only in the neighbourhood of the transition 
point 7’, where there is a region for which »g<A. In the main region the 


} Strictly “speaking, the strong inequality vg>A is not valid. The applic- 
ability of the formulae which are obtained in this limiting case is however 
guaranteed by the presence of considerable numerical factors. 
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electrodynamics are based on the equation suggested by Pippard, while 
in the second region (near 7’,) the London equations are applicable 
(London region)+. 

For the case of an alternating field the London region becomes narrower. 
This is clear if only from the fact that the normal metal must always be 
of the Pippard type, since in it A=0, and so it follows that the London 
region cannot continue right up to the transition point. Leaving aside 
for the present a detailed discussion of this question, we shall now consider 
how (6.6) may be simplified in the Pippard region. 

In (5.9) we make the change of variables ¢,=«,, so that all the terms 
become proportional to tanh(Ae,/27’). Collecting terms, we arrive at 
the following expression : 


ma Ae de Foe Lethe 
w)—Q(0) = tanh 5a) ie Ua! ae 0 | pee See 
O(e)—@t0)= | Pande) eran, Weed 
x ete to) +1 eet )et 2(€,?+ 1) ; (6.10) 
(ey +@+76)?—€.2 (e,—-G—16)?—€.? =e? -€,? 
This expression can easily be integrated with respect to «,. Supposing for 
definiteness that w > 0, and separating real and imaginary parts, we havet 


Layee ete cae el 0) atest 
Q(o)=7 I eyt™ 27 »/[(e?— 1){1 — e(e—&)}] 
amon abate yin ee eis hy ek 
— in6(& — 2) ie tanh oP /e—ie—eP— 1} 
? -) A = Ae e(wte) +1 € 
Bae al eam 


where 6(~)=1(a>0) or O(x<0), and G=a/A. 

First of all we consider the case 7’=0. In this case (6.11) simplifies 
considerably and it is easy to see that it can be expressed in terms of the 
elliptic integrals: 


2n7H(4H) @B<2 
Q(6)= 4 w[GH(2/6) —(@?-—4)@1K(2/w)—i17O(G@—2)]  . . « (6.12) 
x [@H(V 1 — (2/)?) —4@2K (V1 — 2/6?)] @ > 2. 
For low frequencies the limiting value of this expression is 
OG) al (2/4) yi: et a eto ee Ord) 
For > 1 we have 
Q(@) = 77/4 —in{@? —2G(In2G+4$)}.  . . . . (6.14) 


+ Apparently the case for which the superconductor is described by the 
equations of F. and H. London over the whole temperature region at w=0 
is much rarer, and we shall not consider it. : eee 

t In our paper (Abrikosov et al. (1958)) there was an error In the derivation 
of this expression, and this also affected the equations which followed. The 
correct equation appears in the recently published work of Bardeen and 


Mattis (1958). 
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We also quote the expression for Q(w) in the neighbourhood of & = 2 


i | 2 1) | — in?6(a® —2)($6—1). (6.15) 
2-—@ 

At 7'=0 the imaginary part of Q(@) is non-zero only for @>2. This is 
quite understandable from physical considerations as the imaginary part 
of the current determines the absorption, and in the absence of excitations 
this absorption can take place only by the break-up of pairs. The quantity 
2A, taking into account the statistics, defines the amount of energy necessary 
for this. 

For 7'40 the integrals depend considerably on the relation between 
A, wand 7’. It should be noted that the condition A<T is only satisfied 
in the immediate neighbourhood of the critical temperature (for example 
A/T <0-1 only for 1—7'/7,,< 0-001), and thus does not correspond to a 
useful temperature interval. Therefore we only consider the cases 
A>T and A~T. 


(a) The case T~w<A 


The imaginary part of Q@(w) is exponentially small, as also is the change 
in the real part with temperature. Transforming the integrals we find 
Q(w, 1) =7[1 — (w/4A)?— 2 exp (— w/27)Iy(w/27) exp (— A/T)] 

— 471 sinh (w/2T)Ky(w/2T)exp(—A/T) . . . . (6.16) 
where A, and J, are Bessel functions ofimaginary argument. In particular 
for T'<w<A we have 

Q(w, 1) =72[1 —(w/4A)?— 2(7/mw) — exp (— A/T)] 
= die] (aT eexp (AT 2 ee be Lake 


and forw<7T'<A 
Q(w, T) =7?[1 — (w/4A)? — 2 exp (— A/T) ]— 2in(w/T) 
xin (47 pulexpt— AIL oh ke. Gos eo eens 
where y= 1-78. 
(b) The casew<T~A 
After a series of transformations we find 


“ A w \2 sh eee oe ae 2A 
O(a, 1) =n° tanh 5741 (4) i —in | Foosh pina, [= 


w A w A 
+ (1-tonh 55) -25P(5) | 5 4 est Beene Cuba 


where the function P(x) is the integral 


Pa =|" de cosh ve — cosh x 6.20 
1 €°—1 (cosh we + 1)(cosh a+ 1) Coe 
whose limiting values are 
forae le P(x) =exp (— 2a) In 2ya | 
7 (6.21) 


 forx<1: P(x) = — &(8)z. | 
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From (6.19) we find that for w<7<A the result is the same as with 
(6.18): forw<A<T: 


Q(o,T)=7 opin | St on (B= — i) |. AGS) 


(c) The case T’<w~A 


Here there are two possibilities. If w<2A the imaginary part is 
exponentially small. The expression for Q(w) is 


Q(w,1)=2n| B(o/2A)- ul {ar (= - sq)exv(- ajay} | 
-2in] [an (= + xa) exp(—A/r) |. 603) 


In particular, we obtain (6.17) for ’<w<A. Ifw>2A,anon-exponen- 
tial term appears in the imaginary part 


Q(w, T) = Q(a, 0) + 27 exp [—(w— A)/T] far (a-3)} 


Berexn (AID) a/ aT) RG 3 sa) - Is 4 ~) | (6.24) 


In particular, for 7’<A<w we have 
Q(w, T) =Q(w, 0) — 2ais/(27T A) exp(—A/T)/w. . . (6.25) 


(d) The case 7~A<w 


The expression for Q(w, 7’) becomes 


4miA [* dd 
ry 0) ee —_____——.._.. 2 
Sere Y) w I. exp{(Acosh ¢)/T} +1 tly) 
Here Q(w,0) denotes the asymptotic value for wsA, (6.14). For 
T<A<w we again obtain (6.25), but for A<7'<w we find 


Q(w, T)=Q(w, 0) —(271Aw) In (wT /yA). . . . (6.27) 


We note that for all cases in which ws A, the real part of Q is small 
compared to the imaginary part, which is —7w/A independently of the 
value of 7. Thus for ws A the relation between current and the vector 
potential is: 

j(q, w) = (37ie?Nw/4mcevg)A(q,w). ©. . . « (6.28) 
This relation does not contain A, and is identical with that obtained from 
the theory of the anomalous skin effect for a normal metal as is to be 
expected. However, in this case we cannot limit ourselves to the basic 
term, since to find the difference between the superconducting and the 
normal states we must also take into account higher order terms. This is 
the reason for the inclusion of a small additional term in the eqns. (6.25), 
(6.26), and (6.27). 
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6.2. The London Region 


We now calculate the conditions for the appearance of the London 
region vg <A, and, as already stated, we shall not consider the possible but 
apparently rarer case when under static conditions the superconductor is of 
the London type over the whole temperature region, 

As has already been noted, the London region can only appear near 7’,. 
On the other hand, however, near to 7’, the energy level separation A 
becomes very small, and since for w/A>1 the metal differs little from the 
normal state, the London region cannot continue right up to 7, and is in 
any case limited at the upper end by the condition w < A. 

To find more accurately the position of the London region we first find 
a relation between the current and the vector potential for ASv,. We 
take into account the previously noted fact that over the whole interesting 
frequency range we may take w <v,, and also that the London region arises 
only very near to 7',, i.e. wemay take A<7. The necessary equation then 
follows from (6.6), in which only the integral with the difference of the two 
tangents and the last term are important. 

Thus we find 


: _  @NA,(w) { 7 A? = 37m 
j(a,0) = - SR) | 08) 7a - Gee} 
If we limit ourselves to the real part of this relation, we get the London 
equation, and the expression in the curly bracket plays the réle of the ratio 
of the number of superconducting electrons NV, to the total number of 
electrons V. The imaginary part corresponds to absorption. We note 
that this expression is obtained under the assumptions w<A<7' and 
w <vq, but there can be an arbitrary relationship between vg, and A and 
between vgand 7’. This explains the fact that apart from a small correction 
form this expression coincides with the corresponding one for the Pippard 
case. 
As a criterion for the transition between the Pippard and London regions 
we may take the equality of the coefficients of A in the expression for the 
current. For the Pippard case, taking w<A<T, from (6.13), (6.15), 


(6.25) we find 
q e2N 372 A? = 3t77w 
=— A - stapes 5. “ 
: me a) { 5 vql, 4uq \ A 
We compare this formula with (6.29). If we assume that the imaginary 
parts of both formulae are small compared to the real, the criterion coincides 
with that for the static case, namely that the London region occurs for 
0 K Ugg = 3774 7',/14¢(3) ~ 187. ty ate tx eet tea 
Taking q~ 1/6, where 6 is the penetration depth, and using formula (6.29), 
and the London formula for the penetration depth 8, = y/{(me?)/42N eh 
it can be shown that this corresponds to the condition 
A <Dar9l8 (me N ahy* NAT ons iis ge eee 0:92) 
From this inequality it follows that, for example, for aluminium the 
London region begins at 7’— 7',~ 4 x 10-4 7, and for tin at 7’ —7',~ 1071 Th 
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As the temperature increases the real part of (6.29) becomes small 
compared to the imaginary. But since in this case the superconductor 
differs little from the normal metal, we necessarily get the Pippard situation 
(we note that the imaginary part of (6.29) refers to this case also). Thus 
we may consider that the London region is limited by the condition 


04 KVG1 = (87rw/40A?) 72, «= (7/477)C(8) 20-2, . . (6.33) 
or, after substitution of the London penetration depth 
Bo me (Aa (meee Tom. ie, (6.34) 


Comparing formulae (6.30) and (6.31) we get the condition for the appear- 
ance of the London region 


@) < w, = 94 /(80:)® (me?/Ne*y?)T 3. . . . « (6.35) 


This condition may define completely different frequencies for different 
elements. For example, for aluminium this frequency corresponds to 
2-7x 10-7", or approximately 3 x 10°sec-!, while for tin it is 0-57’, or 
2x10'sec!. Incidentally tin may be considered by convention to belong 
to the superconductors of the Pippard type, since even at T’=0, vg~ vq. 

In view of the fact that the most interesting frequency range is that near 
to A(0), and that for a Pippard metal this range is practically always 
beyond the limits of the region defined by the inequality (6.35), we will not 
in general consider the London region. 

Here it should also be pointed out that besides the limit on the frequency, 
the applicability of the London equations in the form (6.30) is limited on 
the high-temperature side by the condition that the dimensions of a pair, 
which are of the order v/A, must be small compared to the mean free path. 
If this does not hold, the position changes radically. Since on the average 
the mean free path may be considered to be of order 10-3 cm, and v/A ~ 10~# 
em, the region of applicability of (6.30) is, in general, small. For v/As1 
no electrodynamics has been constructed and we shall not examine this 
case. 

In conclusion we note that, as can easily be seen for frequencies greater 
than w,, a direct estimate gives vg> 7’, over the whole temperature region, 
including the neighbourhood of 7',, where the penetration depth is a 
maximum. Such a relation between vg and 7' was important in an earlier 
derivation. 


6.3. The Impedance 


The relation between j and A can be substituted in Maxwell’s equations 
and an expression for the dependence of the vector potential on the co- 
ordinates can then be obtained. In the Pippard case considered above the 
corresponding calculations do not differ from the static case. or the case 
of specular reflection of the electrons from the surface the vector potential 
must be considered as reflected symmetrically across the boundary. 


80 I. M. Khalatnikov and A. A. Abrikosov on the 


Taking the external field at the boundary as H(0) we get 
H(0) cos kz dk ; (6.36) 
9 k2+(37Ne2A/mev) Q (w)/k 


p= | Hae H(0)= ee (sr) Fe Qo) *91 


The impedance is defined as follows 


Z=R+iX =K(0) / | jdz=47E(0)/cH(0) 


i wee 0 (“ea5) {O(a) ae. ee 


For complex Q(w) the value of the root is determined as an analytic 
continuation of the real root for real Q(w). 
In the diffuse scattering case we must take A(z)=0 outside the metal. 
In accordance with the theory of Reuter and Sondheimer (1948) this gives 
Ania 377%e2N =e 
Z=- dk\n| 1 ‘ : 
ai], S| 2+ aeree Of =< (28 
This expression differs from (6.37) only by the factor 8/9. 
It is convenient to relate the magnitude of the impedance to the value R 
of the resistance in the normal state, which for specular reflection is 
R,, = (87? /34/3c*)(ncton/ BNE ee ee 
as can be shown by substituting Q@ = — 7iw/A in (6.37). The ratio Z/R,, is 
the same for both diffuse scattering and specular reflection and is given by 


w)/R, = —214/{rw/AQ(w)}. . . . . . (6.40) 
For 7'=0 the frequency dependence of the impedance is obtained from 
eqns. (6.12)—(6.15) 


TW . i . 
w)/R, = 2 ~/ ATAw)] {sin (4 tan“ Lu) —7 Cos (a tan Al * 
(Oat) 
It is interesting to note that the penetration aries shanees very little over 
the region w=0 to w= 2A, and that 6(2A)/5(0) =84/ (7/2) =1-16. 


For non-zero temperatures we He now analyse the temperature 
dependence of the impedance at various frequencies. 


where 


(a) w<A(0) 

At the lowest temperatures we have successively 7’¥<w<A, then 7 
becomes of order w, and finally we reach w<7’<A. This transition is 
described by (6.16), (6.17), and (6.18). Substituting these into (6.40) and 
taking into account the fact that Q(w) differs little from the value of Ze Q(w) 
for w=0, 7’=0, we have: 

(w)/R, = 24/ (w/7A){(4/377) sinh (w/27’) Ko(w/27’) exp (— A/T) 
—t[1 + $(w/4A)? + 31, (w/27') exp (— A/T — w/2T')]}. (6.42) 
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When the temperature is increased further we reach w<A~7', which 
then becomes the small region w<A<T. This transition is described by 


(6.19)-(6.22). While A/7'ssw/A, Q(w) will differ little from its value at 
7’ =0, and thus we find 


w)/R, = 24/{w/7A tanh (A/27)} (2/37) In 24/(2A/w) + (1/37)(w/A) 

+ {coth (A/27’) — 1} — (2/377)(w/T’) P(A/T) coth (A/27) 
SURAT Gy EUAN DN MEM col ima Naas a ane aa ee (97 
In the immediate neighbourhood of the critical temperature A becomes 
small and we havew<A<T,w~A<T,andA<w<T. But,as hasalready 
been pointed out, the consideration of these cases is of no interest since they 
are all contained in a very small region. Instead of this we can use the 
fact that although (6.22) is not valid for A<w<T it nevertheless gives the 
correct limiting value, corresponding to the normal metal. Thus this 


expression (but not, of course, (6.42)) can be considered as an interpolation 
formula in the region w> A. 


(b) w~ A(0) 

At low temperature we have 7’<w~A. In this region (6.23) is valid for 
w<2A(0), while (6.24) and (6.25) are valid for w>2A(0). In the first case 
Q(w) will differ little from the value Ze Q(w) for T’=0, which is 27H(w/2A), 
and the expression for the impedance is : 


NY {samo | sme N {"? (5 * 35) POS 
—i{1+ sony af (225 - a) e=P(- ar))t |. SOee(e a2) 


For w > 2A(0) both parts of Q(w) are of comparable magnitude and it is 
simpler to use the general formula (6.41). 

As the temperature is increased we pass through w~A~T' and reach 
A<T~w. Here we can use the fact that for ws A the limiting value of the 
impedance is independent of temperature and corresponds to that for 
the normal metal. Therefore the use of (5.24) for interpolation will not 
introduce any significant error. 


(c) o> A(0) 

In this case only the relation between 7’ and A changes, as w is always 
large compared to them. Equations (6-26)-(6.28) are applicable here, and 
taking into account the fact that for this case the real part of Q(w) is small 
compared to the imaginary part, we find 


AO) <14 (2) {50x +5 | -iva[ a+ ei 


« (Fin ge + 5+ 5) |: te cee ee es Le) 


These equations enable us to compare the theory with the experimental 
data. 
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§ 7. THe THERMAL CONDUCTIVITY OF SUPERCONDUCTORS 


In superconductors heat transfer may take place by electronic excitations 
and by phonons. At temperatures which are not too low the electronic 
excitations make the main contribution to the thermal conductivity. In 
this case there are several possible scattering mechanisms, of which the 
scattering of excitations by impurities will usually be the most important. 
The scattering of excitations by phonons will be important only for very 
pure metals and we shall calculate the thermal conductivity of a super- 
conductor for the case in which impurity scattering predominates. We 
write the kinetic equation for the excitations 


see ae = Ue pe 


Here ¢(n) is the collision integral for the scattering of excitations by im- 
purities, and we first find the form of thisintegral. It contains the scattering 
probability of excitations by impurities and the Hamiltonian of the inter- 
action between an electron and the impurity can be written in the form 


Ht — Yat(k)a(k’) V(k,k’). . . . . « (7.2) 


In this expression we now transform from the particle operators to the 
excitation operators, by (3.3) the terms which describe the scattering 
of excitations will have the form: 


He = ¥ at(k) a(k’) V(k, k’) (w,2—0,2). 2... (7.8) 


and by (4.5) the last term is 
£;,/€p- . . . . . . ° . : (7.4) 


The probability of an elastic transition kk’ for the excitations is 


W5' = Ir |V(k, k’) £,/€% |? p(e,). . ° . . . (7.5) 


The density of states p is found by the integration of the energy 5-function 
with respect to dé: p=e;,/|&,|. Thus the probability of the above 
transition is 


Wey’ = 2 | V(k, I |E, Vex. 2... - (7.8) 


For the normal state the ratio |€,|/e, equals unity, since |£,]=e,. In 
the superconducting state this ratio is some function of k (or €). It is well 
known that in the normal state the scattering of electrons may be charac- 


terized by a time 7 which is independent of temperature, and that we may 
write the collision integral in the form: 


E(t) > — (9 — 1g) reamed tae Pet) 


In the superconducting state we must also take into account the factor 
| é |/e by (7.6), i.e. instead of + in (7.7) we must write 7’ which differs from 7 
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by this factor}. Therefore in future we will write the right-hand side of 
(7.1) in the form 


C(x) = —(n—m9) |E,l/re, ~ - . - : . (7.8) 


where 7 is the time constant which appears in the classical theory of 


electrical conductivity. Thus we write the kinetic equation (7.1) in the 
form 


OcOn  dedn _ _— (n— mM) |E | 


ok or Oorok T € 


In non-superconducting metals an electric field is set up in the presence of a 
temperature gradient such that the resultant electric current is zero. How- 
ever, for a superconductor there are two opposite flows of current—the 
superconducting and the normal, which are equal and opposite. The 
presence of a superconducting current with velocity u may be taken into 
account by substitution in (7.1) 


(7.9) 


e>e+ku. 


We will bear this substitution in mind in further calculations, and we now 
turn to the solution of (7.1). 
We substitute on the left-hand side the equilibrium function 


Vee OXON E/E Lo) ve aye ae es (7410) 
On simplification we get 

Ben AT _ (w—no)le| 

Okde T | TE : 


It is interesting that the derivative Ve cancels in the equation, and that in 
the result there are no terms connected with the differentiation of the 
energy level separation A(7'). The energy flux equals 


(7.11) 


Q= | «(@c/ak)ndry. ST ne. Stake) 


Since the equilibrium function n, depends only on e, the integral (7.12) - 
vanishes for n = n, and the flux Q is completely determined by the difference 
N—N, which is given by (7.11). Substituting for n— 7, from this equation 


£ 


+ The replacement of + by 7’ is carried out, as can easily be seen, by multi- 
plication by the ratio of the velocity of excitations to the velocity of electrons w 
at the Fermi surface. ; 


0€/0En = En/e- 
Consequently we have the relation rv=7’fk/e, which shows that the mean 
free paths for electrons in the normal state and for excitations In the super- 
conducting state are equal. 
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We thus find the magnitude of the coefficient of thermal conductivity to be 


K 


= Th e(o\- eee 7.14 
ok (5) exp(—/T)+lexp(—¢/7)+1 Oe 
If we introduce a new variable of integration ¢=cosh ¢A, then formula 


(6.14) reduces to a form which is convenient for the calculation of the 
temperature dependence. 


Fig. 5 
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378 vA iL cosh p sinh} dt a (een aoa re ave 
Equation (7.15) gives the electronic part of the coefficient of thermal 
conductivity of a superconductor when the electrons are scattered mainly 
at impurities. We note that at the transition point (A— 0) (7.15) gives a 
value of « which corresponds to that for the normal metal obtained from the 
usual theory of the electrical conductivity of metals. It can also be 
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rewritten in the following form which is convenient for comparison with 
experiment 


RAL eel, ip ‘galls bev (7.16) 


KL) aT", J ajo, cosh?” 
Zavaritskii (1958) compared the temperature dependence of « given by 
(7.16) (see fig. 5) with the experimental data for Al and Zn anid obtained good 
agreement. The data for zinc deviate from (7.16) only below 0-2°K, when 
the lattice conductivity becomes important. 


§ 8. CONCLUSION 


Thus, one of the most enigmatic of natural phenomena—superconduc- 
tivity—has at last beenexplained. The agreement ofthe theory of Bardeen, 
Cooper and Schrieffer with experiment has proved even better than could 
have been expected from a theory which is based on the model of a 
Hamiltonian with a four-fermion interaction. At the same time, although 
the phenomenon of superconductivity can be considered to have been 

completely explained, there still remain isolated questions which need to be 
investigated on the basis of this theory. The behaviour of superconductors 
in high magnetic fields (near to H;) has not yet been investigated. In 
connection with this it is apparently necessary to develop a gauge-invariant 
method independent of perturbation theory and the work of Gor’kov (1958), 
in which he has applied the gauge-invariant technique of quantum field 
theory to the theory of superconductivity deserves attention and has 
considerable methodological interest. 

There is also the question of the possibility of the existence of Bose type 
excitations in superconductors (Bogolyubov et al. 1958). The Fermi 
excitations considered above had a characteristic spectrum with an energy 
gap. The presence of, for example, sound vibrations in the form of 
phonons (with an energy spectrum without a gap) would not interfere with 
the existence of superconductivity just as phonons do not prevent super- 
fluidity. However, if there were such a branch in the spectrum, this would 
considerably affect the temperature dependence of all the thermodynamic 
variables. The question of Bose oscillations of the system of Fermi- 
particles considered has however no great importance in the theory of 
superconductivity. This is because such oscillations would be related to 
oscillations of the electrical charge density, which, owing to the large 
Coulomb interaction would be in the optical frequency range, so that to 
excite them, energies of the order of an electron volt would beneeded. Thus 
the whole question of a Bose branch in the spectrum of Fermi particles is 
only of academic interest. 
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